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FOREWORD 


This  document  is  the  final  report  for  the  Hardened  Radar 
Component  Materials  Development  Study.  This  work  was  performed 
by  Science  Applications,  Inc.  (SAI)  for  the  Army  Materials  and 
Mechanics  Research  Center  (AMMRC)  under  Contract  DAAG60-78-C-0052 . 
The  period  of  technical  performance  was  May  28,  1978  through  Jan¬ 
uary  20,  1979.  The  AMMRC  program  manager  was  Mr.  John  Dignam. 

The  SAI  Program  Team  included  the  following  key  participants: 

Mr.  Irving  Osofsky  -  Principal  Investigator,  Mater¬ 

ials  and  Requirements 

Mr.  Thomas  Duncan  -  Thermal  Analysis  and  Test 

Planning 

Dr,  Robert  Semelsberger-  Electrical  Property  Test 

Planning  and  ANalysis 

Dr.  Donald  Dudley  -  Electrical  Property  Test 

Planning 

Mr.  Leonard  Porcello  -  Electrical  Property  Test 

Planning 

Mr.  John  Dishon  -  Thermal  Radiation  Simulator 

Test  Conduct 

Mr.  Kenneth  Sites  -  Instrumentation 

Thermal  testing  was  performed  at  the  French  Solar  Furnace  in 
Odeillo,  France  and  at  the  SAI  Thermal  Radiation  Simulator 
Facility  at  the  Nuclear  test  site  in  Nevada.  The  electrical 
properties  tests  were  conducted  at  the  Georgia  Institute  of 
Technology,  Engineering  Experiment  Station,  with  Dr.  J. 

Fuller  as  principal  Investigator. 
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1.0 


INTRODUCTION  AND  SUMMARY 


The  multiple  nuclear  burst  environment  resulting  from 
an  attack  on  ground  based  US  ICBMs  provides  a  severe  environ¬ 
ment  for  the  nearby  ground-based  phased  array  radars  which  are 
part  of  a  ballistic  missile  defense  system.  The  phased  array 
radars  must  survive  the  nuclear  blast  overpressure  and  result¬ 
ing  thermal  environment  with  minimal  performance  degradation 
prior- to  as  well  as  post-attack. 

The  severity  of  the  nuclear  burst  generated  thermal 
radiation  environment  requires  identification  of  resistant 
materials  which  satisfy  the  thermal,  structural  and  electrical 
requirements  with  minimal  phased  array  radar  performance 
degradation.  The  objective  of  the  present  study  was  to  identify 
and  evaluate  materials  for  nuclear  thermal  radiation  and  blast 
hardening  of  BMD  phased  array  radar  components. 

The  scope  of  the  present  study  included  identification 
of  candidate  materials,  compilation  of  existing  structural  and 
thermal  property  data,  design  and  fabrication  of  test  specimens 
for  thermal,  electrical  and  mechanical  property  tests,  design¬ 
ing,  planning,  and  conducting  thermal  and  electrical  property 
tests,  and  simulating  the  performance  of  candidate  materials 
in  nuclear  environments.  A  wide  range  of  candidate  materials 
for  both  radiating  and  non-radiating  phased  array  radar  compon¬ 
ents  was  identified.  Emphasis  was  placed  on  the  use  of  unin¬ 
filtrated  and  infiltrated  porous  ceramic  coatings  and  low- 
carbonaceous  char  ablative  materials.  However,  the  materials 
list  includes  metals,  bulk  and  composite  ceramics  as  well. 
Limited  thermal  radiation  tests  were  conducted  at  the  French 


Solar  Furnace  Facility  in  France  and  at  the  SAI  Thermal  Radia¬ 
tion  Test  Facility  at  the  Nevada  Nuclear  Test  Site.  Electrical 
property  tests  at  both  S-band  and  X-band  were  conducted  at  the 
Georgia  Institute  of  Technology  Engineering  Experiment  Station 
in  Atlanta,  Georgia. 

The  results  of  this  study  showed  that  uninfiltrated 
ceramic  coatings  and  a  silicone  rubber  ablatiave  material  had 
promise  for  thermal  protection  of  radiating  phased  array  radar 
components.  The  use  of  an  infiltrant  in  the  porous  ceramic 
coating  resulted  in  high  surface  thermal  absorptivity  as  the 
infiltrant  decomposed  and  charred  on  heating.  Analyses  show 
that  the  absorptivity  must  be  approximately  0.2  in  order  for 
the  thin  ceramic  coating  to  provide  the  necessary  thermal  pro¬ 
tection  for  the  required  nuclear  criteria  environment.  Emis- 
sivity  values  of  0.2  have  been  obtained  in  previous  program 
bests  of  aluminum  oxide  and  zirconia  ceramic  coatings.  Measure¬ 
ments  of  the  electrical  properties  (dielectric  constant  and 
loss  tangent)  showed  the  ceramics  to  have  adequate  properties 
for  possible  radar  applications,  e.g.,  loss  tangent  less  than 
0.03.  A  major  concern  for  the  ceramics  is  a  potential  thermal 
stress  failure  problem.  However,  analyses  suggest  the  possible 
solution  of  pre-heating  during  processing  to  minimize  tensile 
stresses  during  exposure  to  the  high  heat  flux  environment 
from  nuclear  bursts. 

The  silicone  ablator  shows  much  promise  for  protecting 
radar  components.  This  material  forms  a  low  carbon  content 
char  which  is  reasonably  strong.  The  excellent  insulative  prop¬ 
erties  and  high  temperature  char  characteristics  result  in  good 
thermal  resistance.  Measurements  of  electrical  properties  showed 
low  loss  tangents  and  dielectric  constants  for  the  virgin  mater¬ 
ials.  Although  the  char  properties  are  not  as  good  as  the  vir¬ 
gin  properties ,  analyses  show  that  the  material  could  be  used 
as  a  coating  in  integral  phased  array  radar  designs. 


The  following  sections  present  the  results  of  test  and 
analysis  work  performed  during  this  study.  Section  2.0  des¬ 
cribes  the  phased  array  radar  hardening  requirements  including 
a  definition  of  the  nuclear  environment,  contractor  phased  array 
radar  design  studies  and  material  hardness  requirements.  Sec¬ 
tion  3.0  describes  the  candidate  materials  selected  for  eval¬ 
uation  in  this  study.  In  Section  4.0,  the  thermal  analysis 
and  test  results  are  presented.  The  electrical  properties 
test  approach  and  results  are  presented  in  Section  5.0.  Mech¬ 
anical  properties  tests  and  analyses  are  described  in  Section 
6.0  and  the  conclusions  and  recommendations  resulting  from 
this  study  are  presented  in  Section  7.0.  Detailed  test  and 
analysis  results  are  presented  in  Appendices  A-G, 
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2 . 0  HARDENING  REQUIREMENTS 

2 . 1  Nuclear  Environment 

The  hardened  radar  will  be  designed  to  withstand  the 

nuclear  weapon  environments  associated  with  three  successive 

nuclear  bursts  with  a  time  spacing  sufficient  for  radar  cool 

down  between  bursts.  The  free  field  peak  design  overpressure 

hypothesized  is  150  psi.  The  nuclear  environments  are  based 

on  IDEA  code  calculations  documented  in  References  1  and  2. 

Figures  2.1  through  2.5  present  the  overpressure,  dynamic 

pressure,  local  air  temperature,  particle  velocity  and  density 

for  each  burst.  Note  that  the  nuclear  fireball  sweeps  over 

the  radar  at  about  0.7  seconds  after  burst.  The  thermal 

radiation  flux  is  given  in  Figure  2.6  and  the  fluence  is  given 

in  Figure  2.7.  The  shock  arrival  time  is  0.48  seconds  after 

2 

burst,  and  approximately  250  cal/cm  (^10  percent  of  total 

fluence)  preceeds  the  shock.  The  thermal  flux  at  shock 

2 

arrival  is  about  1500  cal/cm  -sec.  Figure  2.8  presents  the 
fireball  equivalent  blackbody  temperature  useful  for  spectral 
analysis . 

2.2  Phased  Array  Radar  Design  Requirements 

The  phased  array  radar  must  survive  the  nuclear  ther¬ 
mal  and  blast  environment  and  remain  operational.  The  system 
must  still  be  capable  to  electromagnetically  function  in  an 
undegraded  manner  after  exposure  to  the  three  consecutive  nu¬ 
clear  events.  The  damage  after  a  single  nuclear  burst  must 
not  detract  from  the  capability  of  the  radar  to  function  or 
survive  a  second  or  third  burst. 
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Figure  2.4.  Air  Particle  Velocity  History 
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Figure  2.5.  Air  Density  History 
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Figure  2.8.  Effective  Blackbody  Temperature 
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The  thermal  radiation  pulse  from  the  threat  nuclear  burst 
heats  the  radar  emitting  surfaces  and  radar  structure.  Insu¬ 
lation  can  be  used  to  protect  the  radar  surfaces  from  thermal 
radiation  but  the  requirements  imposed  by  the  radiating  sur¬ 
face  are  different  from  those  imposed  by  the  support  structure 
The  insulation  selected  for  protecting  the  emitting  surface 
must  not  degrade  radar  performance  in  its  unheated  state  and 
its  char  and  products  of  decomposition  must  have  minimal 
effect  on  performance  after  being  heated  by  a  nuclear  burst. 
The  loss  tangent  and  dielectric  constants  of  the  insulator  and 
decomposition  products  must  be  selected  with  regard  to  their 
effect  on  performance.  Also,  the  thickness  allowable  is 
determined  by  the  electrical  properties  and  radar  interactions 
Thus,  many  well  proven  ablators  developed  for  missiles  and 
reentry  vehicles  are  not  suitable  for  use  in  protecting  rad¬ 
iating  surfaces  because  their  pre-exposure  or  post  exposure 
electrical  properties  are  incompatible  with  radar  performance. 
Free  carbon  in  the  ablator  or  char  is  detrimental  because  it 
is  a  conductor.  The  ablators  that  can  be  used  to  protect  the 
non-radiating  structural  surfaces  of  the  radar  have  no  such 
restriction.  They  can  be  of  any  suitable  thickness  and  their 
electrical  properties  are  unimportant.  They  must  simply 
survive  and  remain  weather  resistant. 

The  insulators  or  ablators  chosen  to  protect  the  radar 
radiating  surfaces  and  structural  surfaces  must  be  able  to 
withstand  the  environment  (rain,  snow  and  ice)  without  deter¬ 
iorating  and  they  must  be  capable  of  withstanding  the  airblast 
pressures  which  will  tend  to  deform  them.  They  must  be  cap¬ 
able  of  transmitting  the  pressure  to  a  load  carrying  member. 
Furthermore,  the  ablator  should  not  absorb  water  to  any  depth 
because  the  steam  formed  by  the  heat  from  a  nuclear  burst  may 
destroy  the  ablator  or  if  it  collects  in  the  bond  area,  it  may 
blow  the  ablator  from  the  protected  structure.  Furthermore 
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trapped  or  absorbed  water  will  freeze  and  tend  to  breakup 
the  insulator  when  ice  is  formed. 

Thermally,  the  ablator  must  be  adequate  to  prevent 
wraping  of  the  structure  due  to  differential  thermal  expansion. 

Of  course,  it  should  prevent  melting,  cracking,  or  electrical 
property  and  mechanical  property  changes  in  the  protected 
structures.  For  example,  although  most  steels  could  withstand 
heating  to  400°F  for  a  short  time,  the  thermal  stresses 
introduced  by  suddenly  heating  the  surface  to  400°F  even  for 
a  short  time  could  cause  wraping  or  failure  of  a  structure. 

Thus,  the  establishment  of  thermal  performance  require¬ 
ments  and  electrical  performance  requirements  for  materials 
to  protect  radars  from  nuclear  thermal  radiation  and  blast  envi¬ 
ronments  is  not  a  simple  task.  The  requirements  can  only  be 
determined  after  establishing  selected  materials  performance 
levels  and  evaluating  the  impact  of  these  materials  on  the  radar 
performance.  This  analysis  was  performed  during  this  program  and 
is  reported  in  Section  5  of  this  report. 
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3.0 


CANDIDATE  HARDENING  MATERIALS 


3 . 1  Ceramic  Coating 

During  previous  studies* for  BMDATC,  SAI  has  investi¬ 
gated  the  use  of  thin  ceramic  coatings  for  thermal  protection 
under  severe  environments.  The  mission  for  the  ceramics  was 
to  protect  a  metal  missile  structure  from  atmospheric  heating 
in  the  lower  atmosphere  while  flying  through  the  air  at  vel¬ 
ocities  as  high  as  16,000  fps.  Heating  rates  experienced  by 

2 

the  ceramic  coating  would  be  as  high  as  3500  BTU/Sec.Ft  .  (949 
2 

cal/cm  sec.)  Numerous  tests  were  conducted  in  the  SAI  arc jet 

heater  facility  in  which  thin  layers  of  various  ceramics 

attached  to  thin  metal  substructures  were  heated  at  rates  of 
2  2 

2000  BTU/Sec.Ft  (542  cal/cm  Sec) .  The  specimens  were  repre¬ 
sentative  of  missile  structures  in  which  the  substructure  was 
a  MAR-M2 0 0  superalloy  or  TZM  molybdenum  alloy.  Plasma  spray 
and  flame  spray  processes  were  used  to  spray  molten  ceramic 
droplets  onto  the  surface  of  the  base  metal  for  thermal  protec 
tion.  The  thickness  of  ceramic  varied  between  0.005  to  .090 
inches  and  metal  thickness  was  .090  inches  for  the  MAR-M200 
and  .125  inches  for  the  TZM.  Ceramics  evaluated  in  the 
thermal  test  program  included: 


Chromium  cxide 
Aluminum  oxide 
Zirconium  oxide 

Zirconium  oxide  plus  yttrium  oxide 
Zirconium  oxide  plus  magnesium  oxide 
Aluminum  oxide  plus  calcium  oxide 
Aluminum  oxide  plus  magnesium  oxide 

The  thermal  tests  conducted  in  the  SAI  arejet 
successfully  demonstrated  the  thermal  performance  predicted 
for  the  ceramic  coatings.  Some  ceramics  were  infiltrated 
with  epoxy  which  acted  as  an  ablator  and  reduced  the  back 
surface  temperature  of  the  metal  in  a  very  effective 


*  See  References  3  and  4 
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manner.  This  is  illustrated  in  Figure  3.1  which  shows  that 
the  back  surface  temperature  of  the  ZrC>2  -  Y2C>3  ceramic  pro¬ 
tected  specimen  is  1500°F  lower  at  the  end  of  test  (1.8  sec- 

2 

onds  at  2000  BTU/sec.  Ft  )  when  it  is  infiltrated  with  epoxy 
than  when  it  is  uninfiltrated.  The  figure  also  shows  that  at 
0.8  seconds  the  ceramic  protected  specimens  are  heated  to 
approximately  300°F  whereas  a  bare  unprotected  specimen  is  at 
2000°F . 

An  optical  pyrometer  was  used  to  obtain  front  surface 
temperature  of  the  specimens  during  heating.  The  alumina  and 
zirconia  specimens  all  exhibited  an  emissivity  of  approximately 
0.2.  Thus,  when  the  ceramic  coatings  were  postulated  for  use 
as  radar  protection,  the  heating  results  were  compared  to 
radar  environments. 

2  2 

Because  the  2000  BTU/Sec. Ft  (542  cal/cm  Sec.)  is  a 
measure  of  actual  input  flux  it  must  be  compared  to  the  actual 
input  from  thermal  radiation  with  emissivity  of  0.2.  The  test 
conditions  corresponded  to  nuclear  thermal  equivalent  radia¬ 
tions  as  shown  below: 

2 

Absorbed  Flux  »  2700  Cal/cm  Sec. 

2 

Absorbed  Fluence  >  300  Cal/cm 

Actual  nuclear  radiation  threat 

2  2 

Absorbed  Flux  =  1650  cal/cm  Sec  x  0.2  =  3330  cal/cm  Sec. 

Because  the  flame  and  plasma  sprayed  ceramic  coatings  had  suc¬ 
cessfully  withstood  heat  fluxes  greater  than  those  expected 
for  the  radar  environment  and  because  they  were  a  dielectric, 
they  were  proposed  as  candidate  radar  hardening  materials. 

Also,  because  infiltrants  were  so  effective  in  protecting  the 
structure  in  the  arc  jet  tests,  it  was  proposed  to  use  them 
for  radar  protection  with  the  additional  constraint  of  using 
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1/cm  -sec 


an  infiltrant  that  did  not  form  large  amounts  of  free  carbon 
during  decomposition. 

The  ceramic  coating  materials  are  applied  to  a  metal 
surface  by  two  processes.  In  either  process  the  metal  surface 
must  be  prepared  by  abrasive  grit  blasting  to  roughen  the 
surface.  In  the  flame  spray  process  called  the  Norton  Rokide 
process,  solid  bars  of  ceramic  (similar  to  welding  rods)  are 
fed  into  an  oxyacetylene  torch  which  melts  the  rod.  A  com¬ 
pressed  air  nozzle  sprays  the  molten  ceramic  droplets  onto  the 
roughened  metal  surface  where  they  solidify  and  adhere  to  the 
surface  by  mechanical  action.  There  is  no  infusion  of  base 
material  or  ceramic. 

In  the  second  process,  finely  ground  ceramic  powder  is 
fed  into  the  center  chamber  of  an  arc  jet  where  it  melts  and 
flows  along  with  the  superheated  plasma  out  of  the  nozzle  and 
it  strikes  the  workpiece  where  it  solidifies  and  adheres  to  the 
surface. 

The  quality  and  type  of  deposit  is  dependent  on  factors 
such  as  plasma  gas  flow  rate,  electrical  power,  powder  flow 
rate,  distance  from  plasma  gun  to  workpiece,  traverse  speed, 
powder  size,  etc.  Variations  in  the  above  parameters  allow  the 
applied  ceramic  to  become  porous  so  that  an  ablator  can  be 
infiltrated  into  the  ceramic  matrix  after  deposition. 

Another  variable  in  the  process  which  was  used  in  the 
previous  SAI/BMDATC  program  involved  the  deposition  of  a  thin 
layer  of  nichrome  or  molybdenum  as  a  buffer  between  the  plasma 
sprayed  ceramic  and  the  base  metal.  For  some  applications, 
extremely  strong  joints  were  made  between  the  ceramic  and  base 
metal  by  depositing  the  pure  thin  nichrome  buffer  and  then  by 
progressively  changing  the  mix  to  pure  ceramic  powder.  This 
mixture  of  ceramic  and  nichrome  was  then  heated  in  a  hydrogen 
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furnace  to  diffuse  the  nichrome  into  the  base  metal  and  ceramic 
to  form  a  cermet  layer  approximately  .008  inches  thick.  For 
the  radar  experiment  described  in  this  report,  the  use  of  the 
cermet  process  was  not  investigated  because  preliminary  inves¬ 
tigated  showed  it  to  be  bad  from  a  radar  transmission  stand¬ 
point.  Whether  this  initial  assumption  is  still  valid  can  not 
be  answered  at  this  time.  No  thermal  bond  failures  were  ob¬ 
served  with  this  type  of  process.  The  ceramic  coatings  inves¬ 
tigated  are  presented  in  Table  3,1, 

3 . 2  Bulk  Window  Materials 

As  a  part  of  this  program,  SAI  investigated  materials 
that  could  be  used  for  windows  or  plugs  for  several  BMDATC 
radar  concepts.  The  material  had  to  be  strong  enough  to  take 
the  overpressure  and  had  to  be  able  to  resist  the  thermal 
exposure  without  significant  change  in  thickness  or  change  in 
radar  signal  propagation  properties.  Investigation  of  the 
available  window  materials  and  tests  previously  conducted  in 
the  French  Solar  Furnace  showed  that  most  of  the  available 
candidates  failed  by  thermal  shock.  As  a  result,  the  litera¬ 
ture  and  manufacturers  were  consulted  to  see  if  there  were  new 
developments  which  could  be  useful  for  BMDATC  radars.  A 
number  of  candidate  materials  were  obtained.  These  are  listed 
in  Table  3.2. 


3.3 


Bulk  Ablative  Materials 


Unlike  the  bulk  window  materials,  the  bulk  ablative 
materials  need  have  no  special  electrical  properties.  They  are 
to  be  used  to  protect  structural  portions  of  the  radar  from 
thermal  radiation.  They  should  be  effective  ablators  and 
should  be  relatively  inexpensive.  Some  may  be  very  poor  abla¬ 
tors  but  may  have  excellent  insulation  properties  and  resis¬ 
tance  to  thermal  shock.  They  must  be  weather  resistant  and 
moisture  proof.  Structurally,  the  bulk  ablators  must  be  cap¬ 
able  of  transmitting,  by  compression  and  shear,  the  blast  load 
imparted  by  the  shock  wave.  This  will  be  less  than  lOOOpsi  and 
most  materials  are  capable  of  doing  this.  Some  of  the  candidate 
bulk  ablators  were  rubbers,  silicone  rubbers,  firebrick  and 
commercial  insulators.  Table  3.3  lists  the  bulk  ablators. 
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Table  3.1.  Ceramic  Coatings 


Coating  Materials 
A1203 

AI2O3  *  Ti02 
Zr02  *  CaO 
Zr02  *  MgO 
Zr02  *  Y203 
Spinel 

Infiltrant 

Sylgard  184  (Silicone  Resin) 


Table  3.2.  Bulk  Window  Materials 

Boron  Nitride  (NBR) 

Boron  Nitride  (HP) 

Boron  Nitride  (M) 

Alumina  (AD-995) 

Alumina  (AD-995)  with  porous  alumina  coating 
Alumina  (Durafrax) 

Alumina  (M) 

Alumina  (A) 

Beryllia 

Fused  Silica  7941 


Table  3.3  Ablator/ Insulator  Materials 


Silicone  Rubber  (Down  Corning  Silastic  E) 

Alumina  Firebrick 

Zirconia  Firebrick 

Silfrax  Firebrick 

Teflon 
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4.0 


THERMAL  PROPERTIES  ANALYSIS  AND  TEST 


Thermal  models  were  developed  for  the  ceramic  coatings , 
with  and  without  infiltrant,  and  for  the  bulk  "window"  materials 
without  ablation.  The  bulk  silicone  rubber  model  has  been 
formulated;  however,  it  must  be  updated  with  experimental  data 
before  it  will  give  satisfactory  results.  The  models  were  used 
to  predict  the  performance  of  the  ceramic  coatings  in  the  nuclear 
environment,  to  predict  the  material  sample's  response  in  the 
thermal  tests  and  hence  to  size  the  samples  and  define  instru¬ 
mentation  requirements,  and  to  parametrically  investigate  the  in¬ 
fluence  of  various  thermophysical  properties. 

All  of  the  materials  investigated  were  tested  in  the 
Thermal  Radiation  Simulator  (TRS)  and  a  limited  number  of  materials 
were  tested  in  the  French  Solar  Furnace  (FSF) .  The  following 
paragraphs  discuss  the  thermal  modeling,  the  thermal  tests  and 
the  results  obtained. 

4. 1  Thermal  Analysis 

An  analytical  model  was  used  to  predict  the  transient 
thermal  response  of  ablating  ceramic  coatings  and  silicone 
rubber.  The  method  used  and  the  models  developed  are  as  follows: 

4.1.1  Ceramic  Coating  Thermal  Analysis 

The  ceramic  coating  being  modeled  consists  of  a  porous 
ceramic  bonded  to  a  steel  substrate.  The  porous  ceramic  can 
either  contain  or  not  contain  a  silicone  rubber  infiltrant.  The 
basic  model  to  be  described  contains  the  infiltrant.  If  the 
infiltrant  is  not  included,  the  model  is  simplified  accordingly 
by  removing  the  infiltrant' s  contribution  to  the  thermal  response. 

A  physical  description  of  the  material  and  the  thermal 
model  is  shown  in  Figure  4.1.  The  primary  heat  input  in  the 
nuclear  environment  is  time  dependant  radiation  and  it  is  the 
only  heat  input  considered.  The  infiltrated  ceramic  coating 
acts  as  a  thermally  protective  material  through  the  following 
heat  transfer  mechanisms: 


1.  The  majority  of  the  incoming  heat  is  reflected  due 
to  a  low  surface  absorptivity. 

2.  The  surface  reradiates  at  a  high  temperature. 

3.  Heat  is  conducted  into  the  coating  and  absorbed. 

4.  The  infiltrant  decomposes  to  its  product  gases, 
absorbing  its  heat  of  decomposition  ar d  also  absorb¬ 
ing  heat  from  the  ceramic  as  it  escapes  to  the 
surface. 

5.  The  ceramic  melts,  absorbing  its  change  of  phase 
heat. 

In  the  current  model,  the  melt  layer  of  ceramic  is  assumed  to  be 
removed  from  the  surface  immediately  due  to  the  nuclear  environment 
air  flow.  This  is  a  conservative  model,  because  if  the  melt 
layer  did  remain,  it  would  offer  additional  thermal  protection 
due  to  its  thermal  resistance,  thermal  mass,  and  possible  vapor¬ 
ization. 

Another  parameter  not  included  is  possible  radiation 
blockage  due  to  a  vapor  cloud  formed  by  the  reactant  gases.  If 
this  cloud  were  able  to  remain  in  front  of  the  venting  surface 
during  a  portion  of  the  heating,  it  could  absorb  some  of  the 
incident  radiation. 

The  analytical  model  was  developed  using  the  MRCAPER 
computer  code,  of  which  the  basic  theory  is  presented  in  Ref¬ 
erence  5.  A  schematic  of  the  model  is  shown  in  Figure  4.1. 

As  the  outer  surface  of  the  ceramic  is  heated,  heat  is  conducted 
inward  and  the  infiltrant  decomposes  as  its  temperature  of  de¬ 
composition  is  reached.  The  decomposing  infiltrant  evolves 
gases  which  perculate  to  the  ceramic  surface  through  the  pores. 
Thus,  the  infiltrant  is  in  a  steady  state  condition  of  decom¬ 
position  below  the  surface  of  the  ceramic.  As  the  ceramic  sur¬ 
face  itself  reaches  its  melt  temperature,  it  is  removed,  thus 
causing  recession  of  the  surface. 

The  governing  equations  for  energy  conservation,  infiltrant 
decomposition  rate,  and  pyrolysis  gas  mass  flux  are 


/ 
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IE  =  -  Z  (p~pc)ni  A.e“Ei/RT 


(4.2) 


r°  3P 

mg  =  -  4  9t  d* 


(4.3) 


where  A.  =  frequency  factor  for  Pyrolyzing  reaction 


specific  heat 

activation  energy  for  the  i—  pyrolyzing  state 
enthalpy 

thermal  conductivity 
surface  recession 
mass  flux 

reaction  order  for  i—  pyroloyzing  reaction  state 

universal  gas  constant 

time 

coordinate  normal  to  surface 
density 


Subscripts 


c  =  of  ceramic  without  infiltrant  or  of  char 

g  =  of  gas 

s  =  of  solid 

v  =  virgin  material 

The  heat  removed  by  the  melting  of  the  ceramic  is  handled  as  a 
boundary  condition  and  removed  from  the  incoming  heat. 

These  equations  then  form  the  basis  of  the  theoretical 
model  that  is  solved  by  the  MRCAPER  computer  code.  The  model 
has  been  developed  for  porous  alumina  and  zirconia  coatings  with 
and  without  the  silicon  resin  infiltrant,  Dow  Corning  Sylgard 
184.  The  material  thermal  properties  are  given  in  appendix  A. 

The  model  has  been  used  to  predict  the  performance  of  alumina 
with  infiltrant  in  the  nuclear  environment,  Section  4. 1.1. 2,  to 
predict  the  performance  and  size  the  samples  and  instrumentation 
in  the  thermal  tests.  Section  4.3.2  and  4.4.3,  and  used  in  a 
parametric  investigation  of  the  effects  of  infiltrant  and  various 
th'jrma  1/physical  parameters,  Section  4. 1.1.1  and  Appendix  F. 


-26- 


4. 1.1.1  Ceramic  Coating  Thermal  Analysis 

A  parametric  analysis  was  made  of  porous  alumina  and 
zirconia  coatings,  with  and  without  silicone  resin  infiltrant 
bonded  on  a  stainless  steel  substrate.  The  purpose  of  the 
analysis  was  to  size  the  material  specimen  for  the  thermal 
tests,  to  provide  temperature  profiles  for  the  structural 
analysis,  and  to  investigate  the  effect  of  various  material 
property  and  design  variables.  The  heat  flux  utilized  was 
the  TRS  expected  maximum  heating  environment.  This  environ¬ 
ment  is  similar  to  the  nuclear  environment  for  a  material  with 
surface  absorptivity  of  0.2. 

The  thermal  environment  description  and  the  analysis 
results  are  presented  in  Appendix  P.  The  results  of  the  parametric 
investigation  of  variables  are  summarized  as  follows: 


1.  The  ceramic  acts  as  an  insulator  and  its  surface 
heats  to  its  melting  temperature  almost  immediately. 

2.  The  alumina  and  zirconia  perform  simlarly  thermally 
except  that  the  zirconia  has  a  higher  melt  temperature. 

3.  A  thicker  coating  exhibits  more  ablation  than  a  thin 
one;  however,  the  thicker  coating  also  has  more 
material  remaining  after  the  exposure. 

4.  Decreasing  the  porosity  of  a  non- inf iltrated  ceramic 
increases  its  conductivity  and  thermal  mass,  resulting 
in  a  lower  ceramic  temperature. 

5.  A  silicone  resin  infiltrant  in  the  porous  ceramic 
reduces  the  amount  of  ablation  of  the  ceramic,  as 
long  as  it  doesn't  increase  the  ceramic  absorptivity. 


4. 1.1. 2  Ceramic  Coating  Predicted  Performance  In  The  Nuclear 

Environment 


The  porous  ceramic  coating  model  was  used  with  the  design 
nuclear  radiant  flux  profile  to  predict  the  performance  of  15 
percent  porous  alumina  coating  on  1.0  inch  steel.  An  absorptivity 
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of  0.2  was  used.  The  resulting  surface  recession,  surface 
temperatures,  and  indepth  temperature  profiles  are  shown  for  a 
single  burst  in  Figure  4.2.  The  resulting  surface  recession  is 
approximately  0.02  inches  for  the  0.125  inch  sample.  Figure  4.3 
shows  the  surface  temperature  profile  and  cumulative  surface 
recession  for  an  exposure  to  3  bursts,  with  sufficient  time 
between  bursts  to  return  to  ambient  temperature.  This  coating 
would  be  able  to  withstand  3  nuclear  bursts. 

A  low  surface  absorptivity  is  critical.  Figure  4.4 
compares  the  results  of  the  0.2  absorptivity  ceramic  with  a  0.5 
absorptivity  ceramic.  The  ceramic  with  an  absorptivity  of  0.5 
has  a  single  burst  recession  of  0.09  inches.  Thus,  it  could  not 
withstand  exposure  to  multiple  bursts. 

4.1.2  Silicone  Rubber  Ablation  Model 

The  physical  description  of  the  silicone  rubber  coating 
in  the  nuclear  environment  is  shown  in  Figure  4.5.  The  incoming 
heat  decomposes  the  silicone  rubber  and  forms  a  char  layer.  If 
this  char  layer  stays  intact,  it  forms  a  low  conductivity  layer 
that  protects  the  remaining  virgin  silicone  rubber,  which  is 
still  decomposing,  but  at  a  much  slower  rate.  The  char  layer 
has  a  high  absorptivity  so  that  it  absorbs  most  of  the  incoming 
energy;  but  likewise,  it  has  a  high  emissivity  and  high  surface 
temperature  such  that  a  large  amount  of  the  absorbed  energy  is 
reradiated  away. 

The  thermal  model  again  uses  the  MRCAPER  code  and  the 
theory  of  Section  4.1.1.  The  portion  of  the  material  at  the 
lower  density,  P  ,  is  now  the  char  layer  (Figure  4.1).  The 
necessary  material  thermal  and  kinetics  data  were  evaluated  and 
are  presented  in  Appendix  A.  The  model  has  not  yet  been  checked 
out  and  made  operational,  although  preliminary  investigations 
indicated  that  the  char  layer  which  would  develop  in  the  thermal 
tests  would  remain  relatively  thin  with  the  remainder  of  the 
silicone  rubber  remaining  in  its  virgin  state. 
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4.1.3  Bulk  "Window"  Material  Thermal  Modeling 


The  thermal  response  of  the  bulk  materials  in  the  TRS 
tests  was  predicted  to  (1)  indicate  the  backface  temperatures 
expected  in  order  to  specify  the  temperature  indicating  paint 
ranges  required,  and  (2)  indicate  the  type  of  thermal  response 
expected,  i.e.,  surface  melting  or  interior  absorption  of  all  of 
the  heat.  The  response  was  calculated  using  a  transient  conduc¬ 
tion  code  which  does  not  take  surface  decomposition  into  account. 
This  code  is  described  in  Reference  6. 

4.2  Comparison  Of  Thermal  Simulators  Facilities  Capabilities 

In  materials  research,  there  are  a  number  of  different 
types  of  heating  schemes  in  use.  These  types  are  briefly  discussed 
to  show  their  characteristics. 

4.2.1  Solar  Furnace 

The  solar  furnace  uses  radiant  energy  from  the  sun.  It 
collects  the  radiation  in  the  form  of  visible  light,  infra  red 
and  ultra  violet  on  a  large  area  covered  by  mirrors.  The 
mirrors  are  then  focussed  onto  the  target  by  curved  mirrors 
acting  as  a  lens.  Some  of  these  systems  use  first  surface 
mirrors  which  reflect  all  of  the  energy  to  the  target  but  most 
of  them  use  second  surface  mirrors  which  internally  absorb  ultra 
violet  and  other  radiation  wavelengths  depending  on  the  glass 
characteristics.  These  simulators  are  useful  but  their  perfor¬ 
mance  and  characteristics  depend  on  time  and  aging. 

The  French  Solar  Furnace  facility  at  Odeillo,  France  is 
the  most  powerful  thermal  radiation  source  available.  It  has 
a  second  surface  mirror  system  which  means  that  the  aluminized 
reflecting  surface  is  behind  the  glass  as  in  a  common  household 
mirror.  Thus  the  rays  of  the  sun  must  pass  through  the  glass  to 
strike  the  rear  reflecting  surface  and  then  pass  out  of  the 
glass  to  exit.  This  causes  ultra  violet  to  be  absorbed  in  the 
glass  and  also  refracts  the  light  rays  to  some  degree.  In 
addition,  as  the  mirror  system  ages,  its  glass  cover  becomes 


less  transparent  due  to  sandblasting  by  wind  driven  particles 
and  deterioration.  Therefore,  although  the  French  solar  furnace 
had  a  maximum  flux  capability  of  400  cal/cm  sec  on  a  bright 
sunny  day,  it  has  deteriorated  to  approximately  50%  of  that 
value;  however,  it  is  being  refurbished.  Fluence  available 
depends  on  the  quality  and  time  exposure  of  sunlight. 


Testing  is  accomplished  by  focussing  the  hot  spot  to  a 
diameter  in  the  order  of  three  inches.  The  radiation  is  blocked 
from  the  specimens  by  shutters  that  can  open  or  close  upon 
demand.  For  precise  timing,  the  specimen  (generally  2  inches  in 
diameter)  is  exposed  to  the  hot  spot  by  rotating  a  disk  holding 
the  specimen  past  a  water  cooled  iris  by  means  of  a  geneva 
mechanism.  It  is  possible  to  instrument  the  specimens  with 
thermocouples  and  to  observe  the  specimens  with  cameras  and 
optical  pyrometers. 


A  second  solar  furnace  designed  for  materials  research  is 

located  at  the  White  Sands  Missile  Range.  It  has  a  flux  capability 

2 

on  the  order  of  80  cal/cm  sec.  It  is  similar  to  the  French 

facility  but  is  smaller  and  in  a  poor  state  of  repair.  A  visit 

was  made  to  WSMR  to  inspect  the  facility  and  its  instrumentation. 

2 

A  400  cal/cm  sec  facility  is  in  the  design  stage  at  WSMR  and  is 
scheduled  for  construction  in  1979.  It  consists  of  a  large 
parabolic  radar  dish  lined  with  mirrors  which  is  capable  of 
tracking  the  sun  directly. 

4.2.2  Plasma  Jet 

Much  materials  heating  research  is  conducted  in  a  plasma 
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jet  or  arc  jet  heating  facility.  Heat  is  generated  by  striking 
an  arc  between  an  internal  and  external  electrode  and  passing  a 
gas  between  the  electrodes.  The  gas  is  heated  by  the  electrical 
discharge  and  becomes  a  plasma  which  exits  the  nozzle  and  impinges 
on  the  test  specimen.  Very  high  heating  rates  can  be  generated 
in  these  arc  jet  facilities, but  the  heating  action  does  not 
simulate  radiation  heating  from  a  nuclear  burst. 

The  plasma  jet  has  sufficient  kinetic  energy  to  sweep 
away  melt  as  soon  as  it  is  formed.  Thus,  the  heat  absorbing 
capacity  of  the  melt  on  the  surface  is  lost  in  that  no  account 
is  taken  of  the  heat  required  to  raise  the  melt  to  boiling 
temperature  and  the  heat  required  to  actually  boil  the  melt. 

Also,  the  effect  of  effluent  gases  and  vaporized  material  in 
screening  radiation  is  not  measurable  because  the  gases  are 
driven  away  by  the  plasma  and  there  is  no  radiation  heating. 
Furthermore,  because  many  materials  of  interest  may  be  or  are 
transparent  to  radiation  to  a  degree,  heating  in  depth  cannot  be 
simulated.  Therefore,  the  plasma  jet  facilities  were  not  recom¬ 
mended  for  the  tests  in  this  study.  However,  the  plasma  jet  may 
be  a  useful  way  to  screen  materials  for  susceptability  or  resis¬ 
tance  to  thermal  shock. 

4.2.3  Glow-Bar  Heaters 

Considerable  radiation  heating  can  be  obtained  by  electrically 
heating  bars  or  plates  of  tungsten  or  graphite  near  a  specimen. 

The  heating  must  be  in  a  vacuum  to  prevent  oxidation  of  the 
radiator.  This  type  of  facility  was  not  considered  because  the 
radiating  surfaces  evaporate  large  quantities  of  graphite  or  tung¬ 
sten  which  will  contaminate  the  specimen  being  irradiated  and 
change  its  absorbtivity  and  electrical  properties.  The  radiation 
heaters  generally  work  in  a  vacuum  to  reduce  oxidation  of  the 
radiating  element.  Maximum  temperatures  are  limited  to  values 
well  below  the  melting  point  of  the  radiator  because  evaporation 
or  sublimation  near  the  melting  point  becomes  very  severe. 

Generally,  tungsten  radiators  may  operate  below  5000°F  and 
graphite  radiators  below  6000°F.  This  makes  the  sources  marginal 


(even  if  deposits  of  radiator  vapor  are  no  problem)  because  the 
ceramics  melt  and  boil  at  high  temperatures.  Alumina,  A120-j, 
melts  at  3700°F  and  boils  at  5396°F  while  zirconia,  Zr02  melts 
at  4800°F  and  boils  at  5890°F.  Unless  emissivity  and  absorbtivity 
approach  one,  these  temperatures  will  never  be  reached  in  a  glow 
heater  type  of  facility.  Other  variations  of  the  glow  heater 
facilities  are  banks  of  quartz  lamps  with  radiating  filaments 
all  focussed  on  the  test  specimen.  Their  maximum  temperature  is 
approximately  3000°F  which  is  too  low  to  melt  and  ablate  many  of 
the  materials  of  interest  in  this  program,  hence,  this  limits  the 
usefulness  of  such  a  facility. 

4.2.4  SAI  Thermal  Radiation  Simulator 

Under  contract  to  DNA,  SAI  has  developed  and  tested  a 
thermal  radiation  simulator  based  on  the  flashbulb  principal. 

Thin  plastic  bags  are  filled  with  gaseous  oxygen.  Fine  aluminum 
powder  is  dispersed  throughout  the  interior  of  the  oxygen  filled 
bag  and  is  ignited  by  an  electrically  initiated  pyrotechnic.  The 
combustion  of  aluminum  in  oxygen  generates  a  high  temperature 
gas  which  glows  and  radiates.  A  large  fireball  is  formed  which 
expands  and  cools  down  as  time  progresses.  The  fireball  rises 
and  leaves  the  ground  within  a  few  seconds. 

Although  the  aluminum  powder  combustion  type  of  thermal 
simulator  was  fully  developed  for  DNA  in  small  sizes,  a  large 
24  bag  array  was  designed  and  fabricated  for  this  project. 

Figure  4.6  schematically  illustrates  the  bag  array.  The  framework 
was  fabricated  from  heavy  steel  structural  members.  The  array 
consisted  of  three  rows  of  8  bags  each  as  shown  in  the  figure. 

Oxygen  was  supplied  in  liquid  form  and  was  stored  in  a 
reservoir  as  illustrated  in  figure  4.7.  An  air- liquid  heat 
exchanger  was  used  to  convert  the  liquid  oxygen  to  gaseous 
oxygen.  It  was  piped  from  the  storage  and  converter  area  to  the 
test  area.  Instrumentation  consisting  of  calorimeters  and 
thermocouples  were  located  on  the  specimen  mount  and  in  several 
points  outside  of  the  fireball. 


SAI  Thecal  Radiation  Simulator 

—  - - 


Figure  4,6 


Based  on  previous  tests  with  smaller  arrays  (fewer  bags) , 
the  predicted  total  maximum  flux  was  400-1000  cal/cm  -sec  within 
the  fireball,  of  which  350  cal/cm2-sec  is  due  to  radiation.  Flu- 
ence  was  predicted  to  be  403  cal/cm2  outside  the  fireball.  Act¬ 
ual  performance  was  lower  than  expected  and  is  described  in  Sec¬ 
tion  4.3. 

The  primary  advantages  attributed  to  the  flashbulb  type 
thermal  radiation  simulator  are  its  very  high  heat  flux,  its 
large  available  area  for  testing,  and  its  ability  to  radiate 
over  the  desired  frequency  band. 

4. 3  Thermal  Radiation  Simulator  Tests 

Two  thermal  tests  were  attempted  with  the  SAI  Thermal 
Radiation  Simulator  (TRS) .  The  first  test  was  considered  a  no¬ 
test  because  of  a  failure  in  the  simulator  ignition  circuit  which 
resulted  in  very  low  heating.  The  second  test  generated  a 
successful  fireball;  however,  the  heat  in  the  vicinity  of  the 
sampler  was  less  than  10  percent  of  what  was  expected,  and  the 
results  are  quite  limited.  This  section  describes  the  test 
facility,  the  test  procedures,  the  pre-test  analysis,  and  the 
results. 

4.3.1  Description  Of  The  TRS  Facility 

The  TRS  simulates  the  nuclear  event  thermal  radiation  by 
generating  a  large  fireball  through  combustion  of  aluminum 
powder  in  oxygen.  The  facility  used  consisted  of  a  set  of  24 
plastic  bags,  each  5  feet  in  diameter  by  20  feet  tall,  shown  in 
Figure  4.6.  The  bags  are  filled  with  oxygen  prior  to  the  test. 
Aluminum  is  then  injected  into  the  bags  by  nitrogen  gas  prior  to 
ignition.  At  the  time  required  for  the  aluminum  powder  to  fill 
the  bags,  the  mixture  is  ignited,  generating  a  large  fireball. 

A  TRS  calorimeter  printout  from  a  prior  test  is  shown  in 
Figure  4.8.  This  is  the  model  that  was  used  for  the  expected 

TRS  heat  output.  The  largest  prior  measured  maximum  heat  flux 

2  2 
was  225  cal/cm  sec  and  a  fluence  of  109  cal/cm  ,  measured  5 

feet  from  a  4  bag  array. 
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In  order  to  obtain  a  more  severe  heating  environment ,  a 
24  baq  array  was  selected.  The  mechanics  of  this  fireball 
are  unknown  and  the  facility  should  be  considered  very  exploratory. 
The  resulting  fireball  heat  output  is  a  function  of  the  fuel- 
oxidation  mixture,  the  array  layout,  the  calorimeter/sample 
location,  the  ignition 'timing,  and  other  factors  as  yet  unknown. 

For  the  24  bag  array,  the  TRS  facility  operator  predicted  a  flux 

2  2 
of  350  cal/  cm  sec  and  a  fluence  of  400  cal/cm  sec  due  to 

radiation  alone  (no  fireball  convective  heating) .  Within  the 

2 

fireball,  the  total  flux  was  predicted  to  be  400  to  1000  cal/cm 
sec  and  a  fluence  of  460  to  1150  cal/cm  .  This  very  large 
heating  within  the  fireball  was  assumed  to  partially  be  due  to 
condensation  of  aluminum  oxide  vapor.  This  expected  flux  was 
applied  to  the  heat  profile  of  Figure  4.8  with  the  resulting 
predicted  heat  profiles  shown  in  Figure  4.9. 

4.3.2  Pre-Test  Analysis 

Pre-test  analyses  were  run  on  the  materials  for  the 

purpose  of  sizing  the  samples,  specifying  the  instrumentation , 

and  to  predict  the  expected  material  thermal  response.  Radiation 

only  heating  utilized  the  lower  curve  of  Figure  4.9  giving  a 

2 

peak  flux  of  350  cal/cm  sec.  For  samples  within  the  fireball, 

2 

the  upper  curve  with  a  peak  flux  of  1000  cal/cm  sec  was  used. 

2  2  .  . 

Of  this  1000  cal/cm  sec,  350  cal/cm  sec  was  treated  as  radiation 

2  2 
and  650  cal/cm  sec  as  convection  or  condensation  (all  650  cal/cm 

sec  is  absorbed,  which  is  the  maximum  possible  heating) . 

The  radiation  only  heating,  with  the  sample  behind  a 
pyrex  glass  cover  was  used  so  that  the  effective  absorptivity  of 
the  samples  could  be  obtained.  The  full  fireball  heat  flux, 

1000  cal/cm  sec,  represents  the  maximum  heat  the  samples  could 
receive.  The  actual  heat  received  would  be  somewhat  less  depending 
on  the  sample  surface  temperature  and  that  part  of  the  flux 
due  to  convection  or  condensation.  However,  this  maximum  flux 
was  used  for  sizing  the  samples;  and,  it  also  more  nearly  repre¬ 
sents  the  nuclear  flux  for  the  parametric  study  in  Appendix  F. 
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Figure  4.9.  Predicted  24  Bag  TRS  Heat  Flux 
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The  results  of  the  pre-test  analysis  on  alumina  and 

zirconia  coatings  are  presented  in  Appendix  F.  Based  on  these 

results,  plus  a  minimum  thickness  requirement  for  the  plasma 

spraying  process,  a  0.125  inch  thick  stainless  steel  substrate 

was  selected  for  the  radiation  environment  only  samples ,  and 

0.25  inch  thick  stainless  steel  substrate  for  the  fireball 

exposure  samples.  At  maximum  heating,  the  radiation  only  samples 

would  attain  a  surface  temperature  of  about  1500°R  (below  melting) , 

and  a  backface  temperature  rise  of  about  300°F.  Those  exposed 

2 

to  the  full  fireball  flux  of  1000  cal/cm  sec  would  have  melting 
at  the  surface,  and  backface  temperature  rises  of  about  500  to 
1000°F . 

The  thermal  model  discussed  in  Section  4.1.2  is  not  yet 
capable  of  predicting  the  true  thermal  response  of  the  silicone 
rubber  samples.  However,  from  the  chemical  kinetics  and  heat  of 
decomposition  data  developed  and  presented  in  Appendix  A,  it  was 
estimated  that  0.25  inch  silicone  rubber  on  a  0.125  inch  thick 
stainless  steel  substrate  would  develop  a  thin  char  layer,  with 
most  of  the  silicone  rubber  remaining  in  its  virgin  state.  A 
very  small  temperature  rise  is  expected  for  the  steel  backface. 

4.3.3  Test  Procedures 

All  of  the  materials  discussed  in  Section  3.0  were  tested 
in  the  TRS  test.  The  specimens  consisted  of  six  basic  types: 

1.  3  inch  diameter  plasma  coated  ceramic  discs  with 
stainless  steel  backing  plates  and  thermocouples  on 
the  rear  metal  face.  (S  Band) 

2.  4  inch  square  plasma  coated  ceramic  plates  with 
steel  backing  plates  (X  band) 

3.  1  inch  diameter  plasma  coated  ceramic  discs  with 
steel  backing  plates  (X  band) 

4.  2  and  3  inch  diameter  and  2  inch  square  bulk  material 
samples 

5.  2  inch  square  tantalum,  columbium,  aluminum,  and 
stainless  steel  samples 
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6.  2  inch  diameter  teflon  cylinder  with  thermocouples 

imbedded  at  1/8,  1/4,  3/4,  and  1  inch  from  the 
heated  surface. 

The  samples  were  all  mounted  in  transite  holders.  These 
mountings  are  shown  in  Figures  4.10  and  4.11  for  the  3  inch 
diameter  and  4  inch  square  ceramic  coated  steel  substrate. 

Similar  holders  were  used  for  the  bulk  material  samples.  The 
holders  were  designed  to  restrict  heat  conduction  out  of  the 
samples  as  much  as  possible. 

The  test  plan  was  to  expose  4  samples  of  each  material 
during  the  first  test,  with  two  of  the  samples  mounted  behind 
pyrex  glass.  This  is  shown  in  Figure  4.12.  In  the  second  test, 
one  of  the  samples  behind  the  pyrex  glass  would  be  removed  as 
well  as  one  of  the  unprotected  samples.  Thus,  one  sample  would 
be  exposed  to  the  fireball  of  one  test  as  well  as  one  sample 
exposed  to  radiation  only,  and  one  sample  exposed  to  the  fireball 
of  2  tests,  and  similarly  for  one  sample  exposed  to  radiation 
only  in  2  tests.  Because  the  first  test  was  considered  a  no 
test,  only  the  second  test  setup  was  used,  and  one  sample  was 
exposed  to  the  fireball  of  one  effective  test,  and  one  sample 
was  exposed  to  the  radiation  of  this  test. 

The  samples,  in  their  transite  holders,  were  mounted  on  a 
steel  frame  5  feet  from  the  TRS  bags  at  a  height  of  10  feet. 

(See  Figure  4*6.)  Three  calorimeters  were  also  mounted  in  this 
frame.  Instrumentation  included  chrome 1- a lumel  thermocouples  on 
the  backface  of  all  the  3  inch  diameter  ceramic  coated  disks, 
and  temperature  indicating  paint  on  the  backface  of  all  samples. 
Standard  and  high  speed  cameras  were  used  to  take  motion  pictures 
of  the  samples  and  fireball. 

4.3.4  TRS  Test  Results 

As  discussed,  the  first  test  was  considered  a  no-test  and 
only  the  second  test  is  considered  in  these  results.  Although 
better  than  the  first  test,  the  second  test  gave  much  less  heat 
than  expected.  The  peak  flux  and  fluence  measured  by  4  calorimeters 
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TYPICAL  TRANSITS  TRS  SPECIMEN  MOUNT 


4.10.  Three  Inch  Diameter  Ceramic  Coatings  Test  Fixture 
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TYPICAL  TRS  INSTALLATION  -  ROUND  SPECIMEN 


TYPICAL  TRS  INSTALLATION  -  SQUARE  SPECIMEN 


Figure  4.12.  Test  Sample  Installation 


-47- 


j 


are  presented  in  Table  4  1.  The  calorimeters  at  the  5  feet 
distance  were  in  the  vicinity  of  the  samples,  well  within  the 
fireball.  The  calorimeter  located  at  30  feet  was  just  at  the 
edge  of  the  fireball.  The  calorimeter  output  plots  are  shown  in 
Figure  4  13.  The  calorimeter  voltage  outputs  shown  give  an  idea 
of  the  time  history  of  the  heating.  Calimeter  FX5  has  a  different 
sensitivity  than  FX1,  FX2,  and  FX3  have,  hence  the  voltages 
cannot  be  compared  directly. 


Table  4.1.  T?.S  Test  Calorimeter  .Measurements 


CALORIMETER 

NUMBER 

ENVIRONMENT 

CALORIMETER 

DISTANCE 

(FT) 

MAXIMUM  FLUX 
(CAL/CM2-SEC) 

FLUENCE 

(CAL/CM2) 

FX1 

FIREBALL 

5 

105 

45,6 

FX3 

FIREBALL 

5 

114 

41.3 

FX2 

RADIATION 
(BEHIND  GLASS) 

5 

43,9 

3.0 

FX5 

RADIATION, 

SOME  FIREBALL 

30 

68.2 

82,7 

The  fluence  received  by  the  calorimeter  behind  the  pyrex 
glass  shield  was  particularly  low.  The  pyrex  glass  (pyrex 
transmissivity  is  about  0.9)  is  almost  immediately  covered  by  an 
aluminum  oxide  film,  blocking  the  energy  transmitted  after  about 
0.15  seconds.  Thus,  the  glass  shield  proved  not  to  be  a  good 
way  to  prevent  the  samples  from  being  heated  by  convection/conden¬ 
sation  mechanisms.  Similarly,  the  samples  exposed  to  the  fireball 
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gases  become  very  contaminated  with  aluminum  oxide  and  carbon. 

Most  of  the  ceramic  samples  were  essentially  unaffected 
by  the  fluxes  obtained.  The  backface  temperature  rises  were  on 
the  order  of  5  to  80°F,  much  less  than  the  300  to  1000 °F  expected. 

The  specific  results  obtained  are  summarized  as  follows: 

1.  Two  3  inch  diameter  Al203-Ti02  coated  discs  showed 
thermal  stress  failure;  however,  the  4  inch  square 
samples  of  this  material  survived. 

2.  One  3  inch  diameter  A1203  coating  partially  separated 
from  the  steel  substrate. 

3.  The  alumina  coating  on  bulk  alumina  substrates 
spalled,  blistered,  and  cracked  on  all  3  samples 
tested. 

4.  Silicone  resin  infiltrant  formed  char  layers  on  the 
Dorous  ceramic  coating  surfaces. 

5.  Silastic  184  (silicone  resin  infiltrant  material)  and 
Sylgard  E  (silicone  rubber)  bulk  samples  developed 
char  layers  on  their  surfaces. 

6.  Indepth  thermocouples  indicated  radiation  transparency 
of  teflon,  at  least  to  a  depth  of  1  inch. 

Backface  temperatures  were  measured  on  the  ceramic  coated 
steel  samples.  However,  the  instrumentation  sensitivity  was 
sized  for  temperature  rises  on  the  order  of  300  to  1000°F,  and, 
with  the  small  amount  of  heat  received,  the  temperature  rises 
are  only  on  the  order  of  5  to  80°F.  These  temperatures  are  low 
enough  that  the  instrumentation  drift  greatly  affects  the  results. 
Because  the  backface  temperatures  are  of  no  practical  value  due 
to  the  low  heating  which  caused  questionable  accuracy,  they  are 
not  presented  here. 

All  samples  were  measured  and  weighed  before  and  after 
the  test.  However,  only  the  silicone  rubber  and  silicone  resin 
showed  any  significant  surface  decomposition.  Even  the  mass 
loss  in  this  case  was  less  than  what  could  be  accurately  deter¬ 
mined,  considering  the  alumina  contaminated  condition  of  the 
surface. 
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4.4 


French  Solar  Furnace  Test  Program 


A  series  of  tests  were  conducted  on  selected  material 
candidates  in  the  French  Solar  Furnace  at  Odeillo-Font-Romeu , 

France.  A  description  of  this  test  program  is  provided  on  the 
following  pages. 

4.4.1  Facility  Description 

The  solar  furnace  has  a  parabolic  reflector  with  a  focal 
length  of  59  feet,  is  130  feet  high  and  175  feet  wide,  and  is 
composed  of  9500  mirrors,  each  127  by  127  inches.  Because  the 
parabolic  reflector  is  too  large  to  track  the  sun,  63  smaller 
mirrors  (heliostats)  set  in  eight  tiers,  follow  the  sun  and 
reflect  its  rays  in  parallel  beams  onto  the  parabola.  The 
heliostats  are  each  24.6  by  19.7  feet. 

The  solar  energy  incident  on  an  area  of  about  23000  square 
feet  is  concentrated  by  the  parabolic  reflector  into  an  area 
about  2  feet  in  diameter.  The  diameter  of  the  image  of  the  sun 
at  the  focal  point  is  6.61  inches  and  27  percent  of  the  thermal 
energy  is  concentrated  in  this  area.  Reference  7  contains  a 
complete  description  of  the  distribution  of  energy  at  this 
point. 

The  SAI  samples  were  2  inch  diameter.  For  an  area  defined  by 

a  2  inch  diameter  at  the  focal  point,  the  minimum  design  heat 

2 

flux,  with  all  mirrors  operating,  is  365  +  15  cal/cm  sec  for  an 

2 

incident  radiation  of  959  watts/cm  ,  at  an  effective  radiation 
temperature  of  6840°F.  The  flux  is  supposed  to  vary  less  than  + 

2.5  percent  over  the  2  inch  diameter  area. 

In  this  test,  because  of  the  condition  of  the  mirrors, 

2 

the  flux  obtained  was  about  180  to  200  cal/cm  sec  and  varied 
less  than  +  5  percent  over  the  2  inch  diameter.  Exposure  times 
on  the  order  of  16  seconds  were  used  in  order  to  match  the  total 
fluence  for  the  nuclear  burst  model,  approximately  3100  cal/cm2. 
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4.4.2  Test  Procedure  and  Approach 

4. 4.2.1  Samples  Tested 

A  series  of  14  test  specimens,  consisting  of  8  different 
materials,  were  submitted.  These  samples  are  listed  in  Table  4. 

2.  Six  of  the  samples  were  tested,  as  noted  in  Table  4.2.  All 
samples  were  2  inch  diameter  disks. 

4 . 4 . 2 . 2  Instrumentation 

Samples  were  instrumented  with  chromel-alumel  thermocouples 
on  the  backface.  Savereisen  cement  was  used  to  fasten  thermo¬ 
couples  to  the  non  metallic  specimens  and  in  some  cases  the  bond 
failed  due  to  insufficient  cure  time  before  shipping  to  France. 
Metal  substrates  had  thermocouples  spot  welded  to  the  rear  face. 
Front  face  temperatures  were  obtained  using  an  optical  pyrometer, 
which  operated  at  a  wavelength  at  which  the  surface  emissivity 
is  approximately  1,  according  to  Reference  8.  The  pyrometer 
field  of  view  covered  a  1  inch  diameter  area  in  the  center  of 
the  2.0  inch  diameter  sample  face. 

Motion  pictures  were  taken  of  the  specimen  front  surface 
during  the  test  with  a  20  frame/second  camera.  All  samples  were 
measured,  weighed,  and  photographed  before  and  after  the  test. 

4.4.3  Pre-Test  Analysis 

The  thermal  response  of  the  alumina  coating  with  the 
silicone  infiltrant  exposed  to  the  heat  flux  of  the  solar  furnace 
was  predicted  using  the  thermal  model  described  in  Section  4.1. 

The  physical  parameters  are  as  follows: 

2 

Flux  =  350  cal/cm  sec 

Exposure  Time  =  9  seconds 

2 

Fluence  =  3150  cal/cm 
Alumina  Thickness  =  0.125  in 
Steel  Thickness  =  0.25  in. 

Porosity  =  15% 

Alumina  Absorptivity  =  0.2 
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The  resulting  temperature  and  recession  predictions  are  shown  in 
Figure  4.14.  Melting  of  the  surface  begins  at  6.7  seconds  and 
the  total  predicted  recession  is  only  .007  inches. 

In  the  actual  tests,  discussed  in  the  following  paragraphs, 

2 

the  flux  was  195  cal/cm  sec  for  16  seconds.  Also,  the  absorp¬ 
tivity  appeared  to  be  much  larger,  so  that  3  to  4  times  as  much 
heat  would  be  absorbed  as  used  in  the  above  prediction.  Therefore, 
the  prediction  model  was  updated  using  these  data.  Absorptivities 
of  0.2  and  0.5  were  used.  Even  an  absorptivity  of  0.5  is  probably 
less  than  what  it  actually  was  in  the  test. 

The  resulting  temperature  predictions  are  shown  in  Figure 
4.15  and  the  recession  predictions  are  presented  in  Figure  4.16. 
With  an  absorptivity  of  0.2,  the  melt  temperature  is  not  obtained 
and  there  is  no  recession.  However,  with  an  absorptivity  of 
0.5,  the  melt  temperature  is  obtained  at  the  surface  in  4  seconds, 
and  at  16  seconds  the  total  recession  is  .068  inches,  which  is 
more  than  half  of  the  material. 

4.4.4  Test  Results 

A  description  of  the  material  surface  effects  seen  on  the 
film  are  shown  in  Table  4.3.  The  measured  front  face  and  back 
face  temperatures,  pre-test  and  post-test  pictures  and  sample 
dimension  and  mass  data  are  presented  in  Appendix  H.  A  description 
of  the  response  of  each  material  is  given  in  the  following 
paragraphs. 

4. 4. 4.1  Alumina/Silicone  Inf iltrant/Steel 

A  dark  brown  char  layer,  caused  by  the  silicone  resin 
charring,  develops  immediately,  increasing  the  surface  absorp¬ 
tivity.  Also,  the  alumina  coating  was  contaminated  (dirty 
appearance)  due  to  the  processing.  These  factors  caused  the 
heat  absorbed  to  probably  be  3  to  4  times  that  anticipated  for 
the  expected  absorptivity  of  0.2.  Melting  began  at  3.5  seconds 
and  spalling  of  the  coating  begins  at  4.5  seconds.  Very  high 
surface  temperatures  were  measured  (Appendix  H)  about  2500  to 
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Figure  4,14.  Predicted  Thermal  Response  of  Alumina  Coating  With 
Silicone  Infiltrant  in  Solar  Furnace 
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Effect  of  Absorptivity  on  Predicted  Aluminum  Temperatures 


Table  4.3.  Material  Visual  Response 


Sample 

and 

Flux 

Exposure 

Time 

(Sec) 

Sample  Response 

F2  -  Alumina/Sili¬ 
cone/Steel 

0.1 

Brown  char  develops,  smokes 

1.0 

Char  layer  lightens 

Flux  =2197.5 
cal/cm  -sec 

Fluence  -  3152 
cal/cin 

3.5 

Alumina  melting  begins 

: 

4.5 

Upper  1/3  of  alumina  spalls 
off 

13.1 

Steel  melt  begins 

14.5 

Most  of  remaining  lower  part 
of  alumina  spalls  off 

16.0 

End 

F3  -  Silastic  E/Steel 

0.15 

Char  fully  formed 

Flux  --191.3 
cal/cm  -sec 

2.2 

Outgassing,  smoking  visible 

Fluency  -  3061 
cal/cm 

.  16.0 

End 

F4  -  Silastic  E/Steel 

No  film  data 

Flux  =2187.3 
cal/cm  -sec 

Fluence  =  2992 
cal/cin 

F8  -  Boron  Nitride, 
Grade  HBR 

9.1 

Sides  begin  to  spall 

9.6 

Dark  flaw  appears  in  center, 
outgassing  visible 

Flux  =2182.9 
cal/cm  -sec 

17.0 

End 

Fluency  =  2962 
cal/cm 


Table  4,3,  Material  Visual  Response  (Cont'd) 


Sample 

and 

Flux 

Exposure 

Time 

(Sec) 

Sample  Response 

F9  -  AS-3DX 

1.0 

Outgassing  begins 

Flux  =2184.8 
cal/cm  -sec 

10.5 

Small  bubble  forms 

Fluence  =  2957 
cal/cnr 

Large  bubbles,  discoloration 

IS 

End 

F16  -  Spinel/Sili¬ 
cone/Steel 

0.1 

Brown  char  layer  develops, 
smokes 

Flux  =.193.6 
cal/cm  -sec 

0.8 

1/8  of  coating  spalls  off 

Fluency  =  3098 
cal/crn 

1.0 

Char  lightens 

1.2 

^  1/3  of  coating  spalls  off 

1.8 

Remaining  coating  is  clean 
looking 

Coating  melting  in  center 

^  1/4  of  coating  melts,  falls, 
expases  center 

3/4  of  coating  gone 


Steel  begins  to  melt 
End 


4.9 

8.6 

16.0 


2800°C,  which  is  above  the  alumina  melt  temperature  (2100°C) . 

This  may  be  the  silicone  char  temperature.  Backface  temperatures 
were  not  measured  on  this  sample. 

4. 4. 4. 2  Silastic  E  on  Steel 

The  silicone  rubber  Silastic  E  samples  performed  very 
well  during  the  test.  A  char  layer  forms  immediately  and  protects 
the  remaining  material  by  both  insulating  the  virgin  material 
and  reradiating  at  a  high  temperature  (3000°C,  see  Appendix  G) . 

The  decomposing  hydrocarbons  in  the  silicone  rubber  and  the 
pyrolysis  gases  passing  through  the  char  layer  remove  additional 
heat.  As  a  result,  a  .07  to  .09  inch  thick  char  layer  formed 
and  the  steel  backface  temperature  increases  were  only  180°C  and 
50°C  for  the  0.25  inch  and  0.5  inch  thick  samples  respectively. 
There  was  some  dark  brown  smoke  emitted  which  left  a  small 
deposit  on  the  holder,  and  a  considerable  amount  of  whitish  gas 
was  emitted.  A  chemical  analysis  of  the  post-test  sample 
showed  that  the  char  layer  was  7.0  percent  carbon. 

4. 4. 4. 3  Boron  Nitride.  Grade  HBR 

The  boron  nitride  began  to  spall  at  7.1  seconds  and  had 
severe  cracks  at  the  end  of  exposure.  The  frontface  temperature 
profile  is  shown  in  Appendix  G.  The  temperature  profile  is 
quite  non-uniform  which  can  be  due  to  either  changing  emissivity 
in  the  pyrometer  wavelength,  changing  absorptivity,  and/or 
conductivity  decreasing  with  temperature.  Backface  temperature 
data  was  not  available. 

4. 4. 4. 4  AS-3DX 

The  AS-3DX  exhibited  outgassing  immediately;  however,  the 
surface  did  not  exhibit  deterioration  until  after  10  seconds  of 
exposure.  The  surface  temperature  is  shown  in  Appendix  G.  The 
surface  temperature  rise  rate  indicates  a  very  low  absorptivity. 

A  non-decomposing  material  thermal  model  was  made  of  the  AS-3DX 
and  the  results  are  shown  in  Figure  4.17.  A  surface  absorptivity 
of  0.1  matches  the  experimental  surface  temperatures  quite  well. 

An  absorptivity  of  0.2  was  also  considered,  and  the  melt  temper- 
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gure  4,] 7.  AS-3DX  Front  Face  Temperature 


ature,  about  3000°F,  is  reached  in  about  2  seconds.  This  is 
very  similar  to  what  would  happen  if  the  heat  flux  were  doubled. 
Thus,  melting  and  bubbling  would  probably  occur  immediately  in 
the  nuclear  thermal  environment;  however,  the  melting  would  be 
over  a  shorter  time  span. 

The  backface  temperature  in  Appendix  G  for  sample  F9  may 
not  be  accurate.  There  was  difficulty  in  establishing  the  time 
scale  and  also,  the  thermocouple  came  loose  at  some  time  that  is 
not  known.  If  the  temperatures  are  right,  they  are  much  higher 
than  expected  in  this  highly  insulative  material,  indicating 
that  AS-3DX  may  be  transparent  to  the  radiation  energy. 


/ 
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4. 4. 4.1  Spinel/Silicone  Infiltrate/Steel 

This  material  behaves  similarly  to  the  alumina/infil¬ 
trate  sample.  The  infiltrare  charred  and  increased  the 
absorptivity  immdiately.  Also,  spalling  began  at  an  exposure 
of  less  than  0.8  seconds.  Melting  is  evident  at  2.9  seconds, 
and  three-fourths  of  the  coating  has  spalled  off  at  4.9 
seconds . 

4.5  Conclusions 

Material  response  thermal  models  of  the  porous  cer¬ 
amic  coatings,  with  and  without  silicone  resin  infiltrant, 
were  developed.  Similar  models  for  silicone  rubber  were 
partially  developed.  All  of  the  materials  considered  were 
subjected  to  the  Thermal  Radiation  Simulator  tests,  and  a 
limited  number  of  materials  were  subjected  to  the  French 
Solar  Furnace  test.  The  conclusions  obtained  are  as  follows: 

1.  Test  Facilities 

The  large  scale  TRS  test  facility  did  not  produce 
the  flux  and  fluence  needed.  For  even  moderately  high  heat¬ 
ing  rates,  the  sample  must  be  in  the  fireball,  which  con¬ 
taminates  the  surface  with  aluminum  oxide  and  carbon  residue. 
Further  calibration  of  the  TRS  fireball  heat  flux  is  needed 
as  well  as  a  better  understanding  of  the  kinetics  and  dynamics 
of  the  fireball.  It  currently  could  be  used  for  testing  at 
lower  fluxes  modeling  lesser  nuclear  environments;  but,  it 
apparently  does  not  have  the  capability  to  give  the  environ¬ 
ment  conditions  required  in  this  study. 

The  FSF  does  supply  sufficient  flux  and  fluence  to 
severly  test  the  materials.  However,  since  its  flux  is 
much  less  than  the  nuclear  environment  being  modeled  and 
long  heating  time  durations  are  required,  the  material  re- 
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sponse  is  different  than  what  would  be  expected  in  a  nuclear 
environment.  However,  it  is  the  best  testing  method  avail¬ 
able. 


2.  Porous  Ceramic  Coatings 

Based  on  the  thermal  tests,  primarily  the  FSF  test, 
the  porous  ceramic  coatings  conclusions  are: 

a)  The  silicone  resin  infiltrant  forms  a  thin 
char  layer  which  increases  the  absorptivity 
of  the  ceramic  coating.  A  low  absorptivity 
of  about  0.2  is  required,  hence,  the  infil¬ 
trant  should  not  be  used. 

b)  Process  methods  must  be  carefully  controlled 
to  insure  a  non-contaminated,  low  absorptivity 
coating. 

c)  Ceramic  coatings  are  very  susceptable  to 
thermal  stress.  Pre-stress  durinq  processing, 
material  expansion  compatibility,  etc.  should 
be  investigated  to  alleviate  this  problem. 

The  thermal  modeling  gave  the  following  significant  conclu¬ 
sions  : 

a)  Alumina  coating,  0.125  inches  thick,  can  survive 
three  nuclear  bursts,  if  it  has  an  absorptivity 
on  the  order  of  0.2.  The  other  ceramics  may 
have  similar  performance,  but  only  the  alumina 
was  investigated.  This  is  based  on  material 
loss  due  to  melting  only,  thermal  stress  is  not 
considered. 

b)  A  thick  coating  exhibits  more  ablation  than  a 
thin  one;  however,  the  thicker  coating  also 
has  more  material  remaining  after  the  exposure. 

3.  Silicone  Rubber 


The  conclusions  obtained  about  using  silicone  rubber 
as  an  ablative  thermal  protection  coating  in  the  nuclear  en¬ 
vironment  are  as  follows : 
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a)  A  highly  insulative,  low  carbon  content 
char  layer  develops  and  protects  the  re¬ 
maining  virgin  material  and  structure 
from  severe  heating.  -This  material  sur¬ 
vived  the  3150  cal/cni  fluence  of  the 
FSF  test  with  very  little  material  loss 
and  low  backface  temperature  rise. 
Different  types  of  silicone  rubber  have 
different  types  of  chars,  ranging  from 
very  fragile  to  very  tough. 

b)  The  thermal  protection  capability  is  in¬ 
dependent  of  contaminated  surface  condi¬ 
tions  . 

c)  It  is  a  very  practical  material,  in  that 
it  is  inexpensive,  easy  to  apply,  adheres 
well,  resists  weathering,  and  acts  as  a 
moisture  barrier. 

4.  Miscellaneous  Bulk  Materials 


Results  of  the  thermal  tests  produced  the  following 
pertinent  results  on  a  few  of  the  materials  tested: 

a)  Teflon  is  semi-transparent  to  thermal 
radiation  in  the  short  IR  wavelengths 
being  investigated. 

b)  AS-3DX  exhibited  very  low  absorptivity  and 
held  up  well  at  the  test  flux  of  185  cal/cnT 
sec. 

c)  Union  Carbide  boron  nitride  HBR  showed  severe 
spalling  and  cracking  in  the  FSF  tests . 


■ 
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5.0  ELECTRICAL  PROPERTIES  TEST  AND  ANALYSIS 

5.1  Introduction  And  Summary 

This  chapter  presents  the  electromagnetic  aspects  of  the 
array  hardening  problem.  The  electrical  properties  of  the  ma¬ 
terials  are  tabulated  and  discussed,  while  the  array  implications 
of  the  use  of  these  materials  are  briefly  considered. 

The  electrical  property  values  of  the  candidate  materials 
were  determined  by  the  Engineering  Experiment  Station  at  the 

Georgia  Institute  of  Technology  and  are  presented  later  in  this 

(9) 

section.  The  complete  data  base  from  which  these  average 

values  were  calculated  are  given  in  Appendix  B.  A  brief  dis¬ 
cussion  of  the  theoretical  basis  for  the  cavity  measurement  pro¬ 
cedure  used  for  the  plasma-sprayed  ceramic  materials  is  included 
later  in  this  chapter.  The  complete  analysis  which  supports  this 
discussion  is  given  in  Appendix  C.  The  conclusions  of  the  array- 
material  interaction  question  are  summarized  in  Section  5.1,  while 
the  third  section  of  this  chapter  considers  the  problem  in  some¬ 
what  more  detail.  The  theoretical  analysis  associated  with  the 
interaction  question  can  be  found  in  Appendix  D. 

5.1.1  Plasma-Sprayed  Ceramic  Coatings 

The  dielectric  properties  of  the  ceramic  coatings  are 
summarized  in  Tables  5.1  and  5.2.  These  two  tables  present  the 
results  of  S-band  and  X-band  measurements,  respectively.  For 
each  material,  an  average  (Avg.)  and  a  standard  deviation  (a) 
are  shown  for  the  dielectric  constant  (er)  and  loss  tangent 
(tan  6) .  The  statistics  are  based  on  repeated  test  runs  for 
each  sample,  with  the  cavity  completely  disassembled  between 
runs.  The  S-band  tests  were  repeated  at  least  six  times,  and 
the  X-band  tests  of  unburned  samples  were  repeated  four  times. 
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These  samples  are  believed  to  more  porous  than  the  other  samples. 


Note  that  these  runs  are  all  for  the  same  sample.  Thus  the 
statistics  reflect  the  consistency  of  the  measurement  process 
rather  than  the  material  variability  from  sample  to  sample.  A 
complete  tabulation  of  the  results  of  individual  tests  is  pro¬ 
vided  in  Appendix  B.  The  statistical  averages  are  developed 
after  discarding  the  highest  and  lowest  values  of  e and  tan  6. 
The  reason  for  repeated  tests  and  for  discarding  extreme  data 
is  that  the  particular  (TM)  cavity  mode  excited  in  these  tests 
is  known  to  be  sensitive  to  the  electrical  contact  at  the 
assembly  joint  in  the  cavity  wall.  An  alternative  (TE)  mode 
would  have  been  less  sensitive  to  the  conduction  properties  of 
the  joint,  but  it  may  not  have  provided  adequate  measurement 
sensitivity  for  these  electrically  thin  coatings. 

Three  samples  of  burned  coatings  were  provided  for 
X-band  measurements.  The  burns  were  produced  with  an  arc  jet. 
The  burned  samples  had  apparently  expanded  and  were  difficult 
to  insert  in  the  cavity.  One  sample,  #356,  shattered  on  the 
first  attempt  to  fit  it  into  the  cavity,  and  the  other  two 
samples,  #202  and  #333,  cracked  when  they  were  being  pressed 
out  of  the  cavity  after  second  test  runs.  Therefore,  the 
variance  of  the  data  cannot  be  assessed. 

The  unburned  ceramics  generally  have  a  speckled 
appearance.  Contamination  by  metallic  particles  would  con¬ 
tribute  an  "artificial  dielectric"  mechanism  and  could  explain 
the  relatively  high  dielectric  constants  that  were  measured 
for  some  ceramics.  Samples  #2  through  #6  are  known  to  be  less 
porous  than  the  others,  and  this  fact  is  reflected  in  their 
higher  dielectric  constants.  For  the  two  materials  that  were 
tested  before  and  after  burning,  the  burn  increased  the 
dielectric  constant  and  the  loss  tangent  significantly.  The 
surfaces  of  the  burned  samples  were  glazed,  and  it  seems  likely 
that  some  change  of  state  occurred.  Because  the  burned  samples 


were  not  perfectly  homogeneous,  it  is  possible  that  the 
theoretical  model  used  for  data  analysis  may  have  been  in¬ 
adequate,  but  the  results  are  believed  to  be  reasonably 
accurate. 

5.1.2  Bulk  Materials 

The  dielectric  properties  of  the  bulk  materials  are 
summarized  in  Tables  5.3  and  5-4,  which  present  the  results 
of  S-band  and  X-band  measurements,  respectively.  The  solid, 
silicone-based  materials  were  molded  to  fill  short  lengths  of 
rectangular  waveguide,  and  shorted  waveguide  sections  were 
used  as  sample  holders.  For  the  ash  materials,  the  sample 
holder  was  a  female  SMA  coaxial  connector,  from  which  the 
dielectric  plug  had  been  removed. 

The  two  basic  solid  materials  that  were  tested  are 
distinguished  by  color  (i.e.,  clear  RTV  and  white  RTV) .  Two 
sample  lengths,  denoted  -1  and  -2,  were  tested  at  each  of  the 
two  nominal  frequencies,  3  GHz  and  9  GHz.  The  white  material 
was  retested,  after  one  end  of  the  samples  had  been  burned. 

The  X-band  samples  were  burned  with  an  arc  jet,  while  the 
S-band  samples  were  burned  with  an  oxyacetylene  torch.  The  jet 
produced  smooth,  hardened  surfaces,  with  a  fine,  black,  ash 
residue.  The  torch  produced  rough,  hardened  surfaces,  covered 
with  a  black-speckled,  loose,  granular  ash.  The  arc  jet  ash 
was  mostly  lost  during  the  burning  process.  The  loose  ash  from 
the  oxyacetylene  torch  was  brushed  away  and  saved  before  tests 
were  run  on  the  burned  solids.  The  burned  surfaces  were  not 
uniform,  so  measurement  of  the  effective  sample  thicknesses 
involved  some  exercise  of  judgment. 

The  sample  holders  used  for  the  bulk  tests  did  not 
have  mechanical  joints  that  seriously  affected  the  electrical 
repeatability  of  the  measurements.  Therefore,  the  tabulated 
values  are  the  results  of  single  measurements,  rather  than 
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statistical  averages  of  repeated  measurements.  The  solid 
materials  were  supplied  in  two  sample  lengths  for  each  wave¬ 
guide  band.  The  nearly  equal  values  obtained  for  the  two  sample 
lengths  are  a  good  indication  of  the  repeatability  of  the  measure¬ 
ments  . 

Three  types  of  ash  were  tested.  The  granular  ash  that 
resulted  from  the  torch  burn  on  the  white,  S-band  samples  was 
tested  at  S-band.  The  fine  black  ash  that  resulted  from  the 
arc  jet  burn  on  the  white,  X-band  samples  was  tested  at  X-band. 

Ash  from  a  sample  of  white  RTV  that  was  burned  in  the  French 
solar  furnace  was  tested  in  both  frequency  bands.  In  all  cases, 
the  ash  was  packed  as  tightly  as  possible  into  a  short  (0.260 
in.)  length  of  short-circuited,  coaxial  transmission  line.  The 
electrical  properties  of  powders  vary  substantially  with 
density,  although  empirical  evidence indicates  that  the 
dielectric  constant  is  more  sensitive  to  density  than  is  the 
loss  tangent.  The  packing  density  could  not  be  carefully  con¬ 
trolled  in  these  tests. 

5.1.3  Material-Array  Interactions 

The  hardening  concepts  under  consideration  here 
involve  the  use  of  waveguide  plugs  and/or  slabs  which  cover 
the  antenna  array.  The  effects  of  these  structures  on  array 
performance  are  well  known  to  be  deleterious,  and  involve  the 
issues  of  (1)  reflections,  (2)  losses  and  (3)  blind  spots. 

For  hardening  applications  (or  any  other  that  requires 
covering  slabs  or  waveguide  plugs) ,  it  is  desirable  that  any 
material  have  as  low  a  dielectric  constant  and  loss  tangent  as 
possible.  Further,  in  order  to  reduce  the  incidence  of  array 
blindness,  it  is  exceedingly  desirable  that  the  material  be 
as  thin  as  possible.  From  the  interaction  point  of  view,  the 
array  blindness  is  perhaps  the  most  significant  problem  for 
either  the  plasma- sprayed  ceramic  materials  or  the  silicone 
rubber  (RTV)  materials. 
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In  general,  the  most  significant  conclusion  of  the 
interaction  analysis  is  that  the  hardening  of  the  array  must 
be  treated  as  an  integral  part  of  the  array  design  process. 

Given  an  integral  design,  there  is  reasonable  hope  that  the 
silicone  rubber  materials  as  well  as  some  of  the  plasma— sprayed 
materials  may  have  some  applicability  to  the  array  hardening 
problem.  In  particular,  Al2°3  and  Spinel  are  the  most  likely 
ceramic  candidates  from  the  electrical  point  of  view.  Although 
the  silicone  rubber  materials  produce  a  char,  the  loss  associated 
with  the  char  appears  to  be  tolerable  provided  that  the  char  is 
not  too  thick.  For  most  of  the  char  values  in  the  tables,  this 
means  less  than  0.125  inches. 

5.2  Electrical  Properties  Test  Program 

This  section  describes  the  electrical  properties  test 
program.  It  consists  of  three  parts:  (1)  Test  Facility 
Description,  (2)  Test  Procedures  and  Approach,  and  (3)  Pre- 
Test  Analysis.  The  first  two  parts  deal  primarily  with  the 
resonant  cavity  measurement  procedure.  The  bulk  sample  measure¬ 
ment  technique  is  conventional,  and  is  not  discussed  in  detail. 

5.2.1  Test  Facility  Description  (Resonant  Cavity  Measurements) 

The  dielectric  properties  of  the  thin  ceramic  coatings 
were  measured  using  a  resonant  cavity.  The  cavity  was  a  right- 
circular  cylinder  of  radius  a  and  height  d,  with  d<<a.  The 
coatings  were  applied  to  cover  the  surface  of  one  end  of  the 
cavity,  and  the  electrical  coupling  was  made  via  a  short  coaxial 
probe  inserted  through  the  opposite  end  of  the  cavity. 

A  cross-section  of  the  cavity  is  shown  in  Figure  5.1A.  A 

photograph  of  the  disassembled  cavity  is  shown  in  Figure  5.1B 

and  a  photograph  of  the  assembled  cavity  provided  as  Figure  5. 1C. 

The  main  box  (the  upper  section  in  Figure  5. 1A  is  brass,  and  the 

interchangable,  ceramic-coated  endplates  are  stainless  steel.  The 
endplate  is  held  tightly  in  the  cavity  by  a  clamping  plate, 
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7MM  Precision 
Connector 


Figure  5.1A.  Sectional  View  of  the  Resonant  Cavity  Sample  Holder  to  be 
Used  for  Testing  Thin  Coatings 


PRECISION  CONNECTOR 


Figure  5. 1C,  Photograph  of  Assembled  Cavity 


which  is  bolted  through  the  main  body.  The  nominal  dimensions 
of  the  three  cavities  that  were  used  for  these  tests  are  shown 
in  Table  5.5.  The  nominal  coating  thicknesses  were  in  the  range 
0.050  to  0.100  in. 

The  coupling  probe  was  the  center  contact  extension  of 
a  flange -mounted,  7mm  precision  connector.  This  configuration 
minimized  transmission  line  discontinuities  between  the 
cavity  port  and  the  instrumentation.  The  probe  penetration 
was  0.060  in.  for  S-band  measurements  and  0.040  in.  for  X-band 
measurements.  The  probe  diameter  was  0.050  in. 

The  basic  instrumentation  used  in  this  measurement  is 
diagrammed  in  Figure  5.2A  and  photographically  illustrated  in 
Figure  5.2B.  A  list  of  equipment  is  given  in  Table  5~6.  The 
Network  Analyzer  and  the  Reflection/Transmission  Test  Set  provide 
the  capability  for  swept- frequency  measurements  of  complex  reflec¬ 
tion  coefficients.  Amplitude  and  phase  are  plotted  vs.  frequency 
on  an  X-Y  recorder.  The  system  is  calibrated  by  recording  the 

(11/ 

reflection  coefficients  of  several  (at  least  three)  known  loads 
For  these  measurements,  calibration  loads  included  a  matched  ter¬ 
mination,  an  open  circuit,  and  short  circuits  placed  at  three 
different  positions  relative  to  the  connector  reference  plane. 


5.2.2  Test  Procedure 


5. 2.1.1  Resonant  Cavity  Measurements 


Before  each  experiment,  the  equipment  is  allowed  to 
warm  up  for  one  hour.  The  sample  to  be  tested  and  the  brass 
cavity  are  cleaned  with  methanol  and  dried  with  laboratory 
tissues.  The  sample  and  end-plate  are  placed  in  position  and 
the  cavity  body  is  assembled.  Caution  must  be  taken  at  this 
point  to  torque  the  bolts  holding  the  cavity  to  a  uniform 
tightness  since  the  TMq^q  cavity  mode  is  known  to  be  sensitive 
to  the  electrical  conductivity  of  the  mechanical  joint  between 
the  end-plate  and  the  cylinder.  A  problem  of  this  type  would 
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Table  5,6.  Test  Equipment  for  Cavity  Measurements 


Item 

Manuf ac  turer/Model 

Sweep  Oscillator 

Hewlett  Packard/HP-8620C 
*We in  sc  he 1 / 4 30-A 

Frequency  Meter: 

Counter 

Wavemeter 

Hewlett  Packard/HP-5245 
♦Hewlett  Packard/HP-5 37A 

Reflection/Transmission 

Test  Sec 

Hewlett  Packard/KP-8743A 

Network  Analyzer 

Hewlett  Packard/HP-S410B 

Phase /Amplitude  Display 

Hewlett  Packard/HP-8412A 

Recorder 

Hewlett  Packard/HP- 7044A 

*  Alternate  Equipment 
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be  evidenced  by  poor  repeatability  of  the  test.  To  minimize 
this  problem,  each  sample  was  tested  several  times,  and  extreme 
(high  and  low)  values  were  discarded. 

A  Reflection/Transmission  Test  Set  is  used  to  separate 
incident  and  reflected  signals  at  the  cavity  excitation  port, 
and  a  Network  Analyzer  is  used  to  measure  the  ratio  of 
ref lected-to- incident  signals.  The  "Gain"  and  "Phase  Offset" 
controls  of  the  Network  Analyzer  provide  the  basic  amplitude 
and  phase  calibration  information.  The  Test  Set  can  be 
viewed  as  an  imperfect  two-port  network  that  is  placed  between 
the  Network  Analyzer  and  the  load  under  test.  The  system 
calibration  can  be  improved  by  determining  the  scattering 
parameters  of  the  two-port.  Such  a  calibration  is  accomplished 
routinely  with  a  computerized  "Automatic  Network  Analyzer",  and 
the  same  procedures  have  been  adapted  for  these  cavity  measure¬ 
ments  with  the  manual  Network  Analyzer . ^  Basically,  the 
scattering  parameters  of  the  two-port  are  determined  by  recording 
reflection  coefficients  for  several  known  loads.  Typically,  a 
matched  termination,  a  short,  and  an  open  are  used.  Offset 
shorts  may  be  used  instead  of  the  open.  The  procedure  adopted 
for  these  tests  uses  a  matched  termination,  three  shorts,  and 
one  open.  The  three  shorts  are  designated  S-0  (Offset=0 . 000  in), 
S-l  (Of fset=0. 289  in),  and  S-2  (Of f set*0 . 983  in).  The  open 
circuit  was  represented  mathematically  by  the  capacitance  of 
an  open-circuited  7mm  connector.  This  calibration  procedure 
requires  manual  transcription  (digitization)  of  a  large  quantity 
of  recorded  data,  but  a  substantial  improvement  in  measurement 
accuracy  is  achieved. 

A  Frequency  counter  is  used  to  measure  frequencies  at 
discrete  points  in  the  sweep  range  for  each  test  recording .  For 
the  relatively  narrow  sweep  ranges  used  here,  the  microwave 
frequency  varies  linearly  with  the  dc  "Sweep  Output"  signal  from 


the  Sweep  Oscillator.  Therefore,  the  counter  is  used  to 
measure  the  frequency  at  the  start  of  the  sweep  and  at  the 
resonance  peak,  and  intermediate  frequencies  are  determined 
by  linear  interpolation. 

The  Network  Analyzer  measures  the  cavity  reflection 
coefficient  at  the  7mm  connector  interface  on  the  Test  Set. 

The  data  is  digitized,  and  a  computer  program  is  used  to  locate 
the  reference  plane  that  corresponds  to  the  plane  of  the  detuned 
open.  At  the  detuned  open  position,  the  impedance  varies 
like  that  of  a  series  RLC  circuit, 

Z  =  R  +  j  (wL  -  -^?)  (5-1) 

a)  C 

where  u  is  the  source  frequency.  For  real  frequencies,  the 
impedance  locus  is  a  straight  line  that  is  parallel  to  the 
imaginary  axis.  At  other  reference  planes,  the  locus  is  not 
parallel  to  the  imaginary  axis.  The  computer  program  adjusts 
the  reference  plane  to  minimize  the  standard  deviation  of  the 
resistance, 

°R  =  (Rav  ■  Rn)2]  (5"2) 

n=l 

where  R  is  the  resistance  at  the  n-th  frequency,  and  R  is 
n 

the  average  resistance.  When  a  minimum  is  obtained,  the  detuned 
open  position  is  located.  The  plane  of  the  detuned  open  will 
nearly  coincide  with  the  wall  where  the  probe  enters  the  cavity. 

At  the  resonant  frequency 

ui1  =  (LC)_1/2  (5-3) 
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the  locus  intersects  the  real  axis  at  Z=R.  At  the  half-power 
frequencies  u>^+,  the  impedances  are 

Z  (0^+)  -  R(1  +  j)  (5-4) 

Therefore,  when  the  cavity  impedance  is  referred  to  the  plane 
of  the  detuned  open,  the  resonant  frequency  and  the  half -power 
bandwidth  may  be  easily  recognized.  Near  resonance,  the 
reactance  varies  almost  linearly  with  frequency, 

X (at)  =  s (u>  -  u^)  (5-5) 

where  s  is  the  slope  of  the  curve  at  resonance.  The  slope  is 
determined  by  a  least-squares  fit  of  a  linear  function  to  the 
measured  data.  The  half-power  frequencies  are  calculated  in 
terms  of  the  resistance  at  resonance  and  the  slope, 

0)^  +  =  +  R/s  (5—6) 

The  factor  is  a  measure  of  the  cavity  resonance 

quality,  with  effects  of  the  coupling  network  removed.  In 

general,  includes  losses  in  the  cavity  walls  and  losses  in 

any  dielectric  material  in  the  cavity.  If  a  factor  Qw  is 

associated  with  wall  losses,  and  a  factor  is  associated  with 

°  (5) 

dielectric  losses,  the  unloaded  factor  is  given  by  , 


If  a  factor  Qw'  is  determined  from  measurements  of  an  empty 
cavity,  Qw  in  Equation  (5-7)  is  related  to  Qw'  by  the  ratio  of 
resonant  frequencies. 
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15-8) 


°W  = 


where  u)1 '  is  the  resonant  frequency  of  the  empty  cavity.  When 
low-loss  dielectrics  are  measured,  wall  losses  are  a  significant 
contribution  to  the  unloaded  Q,  and  the  measurement  of  an  empty 
cavity  is  necessary  for  an  accurate  calculation  of  Qq.  Because 
this  procedure  incorporates  all  of  the  recorded  data,  it  is 
presumably  less  sub'ject  to  errors  than  one  based  on  a  simple 
observation  or  interpolation  of  the  half-power  points.  However 

It 

the  assumption  that  X(u)  is  linearly  varying  does  imply  that 
the  half-power  bandwidth  must  be  relatively  narrow  and, 
therefore,  that  the  losses  must  be  small.  This  assumption  is 
not  restrictive  for  the  materials  tested  here. 


In  these  tests,  the  walls  of  the  main  body  of  the  brass 
cavity  were  polished  to  reduce  surface  losses.  However,  the 
ceramics  were  sprayed  on  relatively  rough,  molybdenum-coated 
end-plates.  To  determine  Qw',  empty  cavity  measurements  were 
made  using  a  molybdenum-coated  end-plate  with  no  ceramic  coating 


The  objective  of  the  measurement  is  to  determine  the 
complex  wave  number  kQ.  In  the  complex  frequency  plane, 
u>  =  w  +  jv,  and  the  resonance  occurs  at 

C 


ai 


=  k  (y  e  ) 
o  o  o 


-1/2 


(5-9) 


It  can  be  shown  that  uQ  can  be  expressed  in  terms  of  a  real 
frequency  <u  and  the  quality  factor  at  u>^,  as  follows: 


(5-10) 
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The  cavity  reactance  is  zero  at  and  Q1  is  simply  related 
to  the  half-power  bandwidth.  Therefore,  ,  Q^,  u)Q,  and  kQ 
may  be  determined  from  a  swept-frequency  impedance  measurement. 

Table  5.7  shows  dielectric  properties  for  some  known 
materials.  The  measured  values  were  obtained  using  the 
measurement  and  data  analysis  procedures  that  were  developed 
for  these  tests.  The  results  agree  reasonably  well  with  the 
published  values. 

5. 2. 2. 2  Shorted  Waveguide  Measurements 

The  dielectric  parameters  of  a  solid  dielectric  can 
be  determined  from  a  measurement  of  the  input  impedance  of  a 
sample— filled  section  of  shorted  transmission  line.  The 
measurement  is  convenient  because  special  sample  holders  are 
not  necessarily  required  and  because  the  measurement  itself 
can  be  accomplished  by  conventional  microwave  impedance 
measurement  techniques  (e.g.,  with  slotted  lines,  reflectom- 
eters,  or  microwave  network  analyzers) .  The  basic  procedures 
are  outlined  here  and  measurement  results  for  some  known 
dielectrics  are  given  in  Table  5. 3 . 

A  sample  holder  for  this  measurement  can  be  fabricated 
from  standard,  uniform  transmission  line  (e.g.,  rectangular 
waveguide  or  coaxial  line).  As  shown  in  Figure  5.3,  a  section 
of  the  line  is  shorted  at  one  end,  and  the  end  adjacent  to  the 
short  is  completely  filled  with  the  dielectric  to  be  tested. 

The  thickness  t  of  the  dielectric  sample  should  be  measured 
with  an  accuracy  better  than  one  percent  of  the  wavelength  in 
the  sample.  The  standard  transmission  line  is  usually  air- 
filled,  but  in  general  it  may  be  filled  with  a  non-magnetic , 
low-loss  dielectric  having  permittivity  e^.  The  free  space 
permeability  and  permittivity  are  yQ  and  eq,  respectively, 


Table  5.7.  Nominal  Dielectric  Parameters  of  Known  Materials 


Material 

Source* 

Frequency 

(GHz) 

Dielectric 

Constant: 

Loss 

Tangent 

Teflon 

la 

3.000 

2.10 

0.0001 

2 

2.728 

2.21 

0.0006 

la 

10.000- 

2.08 

0.0006 

2 

8.250 

2.33 

0.0110 

Polyolefin 

lb 

0.001 

2.36 

0.0002 

lb 

1.300 

2.36 

- 

2 

2.715 

2.38 

0.0011 

2 

8.227 

2.75 

0.0009 

*  la:  Published  tables:  Dielectric  Materials  and  Applications, 

A.  R.  Von  Hippel,  ed.,  MIT  Press,  1954. 
lb:  Manufacturer's  literature:  "Polyguide",  Electronized 

Chemical  Corp. 

2:  Cavity  Measurements 
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le  Holder  Made  From  a  Shortened  Section  of  Waveguide 
Coaxial  Line) 


and  the  (relative)  dielectric  constant  is  defined  as  eA  =  emA0- 
The  material  under  test  (MUT)  has  dielectric  constant 
e  =  e  /e  ,  and  the  conductivity  o  . 

IT  S  O  ® 

The  dielectric  parameters  of  the  sample  can  be  expressed 
in  terms  of  propagation  constants  and  impedances .  With  an 
assumed  time  dependence  exp(jwt),  the  transmission  line 
propagation  constant  is 


Y 


g 


*CV/2 


(5-11) 


where  k  is  the  wavenumber  in  the  dielectric, 
m 


m 


0,(^0^) 


1/2 


(5-12) 


and  k  is  the  cutoff  wavenumber  for  the  dominant  mode  in  the 
c 

line.  For  coax,  k  =0,  and  for  rectangular  waveguide  of  width 

c 

a,  k  =  ir/a.  If  losses  are  significant,  the  propagation  con- 
c  t  , 

stant  may  be  measured  and  expressed  in  the  form  =  a  + 

where  oi  is  an  attenuation  constant.  For  both  TE  and  TEM  modes, 

the  characteristic  impedance  of  a  transmission  line  can  be 

defined  as 


(5-13) 


which  for  TEM  modes  (k  =0)  reduces  to  the  intrinsic  wave 

c 

impedance  of  the  medium. 


Z 


m 


(vv 


1/2 


(5-14) 
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In  the  sample-filled  section  of  the  transmission  line, 
the  propagation  constant  is  given  by 

Ys  -  j<*s2  -  kcV/2  (5-15) 

with 

ks  *  w2uoes  "  ju)yo°s 

=  (‘)^tJloes  ^  ”  jtanfi)  (5-16) 


where  the  loss  tangent  of  the  sample  is 


tan  6 


(5-17) 


For  most  microwave  dielectrics,  the  conductivity  is  much  more 
nearly  constant  with  frequency  than  is  the  loss  tangent.  The 
useful  aspect  of  the  loss  tangent  parameter  is  that,  for  low- 
loss  dielectrics,  the  power  of  a  TEM  wave  is  attenuated  by  a 
factor  exp(-2-ii  tan<$)  in  a  distance  of  one  wavelength.  The 
characteristic  impedance  of  the  sample-f illed  section  of  line 
is 


(5-18) 


The  input  impedance  Z^n  at  the  sample  interface  (the 
plane  A-A')  in  Figure  5.8  is  given  by  the  familiar  formula  for 
a  shorted  line, 

Z.  =  Z  tanh  y  t  (5-19) 

m  s  s 
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By  normalizing  impedances  with  respect  to  Z^,  and  utilizing 
Equations  (5-13),  (5-18),  and  (5-19),  the  following  transcen¬ 
dental  equation  for  the  unknown  propagation  constant  ys  is 
obtained : 


Z.  y  t 

T»  -  V 

g  '  s 


(5-20) 


The  quantities  t,  y  »  and  Z.  /Z  can  be  determined  from  measured 
’  g  in  g 

data. 

The  impedance  ratio  is  related  to  the  reflection  co¬ 
efficient  T  at  the  interface  by 


2  in  _  1  -  T 


(5-21) 


The  impedance  ratio  or  the  reflection  coefficient  may  be 
determined  with  a  slotted  line,  a  reflectometer ,  or  a  micro- 
wave  network  analyzer. 

When  transmission  line  losses  are  negligible,  y may 
be  calculated  from  a  measured  frequency,  using  Equation  (5-11) . 
Otherwise,  the  input  impedance  ratio  Z^n/Zg  for  the  empty  sample 
holder  may  be  measured;  then,  the  propagation  constant  is  given 
by 


Yg  =  t  arc  tanh  tZin/zg) 


(5-22) 


The  transcendental  equation  (5-21)  can  be  solved  for 

Y  by  numerical  iteration  using  Newton's  method.  Let 
s 


"  =  V 


(5-23a) 


2  =  Zin/tZgV) 


(5-23b) 


f(to)  =  to  tanh  to  -  z 


(5-24) 


Then  the  roots  of  f(to)  =  0  are  required.  By  Newton's  mer.hod. 


successive  iterations  of  <on,  are  generated  according  to 


0)  -  =  co  — 

n+1  n 


(5-25) 


and  it  may  be  shown  that 


f(to)  _  smh  uj  -  zu  cosh  u  _ _ v. 

-s-.  ■  =  - r - - — r -  0)  cosh  CU 

f '  (u)  a)  -  cosh  u  sinh  to 


(5-26) 


Equation  (5-25)  has  been  found  to  converge  rapidly  and  accurately. 
Multiple  roots  may  be  resolved  by  inspection  or  by  taking  measure¬ 
ments  with  more  than  one  sample  thickness. 


The  dielectric  constant  and  loss  tangent  are  determined 


from  the  solution  y  ,  using  Equations  (5-11)  and  (5-12)  as 

s 


follows 


.2  2 
k  -  y 

K  =  (1  -  j  tan  6)  =  — - 5— §— 

S  ko 


(5-27) 


where  k  =  /y  e  is  the  free  space  wavenumber. 
000 


i 


Table  5-8  presents  results  of  dielectric  measurements 
for  a  number  of  low-loss  materials.  Impedances  were  measured 
by  slotted-line  techniques.  For  low-loss  materials,  the  most 
significant  error  is  in  the  measurement  of  voltage  standing 
wave  ratio  (VSWR)  or,  equivalently,  in  the  measurement  of  the 
magnitude  of  the  reflection  coefficient.  The  VSWR  is  very  large 
(typically,  greater  than  50dB) ,  and  the  magnitude  of  the  reflec¬ 
tion  coefficient  is  very  nearly  unity. 

A  standing  wave  pattern  is  caused  by  interference 

between  incident  and  reflected  waves,  and  the  null  depths  are 

very  sensitive  to  slight  differences  in  the  wave  amplitudes. 

However,  for  nearly  equal  wave  amplitudes,  the  nulls  are  deep, 

and  they  may  be  hidden  in  noise.  With  laboratory  sources, 

it  is  not  practical  to  increase  the  power  level  enough  to  raise 

the  nulls  above  the  noise  level.  To  overcome  this  problem,  a 

procedure  has  been  developed  to  fit  a  standing-wave  interference 

2 

curve  (i.e.,  a  curve  of  the  form  A+B  cos  kz)  to  a  large  number 
of  measured  data  points.  A  least-squares  curve-fitting  routine 
is  used.  Typically,  the  positions  and  relative  power  levels  of 
four  points  on  each  side  of  a  null  are  recorded.  Data  is 
recorded  at  two  successive  nulls,  to  permit  line  losses  to  be 
measured  and  taken  into  account.  Using  this  procedure,  it  has 
been  found  that  severly  mismatched  impedances  can  be  measured 
more  accurately  with  a  slotted  line  than  with  a  Network  Analyzer. 

5.2.3  Pre-Test  Analysis 

As  the  resonant  cavity  measurement  procedure  has  not 
been  used  extensively  in  the  same  form  as  that  employed  in  this 
program,  analysis  was  undertaken  to  ensure  the  accuracy  and 
applicability  of  the  technique.  This  part  of  Section  5.? 
briefly  highlights  the  analysis  of  the  TMQl0  cavity  case.  A 
more  complete  treatment  is  provided  in  Appendix  C. 
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Table  5.8,  Nominal  Dielectric  Parameters  of  Known  Materials 


Material 

Source 

* 

Frequency 

(GHz) 

Dielectric 

Constant 

Loss 

Tangent 

Nylon  66 

La 

3.0 

3.03 

0.0128 

2 

9.0i 

2.94 

0.0175 

Rexolite  2200 

lb 

10.0 

2.62 

0.0014 

2 

9.0; 

2.68 

0.0012 

Plexiglass 

la 

10.0 

2.59 

0.0067 

2 

9.0 

2.68 

0.0120 

Polystyrene 

la 

10.0 

2.54 

0.0004 

2 

3.0 

2.52 

0.0008 

2 

9.0 

2.54 

0.0003 

Teflon 

la 

10.0 

2.08 

0.0006 

2 

3.0 

2.06 

0.0010 

2 

9.0 

2.04 

0.0004 

*  la:  Published  tables:  Dielectric- Materials  and  Applications, 

A.R.  Von  Hippel,  ed.,  MIT  Press,  1954. 
lb:  Manufacturer's  literature 
2  :  Slotted  line  measurements 
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The  geometry  of  this  problem  is  illustrated  by  Figure 
5.4.  A  right  circular  cylinder  of  radius  "a"  and  height  "d" 
is  partially  filled  from  the  bottom  to  a  height  "b"  with  a 
material  having  a  permittivity  e=ere0  an<*  a  conductivity  o. 

The  dielectric- filled  portion  of  the  cylinder  is  called  region 
1,  while  the  remaining  free  space  portion  is  called  region  2. 
The  coordinate  system  is  oriented  so  that  the  z  -  axis  is  coin¬ 
cident  with  the  longitudinal  axis  of  the  cylinder . 


Figure  5,4,  Circular  Resonator  With  Dielectric  Sample  Filling 
Region  1  (0  <  z  <  b) 


m 


For  this  geometry,  the  dielectric  constant  determination 
can  be  shown  to  be  a  function  only  of  the  shift  in  the  resonant 
frequency  of  the  cavity,  viz: 


(5-28) 


where 

f  =  resonant  frequency  of  the  empty  cavity 
f  =  resonant  frequency  of  the  materially-loaded  cavity 

Figure  5. 5  provides  a  graphical  representation  of  Equation 
(5-28) .  The  overall  sensitivity  of  the  procedure  can  be  esti¬ 
mated  from  the  slope  of  the  curves.  The  greater  the  slope, 
the  greater  the  (f  -  fQ)/f0  value  for  an  incremental  change  in 

the  dielectric  constant.  Since  (f  -  f  )/f  is  an  experimentally 

o  o 

measurable  quantity,  the  ability  to  measure  the  dielectric 
constant  value  is  tied  directly  to  the  precision  with  which  the 
resonant  frequency  shift  can  be  determined.  The  figure  should 
be  interpreted  in  this  light. 

For  example,  one  of  the  cavities  used  in  the  experimen¬ 
tal  program  had  a  height  of  0.20  inches  and  a  radius  of  1.50 
inches.  The  thin,  S-band,  ceramic  samples  have  a  nominal  thick¬ 
ness  of  0.050  inches.  Hence,  the  fill  factor  is  0.25  for  this 
example,  and  from  Figure  5-5  it  can  be  seen  that  the  sensitivity 
of  the  measurement  procedure  degrades  markedly  as  the  nominal 
value  of  the  dielectric  constant  increases. 


FILL  FACTOR 


DIELECTRIC  CONSTANT 

Figure  5.5.  Normalized  Frequency  Shift  vs.  Dielectric  Constant 


The  error  sensitivity  of  the  dielectric  constant  measure¬ 
ment  can  be  found  by  differentiating  Equation  (5-28)  to  yield 


Aer  =  2er 


2  [1  -  b/d  (1  -  l/er)] 


A(fo  -  f) 


(5-29) 


where  Aer  is  the  error  in  the  dielectric  constant  due  to  an 
error  in  the  resonant  frequency  shift  of  A(fQ  ~  f )  .  The  nominal 
(true)  dielectric  constant  of  the  material  is  e^. 

From  the  experimental  results?  a  reasonable  value  for 
A(f  -  f)  at  the  95  percent  confidence  level  is  3.0  MHz  for  an 
empty-cavity  resonant  frequency  of  3014  MHz.  If  these  values  are 
inserted  into  Equation  (5-29)  and  plotted  graphically.  Figure  5.6 
results.  The  fill  factor  was  set  at  the  representative  value 
of  0.25  for  the  construction  of  the  figure. 

The  horizontal  axis  is  the  nominal  value  of  the 
dielectric  constant  of  the  material.  The  bounds  on  the  dielec¬ 
tric  constant  value  are  given  along  the  vertical  axis.  The 
nominal  value  trace  is  found  along  the  center  of  the  three 
lines.  The  upper  bound  at  the  95  percent  confidence  level  on 
the  value  of  the  dielectric  constant  is  given  by  the  top  line. 

The  lowest  of  the  three  lines  details  the  lower  bound  on  the 
dielectric  constant  value  at  the  95  percent  confidence  level. 

For  example,  a  material  with  a  true  (nominal)  dielectric  con¬ 
stant  of  12  will  have  its  value  experimentally  determined 
to  be  lower  than  13  and  greater  than  11  with  a  probability 
of  95  percent. 

The  loss  tangent  determination  is  more  complicated  than 
the  determination  of  the  dielectric  constant.  The  important 
expression  for  the  Q  of  the  materially-loaded  cavity  QL  is 


JL  =  _L  +  _L 

°L  QD  QW 


(5-30) 
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TM  CAVITY  WITH  B/D  =  0,25  A(F  -  f.)  -  3.0  MHZ  F-  =  30W  MHZ 


where 


tan  6 


1  +  e. 


1  -  b/d\ 

573“  ) 


(5-31) 


1.202  n 


(1  +  a/b) 


(5-32) 


In  the  previous  expression,  nQ  is  the  intrinsic  impedance  of 
free  space  and  Rg  is  the  resistivity  of  the  metal  used  to  con¬ 
struct  the  cavity.  Note  that  QL  is  an  observable  in  the 
measurement. 

Equation  (5-30)  is  used  to  plot  as  a  function  of  the 
logarithm  of  the  loss  tangent,  with  the  dielectric  constant 
value  varied  to  produce  the  family  of  curves.  The  conductivity 
of  the  metal  forming  the  walls  of  the  cavity  is  assumed  to  be 
equal  to  one-half  that  of  brass  (a^rass  =  x  107  mhos/meter). 

This  value  is  consistent  with  the  experimental  Q  found  for  the 
empty  cavity.  The  curves  are  given  as  Figure  5.7 . 

As  the  loss  tangent  expression  is  a  function  of  both 
the  nominal  dielectric  constant  value  as  well  as  the  Q  of  the 
cavity,  the  error  analysis  for  the  loss  tangent  is  more  complex 
than  that  for  the  dielectric  constant.  The  errors  in  the  loss 
tangent  have  been  investigated  as  a  function  of  one  of  these 
parameters  while  the  other  was  assumed  to  have  no  associated 
error.  This  analysis  has  shown  that  for  errors  on  the  order 
of  those  apparent  in  the  experimental  work,  errors  in  correctly 
determining  the  loaded  Q  of  the  cavity  dominate  the  error 
expression  for  the  loss  tangent.  The  expression  which  relates 
the  errors  in  the  two  parameters  is 

Mtan  „  =  [M!  ♦  2-^  +  ^  «L  <5-33> 
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1000 


Figure  5.7.  Loaded  Q  Versus  Loss  Tangent,,  °  =  0.5  Brass 


where 


'  *  1  +  er  Hsr1  (5-34) 

Figure  5.8  graphically  details  the  loss  tangent  error  due  to 
errors  in  the  Q  of  the  cavity.  The  conductivity  of  the  metal  is 
assumed  to  be  0.50  brass.  The  error  in  the  Q  of  the  cavity  is 
assumed  to  be  50.  This  is  consistent  with  nearly  a  6  percent 
error  in  the  normalized  bandwidth  of  the  cavity  resonance, 
and  is  on  the  order  of  the  apparent  error  in  the  experimental 
work.  The  fill  factor  is  again  specified  as  0.250. 

The  nominal  loss  tangent  values  are  given  along  the 
horizontal  axis,  while  the  error  bounds  can  be  found  along  the 
vertical  axis.  The  nominal  values  of  the  dielectric  constant 
(er  =  8,  16,  24)  are  shown.  The  upper  three  curves  (in  descending 
order)  are  the  upper  error  bounds  for  er  =  24,  er  =  16,  and 
tr  =  8  respectively.  The  lower  three  curves  (in  descending 
order)  are  the  lower  bounds  for  loss  tangent  errors  for  nominal 
values  of  the  dielectric  constant  equal  to  8,  16  and  24, 
respectively.  Note  that  this  analysis  has  assumed  that  there 
are  no  errors  made  in  determining  the  dielectric  constant  of 
the  material. 

The  minimum  detectable  loss  can  be  viewed  as  the  value 

at  which  the  upper  bound  curves  intersect  the  vertical  axis. 

With  this  definition,  it  appears  that  for  e  =  8  the  minimum 

_  2  r 

detectable  loss  tangent  is  about  2  x  10  ,  while  for  e  =  24 

-3 

the  value  may  be  on  the  order  of  6  x  10 

Although  these  values  indicate  that  the  cavity  measure¬ 
ment  procedure  is  not  as  sensitive  as  some  other  measurement 
procedures,  the  cavity  approach  is  usable  for  very  thin  samples 
while  the  other  methods  generally  are  not.  In  addition,  since 


-103- 


LOG  NOMINAL  LOSS  TANGENT 

Figure  5.3.  Loss  Tangent  Error  Bounds  Caused  by  Uncertainty  in 
Q  Measurement 


the  samples  are  so  thin,  the  calculated  losses  attributed  to 
the  material  are  as  insignificant  at  tan  6  =  10  3  as  they  are 
for  tan  <5  <  10~4.  Of  course,  this  is  only  true  provided  the 
frequency  is  low  enough  so  that  the  material  thickness  is  much, 
much  smaller  than  a  wavelength.  For  the  S-Band  and  X-Band 
frequencies  of  interest  here,  the  material  thickness  constraint 
is  always  satisfied. 


5.3 


Array-Hardening  Materials  Interaction  Analysis 

The  materials  under  consideration  in  this  study  would 
primarily  be  used  in  either  of  two  ways:  (1)  dielectric  wave¬ 
guide  plugs  or  (2)  dielectric  slabs  covering  an  array.  It  is 
well  known  in  the  array  design  field  that  the  application  of  any 
material  in  either  of  these  ways  has  important  implications 
regarding  the  performance  of  an  array.  In  this  section,  the 
interactions  between  the  hardening  material  and  array  performance 
parameters  are  considered.  The  issues  of  reflection  losses 
and  dissipative  losses  is  addressed  briefly,  while  somewhat 
greater  emphasis  is  placed  on  the  issue  of  array  blindness.  The 
theoretical  analysis  of  the  array  blindness  question  is  pre¬ 
sented  in  Appendix  D  and  only  the  implications  of  that  analysis 
are  presented. 

5.3.1  Mechanisms  of  Array  Performance  Degradations 

As  a  first  order  approximation,  there  are  three  phenomena 
which  contribute  to  the  degradation  of  array  performance  when 
that  array  has  either  a  covering  dielectric  slab  or  dielectric 
waveguide  plugs.  These  are  (1)  reflection  losses,  (2)  dissipa¬ 
tive  losses  and  (3)  array  blindness.  A  schematic  representation 
of  the  degradation  mechanisms  is  given  by  Figure  5«-9  for  a 
dielectric  covered  array.  The  phenomena  are  essentially  the 
same  for  an  array  with  waveguide  plugs. 

The  figure  portrays  an  array  of  apertures  containing 
incident  energy.  A  portion  of  that  energy  is  reflected  at  the 
slab  interface,  while  the  remaining  energy  passes  into  the  slab. 
In  the  slab,  some  of  the  incident  energy  iL  onsumed  in  material 
losses  while  the  rest  of  the  energy  is,  in  general  divided 
between  the  desired  main  antenna  beam  and  perhaps  a  "leaky  wave" 
in  the  slab  structure.  If  the  leaky  wave  is  launched,  the  energy 
in  that  wave  propagates  out  of  the  slab  along  its  entire  extent 
and  interferes  in  a  vector  sense  with  the  main  beam.  At  those 


INCIDENT 

ENERGY 


Figure  5.9,  Schematic  of  Array  Performance  Degredations 
Mechanisms  for  a  Dielectrically  Covered  Array 


■ 
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array  scan  angles  at  which  the  leaky  wave  exactly  cancels  the 
main  beam,  the  array  is  said  to  have  a  blind  spot;  the  array 
will  neither  transmit  nor  receive  energy  in  that  direction. 

5.3.2  Reflection  and  Dissipative  Losses 

A  simple  model  of  the  reflection  and  dissipative  losses 
in  a  radome  is  the  classical  problem  of  a  plane  wave  propagating 
through  a  three-layered  structure .  The  analysis  is  not 
repeated  here,  but  the  results  were  used  to  generate  the  curves 
of  Figures  5.10  and  5.11.  The  two  figures  portray  the  reflec¬ 
tion  and  dissipative  losses  associated  with  slabs  whose  dielec¬ 
tric  constants  are  2.8  and  7.0  respectively.  Curves  for  various 
values  of  the  loss  tangent  are  included  on  each  figure.  Note 
that  the  dielectric  constant  value  of  2.8  is  approximately  that 
of  the  silicone  rubber  class  of  materials,  while  the  value  of 
7.0  is  consistent  with  that  for  the  plasma-sprayed  Alumina 
(A^O^)  and  Spinel  (Al^O^  and  MgO)  . 

The  two  figures  are  important  because  they  depict  the 

attenuation  as  a  function  of  the  material  thickness,  and  for 

these  hardening  applications  the  thicknesses  will  change  as  the 

material  is  boiled  off  or  otherwise  removed.  Thus,  for  the 

-2 

silicone  rubber  materials  whose  loss  tangent  is  about  10  ,.  quite 

substantial  thickness  variations  may  occur  with  the  associated 

losses  ranging  from  0  dB  to  1  dB.  For  the  ceramic  materials  with 
—  2 

a  10  loss  tangent,  on  the  other  hand,  the  losses  will  range 
from  O  dB  to  nearly  4  dB.  Note  that  a  3  dB  figure  corresponds  to 
a  loss  of  one-half  of  the  available  power  upon  transmission, 
and  again  upon  reception.  Such  losses  are  significant  and 
perhaps  unacceptable,  even  in  hardened  radar  applications  in 
which  significant  losses  may  be  accepted  as  a  cost  of  hardening. 

5.3.3  Array  Blindness 

Of  the  issues  which  are  being  considered  in  this  section, 
the  question  of  array  blindness  is  the  most  important.  In  fact. 
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THICKNESS  [freespace  wavelengths] 

Figure  5,10.  Radome  Efficiency  for  a  Flat  Array  Covered 
by  a  Slab  With  a  Dielectric  Constant  of  2.8 


THICKNESS  [fVeespace  wavelengths] 

Figure  5.11.  Radome  Efficiency  for  a  Flat  Array  Covered  by  a 
Slab  With  a  Dielectric  Constant  of  7.0 


for  detection,  tracking  and  fire  control  applications,  the 
absence  of  blindspots  is  essential.  A  great  deal  of  analytical 
effort  has  been  expended  on  this  question (15'16' 17 )  .  Many  of 
the  details  of  complex  antenna  arrays  do  not  lend  themselves, 
however,  to  exact  analysis. 

A  simplified  model  of  the  problem  is  given  by  Figure  5.12, 
in  which  infinitely  long  parallel  plate  waveguides  are  covered 
by  a  slab  of  thickness  T.  The  radiating  structure  is  an  array 
of  thin  slots  with  width  B  and  spacing  D.  A  linear  phase  shift 
along  the  waveguide  structure  allows  the  main  beam  to  be  formed 
and  scanned. 


Figure  5.12.  Covered  Array  Model 

For  narrow  slot  widths,  the  slab  may  be  treated  as  if 
it  were  lying  on  a  ground  plane,  and  it  is  then  possible  to 
show  (see  Appendix  D)  that 


where 


t  =  (sin  6  +  g^)  2  (5-36) 

That  is,  for  a  given  required  scanning  range  +0,  the  maximum 
thickness  in  free  space  wavelengths  of  the  slab  which  avoids 
blind  spots  within  the  internal  (-0,  0)  is  given  by  equation 
(5-35) . 

Note  that  this  is  the  maximum  tolerable  slab  thickness 
which  avoids  blindness  within  the  desired  scan  angles.  Resource¬ 
ful  antenna  designers  can  usually  trade  off  some  other  perfor¬ 
mance  to  increase  the  maximum  thickness.  Since  the  hardening 
material  is  expected  to  recess  during  nuclear  blasts,  the 
effects  of  the  changing  slab  thickness  are  also  important  and 
have  not  been  explicitly  considered  in  this  analysis. 

It  is  informative  to  plot  equation  (5-35)  as  a  function 
of  the  dielectric  constant.  Figure  5.13  provides  such  a  plot 
with  several  maximum  scan  angles  shown.  Consider  the  0  =  45 
degrees  curve,  for  example:  a  maximum  thickness  of  0.165  X  can 
be  used  when  the  slab  has  a  dielectric  constant  of  2.9  (silicone 
rubber).  For  a  7.0  dielectric  constant  value  (Alumina  or  Spinel), 
the  corresponding  maximum  slab  thickness  is  0.082  X. 

Table  5*9  provides  the  maximum  slab  thicknesses  for 

these  two  dielectric  constant  values  when  the  maximum  array 

scan  angle  is  either  +45  degrees  or  +60  degrees.  The  maximum 

allowable  thickness  is  given  both  in  wavelengths  and  in  the 

corresponding  number  of  inches  for  frequencies  of  3.0  GHz  and 

9.0  GHz.  From  the  Table,  it  can  be  seen  that  for  the  silicone 

rubber  class  of  materials  at  X-Band,  the  maximum  thickness  of 

the  slab  is  0.22  inches  for  0  =  +45  degrees  and  0.16  inches 

max  — 

for  0  *  +60  degrees.  For  the  applications  of  greatest  in- 

max  * 

terest,  scan  angles  of  at  least  +45  degrees  are  desired.  For 


1 1 ELECTRIC  CONSTANT 

Slab  Thickness  vs.  Dielectric  Constan 
dspot  Angles  Greater  Than  Specified  S< 


the  ceramics  at  X-Band,  the  maximum  thickness  is  0.12  inches  for 
0  =  +45  degrees  and  0.09  inches  for  9  =  +60  degrees.  Of 

the  ceramic  materials  which  were  subjected  to  the  electrical 
tests,  no  samples  were  found  with  a  thickness  greater  than  0.09 
inches.  Thus,  for  either  type  of  material,  there  is  good  reason 
to  believe  that  these  materials  may  be  useful  for  array  hardening 
from  an  electrical  point  of  view,  provided  that  the  thicknesses 
which  have  been  discussed  are  adequate  for  the  thermal  protec¬ 
tion  of  the  array.  If  greater  thicknesses  are  required,  then 
the  array  must  be  "tuned"  to  tolerate  the  greater  thickness. 
However,  recession  of  the  layer  due  to  thermal  exposure  would 
tend  to  affect  the  performance  of  such  an  array. 

5.4  Conclusion 

The  major  accomplishments  and  conclusions  of  the 
electrical  properties  test  and  analysis  task  are  as  follows: 

1.  A  cavity  measurement  procedure  for  plasma-sprayed 
ceramic  materials  has  been  developed  and  refined. 

2.  Electrical  properties  have  been  measured  for  plasma- 
sprayed  ceramics  at  S-Band  and  at  X-Band. 

3.  Electrical  properties  have  been  measured  for  silicone 
rubber  materials  (including  their  chars)  at  S-Band 
and  at  X-Band. 

4.  For  the  material  thicknesses  considered  in  this 
study,  the  silicone  rubber  materials  appear  promising 
from  an  electrical  point  of  view.  Of  the  plasma- 
sprayed  ceramics,  the  Alumina  (A^OO  and  Spinel 
appear  to  be  the  most  promising,  again  from  an 
electrical  perspective. 

5.  The  hardening  material/array  interaction  question 
must  be  considered  as  an  integral  part  of  the  array 
design  process. 

While  the  results  to  date  show  promise,  a  significant 
amount  of  work  remains  to  be  undertaken  before  the  usefulness 
of  these  materials  may  be  adequately  assessed.  In  particular: 

1.  Larger  numbers  of  samples  of  each  promising 
material  must  be  evaluated. 
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2 .  The  char  structure  of  the  silicone  rubber  materials 
must  be  quantitatively  investigated. 

3.  The  array  interaction  questions  must  be  treated 

in  greater  depth.  This  can  be  accomplished  only  by 
considering  specific  array  configurations  in  detail 

4.  The  electrical  variations  in  the  materials  must  be 
found  for  such  large  structures  as  arrays. 

5.  The  thickness  variations  in  the  materials  must  be 
found  for  such  large  structures  as  arrays. 

6.  The  implications  of  (4)  and  (5)  on  the  array  per¬ 
formance  must  be  determined  in  an  analytical  manner 

7.  A  proof-of-principle  model  would  be  instructive  in 
assessing  the  viability  of  specific  hardening  con¬ 
cepts  . 
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6.0 


MECHANICAL  PROPERTIES  TEST 


6 . 1  Mechanical  Properties  Requirements 

In  the  structural  analysis  of  structures  with  thermal 
protection,  it  is  necessary  to  calculate  and  predict  how  the 
structure  will  survive  the  thermal  pulse  and  the  shock  wave 
pulse  caused  by  a  nuclear  burst.  The  primary  type  of  construc¬ 
tion  under  investigation  consists  of  a  metal  surface  coated 
with  a  plasma  sprayed  ceramic.  Typically,  the  ceramic  is  thin, 
in  the  order  of  0.125  inches  or  less  and  the  metal  substrate  is 
at  least  .125  inches  thick.  It  may  even  be  on  the  order  of  one 
inch  thick  on  parts  of  the  radar  assembly. 

During  the  thermal  pulse,  the  ceramic  heats  up  rapidly, 
especially  near  its  outer  surface,  but  the  metal  substrate 
heats  relatively  slowly  below  the  ceramic.  During  heating, 
the  ceramic  and  metal  expand  or  try  to  expand  according  to  the 
relationship 

6  =  a£At 

6  =  total  expansion,  inches 

a  =  expansion  coefficient,  inches/inch  °F 
At  =  temperature  rise,  °F 

Because  the  outer  layer  of  ceramic  is  very  hot,  on  the  order 
of  its  melting  point,  t  5  5000°F  for  Zr02*  and  the  inner  layers 
next  to  the  steel  are  relatively  cool,  initially  a  differen¬ 
tial  expansion  exists  across  the  ceramic  layer.  This  is  res¬ 
trained  by  shear  forces  in  the  ceramic  and  metal-ceramic  bond. 
Thus  compressive  stresses  are  generated  in  the  outer  layer  and 
tensile  stresses  in  the  inner  layer  as  well  as  shear  stresses 
in  the  lamanae. 

The  elongation  due  to  such  stress  at  constant  uniform 
temperature  is  given  by  the  relationshiD  6  ■ 

E 


5  =  elongation,  inches 

s  =  unit  stress  (compressive  or  tensile),  psi 
E  =  modulus  of  elasticity,  (psi)  at  specific  temperature 

The  relationship  shown  is  for  a  simple  one-dimensional  case. 
Poison's  ratio  is  used  to  calculate  the  effect  of  two  and 
three-dimensional  stress  and  strain. 

The  values  of  o,  E,  S^,  Sc,  Sg  and  V  are  used  to  cal¬ 
culate  thermal  stress  and  failure  criteria  in  the  structures. 
Therefore,  as  a  minimum,  it  is  necessary  to  have  the  follow¬ 
ing  data  for  the  ceramics: 

1.  Modulus  of  Elasticity 

2.  Tensile  Strength 

3.  Compressive  Strength 

4.  Shear  Strength  (Parallel  to  Substrate) 

5.  Bond  Strength  (Shear  and  Tension) 

6.  Poisson's  Ratio 

7.  Linear  Coefficient  of  Expansion  (Three  Axes) 


It  is  desirable  to  obtain  the  following  properties 
for  the  plasma-sprayed  ceramic  as  a  function  of  at  least  three 
temperatures:  Room  temperature,  1500 °F  and  3000 °F. 


Desired  Data 

1.  Modulus  of  elasticity 

a.  Tensile  (parallel  to  layer) 

b.  Compressive  (parallel  to  layer) 

c.  Tensile  (normal  to  layer) 

2.  Ultimate  strength 

a.  Tensile  (parallel  to  layer) 

b.  Compressive  (parallel  to  layer) 

c.  Tensile  (normal  to  layer) 
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3.  Bond  shear  strength  (parallel  to  layers) 

4.  Adhesive  strength  (normal  to  layers) 

5.  Linear  thermal  expansion  coefficient 

6.  Poisson's  Ratio 

7.  Open  and  closed  porosity  sizes  and  ratios 

8.  Density 

9 .  Shear  strength 

Most  of  these  values  must  be  determined  as  a  function  of  tem¬ 
perature  from  room  temperature  to  as  close  to  the  melting  point 
as  is  practical.  Becuase  the  above  mentioned  properties  are 
functions  of  the  processing,  it  is  important  to  determine  their 
values  for  specimens  manufactured  by  the  same  process  that 
would  be  used  to  coat  the  radar  surfaces. 


6 . 2  Mechanical  Properties  Test  Program  Plan 

Mechanical  properties  testing  will  be  conducted  by  the 
BMD  Materials  Program  test  facility  at  AMMRC.  SAI  has  fabricated 
the  bulk  specimens,  exposed  them  to  a  predetermined  thermal  en¬ 
vironment,  and  suggested  required  properties  and  temperatures 
to  be  obtained.  A  standard  plasma  spray  structural  test  specimen 
was  designed  and  fabricated  by  SAI.  It  was  a  thin  metal  plate 
four  inches  square  coated  with  plasma  sprayed  ceramic  approximately 
0.125  inches  thick.  Figure  4.11  illustrates  a  typical  specimen. 


This  specimen  was  designed  to  be  cut  into  smaller  shapes  for 
actual  test  by  means  of  a  diamond  coated  abrasive  wheel.  The 
actual  shape  and  size  of  each  sub  specimen  to  be  determined 
by  AMMRC.  Forty  of  these  specimens  were  fabricated  and  were 
to  be  exposed  to  thermal  environments  as  described  in  Table 
6.1 

TABLE  6.1 

Structural  Test  Coupon  Matrix 
All  Specimens  4"  Square 


Spec. 

Test 

No.  of 

Coating 

No. 

Type 

Exposures 

Material 

Infiltrant 

52 

No  Exposure 

0 

AI2O3 

Silicone 

53 

Radiation 

1 

n 

II 

54 

•I 

2 

n 

II 

55 

Fireball 

1 

11 

II 

56 

It 

2 

11 

II 

58 

No  Exposure 

0 

•1 

0 

59 

Radiation 

1 

n 

0 

60 

II 

2 

11 

0 

61 

Fireball 

1 

ii 

0 

62 

II 

2 

11 

0 

64 

No  Exposure 

0 

Spinel 

Silicone 

65 

Radiation 

1 

II 

II 

66 

II 

2 

II 

II 

67 

Fireball 

1 

II 

II 

68 

II 

2 

II 

II 

70 

No  Exposure 

0 

II 

0 

71 

Radiation 

1 

It 

0 

72 

11 

2 

11 

0 

73 

Fireball 

1 

11 

0 
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Table  6.1  (Cont'd) 


Spec. 

NO. 

Test 

Type 

No.  of 
Exposures 

Coating 

Material 

Inf iltrant 

74 

Fireball 

2 

Spinel 

0 

76 

No  Exposure 

0 

Zr02 

Silicone 

77 

Radiation 

1 

It 

II 

78 

II 

2 

II 

II 

79 

80 

Fireball 

2 

II 

II 

•1 

II 

82 

No  Exposure 

0 

n 

0 

83 

Radiation 

1 

ii 

0 

84 

f« 

2 

w 

0 

85 

Fireball 

1 

n 

0 

86 

tl 

2 

ii 

0 

88 

No  Exposure 

0 

Zr02  MgO 

Silicone 

89 

Radiation 

1 

II 

It 

90 

tl 

2 

II 

II 

91 

Fireball 

1 

II 

II 

92 

•t 

2 

II 

II 

94 

No  Exposure 

0 

II 

0 

95 

Radiation 

1 

II 

0 

96 

It 

2 

ti 

0 

97 

Fireball 

1 

n 

0 

98 

II 

2 

it 

0 

100 

No  Exposure 

0 

Zr02-Y203 

Silicone 

101 

Radiation 

1 

II 

II 

102 

n 

2 

II 

II 

103 

Fireball 

1 

II 

n 

104 

If 

2 

II 

n 

106 

No  Exposure 

0 

II 

0 

107 

Radiation 

1 

It 

0 

108 

It 

2 

II 

0 

109 

Fireball 

1 

II 

0 

110 

If 

2 

If 

0 
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Some  of  the  specimens  were  to  be  tested  before  heating  in  the 
TRS  facility  to  obtain  results  on  virgin  material.  Similar 
structural  tests  were  to  be  performed  on  once  heated  and  twice 
heated  ceramic  to  see  what  if  any  property  damages  were  caused 
by  the  radiation  and  convective  heating.  The  temperature  cycling 
will  also  be  done  at  AMMRC  by  cyclic  heating  the  specimens. 
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7.0  CONCLUSIONS  AND  RECOMMENDATIONS 

The  conclusions  derived  from  this  study  are  presented 
under  the  categories  of  (1)  test  techniques,  (2)  material  per¬ 
formance,  and  (3)  radar  hardening  implications. 

Test  Techniques 

The  nuclear  thermal  radiation  environment  design  cri¬ 
teria  is  so  severe  that  no  currently  available  test  facility 
can  simulate  the  combined  flux  and  fluence  level.  The  Solar 
Furnace  is  a  proven  technique  for  obtaining  the  desired  flu¬ 
ence  level  but  is  limited  in  flux  level  and  size  of  compon¬ 
ents  which  can  be  tested.  The  SAI  Pulse  Type  Thermal  Radia¬ 
tion  Simulator  shows  promise  for  simulating  the  higher  flux 
levels  and  for  testing  larger  components  but  is  limited  in  total 
fluence.  This  facility  is  also  in  the  development  stage  and 
requires  extensive  calibration  before  it  can  be  used  for  rou¬ 
tine  testing  with  high  confidence.  A  continuous  burning  TRS 
with  high  fluence  capability  is  in  development. 

A  cavity  measurement  procedure  was  designed  for  elec¬ 
trical  properties  tests  for  thin  specimens  at  s-band  and  x- 
band  frequencies.  Analyses  confirmed  that  adequate  measure¬ 
ment  sensitivity  was  available  for  obtaining  dielectric  constant 
and  loss  tangent  values  for  the  materials  under  evaluation. 

Material  Performance 

The  infiltrated  (silicone)  porous  ceramic  coatings 
had  poor  thermal  performance.  The  silicone  forms  a  low  car¬ 
bon  char  when  it  decomposes  under  thermal  radiation  heating. 

This  residual  char  raises  the  absorptivity  of  the  surface  to 
high  values  which  result  in  excessive  material  heating  and 
melting.  Uninfiltrated  porous  ceramic  coatings  are  expected 


to  provide  adequate  thermal  performance  if  the  absorptivity 
can  be  maintained  at  0.2  or  less.  This  level  has  been  ob¬ 
tained  for  AI2O3  and  Zr(>2  coatings  in  previous  program  tests. 

The  ceramics  experienced  some  thermal  stress  failures 
during  heating.  This  remains  a  major  concern;  however,  anal¬ 
yses  suggest  that  pre-heating  the  base  metal  during  plasma 
spraying  of  the  coating  may  reduce  tensile  failures  during 
subsequent  thermal  exposure.  The  electrical  property  test  re¬ 
sults  showed  loss  tangent  values  less  than  0.03  but  reasonably 
large  values  of  dielectric  constant  (6-20)  for  the  ceramic 
coating  materials. 

The  silicone  ablator  showed  excellent  thermal  perform¬ 
ance.  The  thick  char  allows  a  high  surface  temperature  and 
thus  high  re-radiation  of  incident  heat  flux  while  the  insu- 
lative  properties  of  the  char  and  virgin  material  protect  the 
underlying  structure.  The  char  was  found  to  contain  a  low 
fraction  7%)  of  carbon  which  degraded  the  electrical  prop¬ 
erties  to  some  extent.  The  virgin  electrical  properties  tests 
showed  low  loss  tangent  and  dielectric  constants. 

Other  materials  were  tested  thermally  with  the  follow¬ 
ing  results.  An  AS-3DX  composite  partially  melted  but  sur¬ 
vived  the  solar  furnace  test  environment.  The  BN  (HBR  Grade) 
experienced  thermal  stress  cracking  in  the  solar  furnace  tests. 
The  teflon  specimens  were  observed  to  absorb  significant  rad¬ 
iation  in-depth  and  free  carbon  forms  on  the  surface  during 
heating.  Insufficient  thermal  data  was  obtained  on  a  number 
of  bulk  ceramic  and  uninfiltrated  porous  ceramic  coating  ma¬ 
terials  due  to  the  lack  of  adequate  heating  during  the  Thermal 
Radiation  Simulator  tests  and  the  limited  number  of  specimens 
which  could  be  tested  in  the  Prench  Solar  Furnace  during  the 
designated  test  period. 
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Radar  Hardening  Implications 


The  uninfiltrated  ceramic  coatings  are  viable  candi¬ 
dates  for  protection  of  radar  components.  Analyses  show  that 
absorptivities  of  0.2  or  less  must  be  maintained  in  order  to 
survive  the  thermal  design  criteria  for  typical  coating  thick¬ 
nesses.  The  electrical  properties  of  the  ceramic  are  not  as 
good  as  desired  but  appear  adequate  if  an  integral  hardened 
radar  were  designed.  A  remaining  problem  which  must  be  ad¬ 
dressed  is  the  potential  thermal  stress  failure  due  to  multi¬ 
ple  cycle  heating.  Pre-heating  during  processing  may  solve 
this  potential  problem. 

The  silicone  rubber  ablative  material  looks  very  prom¬ 
ising  for  radar  hardening.  The  electrical  properties  of  the 
virgin  material  are  very  good.  The  low  carbon-content  char 
shows  some  degraded  electrical  characteristics  but  not  sig¬ 
nificant  enough  to  preclude  use  of  this  material.  Materials 
similar  to  this  have  been  used  for  radome  applications  and 
combustion  chamber  wall  ablative  materials.  In  addition,  the 
material  appears  practical  because  it  is  inexpensive,  easy  to 
apply,  adheres  well,  resists  weathering  and  would  provide  a 
moisture  barrier.  It  also  would  be  less  affected  by  surface 
contaminents  especially  due  to  the  expected  dust  environment 
during  a  nuclear  attack. 

It  is  recommended  that  further  work  be  performed  to 
evaluate  uninfiltrated  ceramic  coatinqs  and  the  silicone  rub¬ 
ber  class  of  materials.  Improved  process  techniques  must  be 
developed  for  the  plasma  sprayed  ceramic  coating  to  obtain 
desired  porosity,  surface  absorptivity  and  to  minimize  thermal 
stress  problems.  Alternative  silicone  rubber  materials  with 
varying  degrees  of  char  strength,  thermal  characteristics  and 
electrical  characteristic  should  be  evaluated.  Preliminary 
thermal  screening  test  are  recommended  during  the  material 


process  development  and  material  selection  phase.  Larger 
scale  thermal  tests  which  would  use  the  French  Solar  Furnace 
or  updated  WSMR  Solar  Furnace  are  recommended  for  evaluation 
of  the  subset  of  materials  which  passes  the  screening  test 
phase.  Electrical  property  tests  at  x-band  frequencies  should 
be  performed  on  these  materials.  In  concert  with  these  tests, 
analyses  of  these  material  concepts  should  be  performed  for  a 
hardened  integral  array  preliminary  design  because  it  does  not 
appear  possible  to  simply  identify  materials  to  harden  existing 
DAR  designs  without  significant  performance  degradation.  By 
incorporating  hardening  materials  into  the  initial  preliminary 
design,  performance  degradation  can  be  minimized. 
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APPENDIX  A  -  Material  Thermal  Properties 


This  appendix  presents  the  thermal  data  used  for  all 
analyses  of  the  materials. 


Al.  Ceramic  Coatings  and  Steel  Substrates 

Table  Al  presents  the  alumina  coating  thermal  data.  For 
porous  alumina,  the  thermal  conductivity  was  found  to  be  a  strong 
function  of  the  material  structure  (Reference  Al) .  For  example, 
according  to  Reference  Al,  plasma  sprayed  alumina  with  15%  porosity 
has  a  thermal  conductivity  only  10  to  25%  that  of  bulk  alumina. 

This  indicates  that  the  conductivity  of  porous  alumina  cannot 
accurately  be  expressed  as : 


(1 


rp)kb 


(Al) 


where  kp  =  conductivity  of  a  porous  layer  of  material 
kfa  =  conductivity  of  non-porous  bulk  material 
np  =  porosity 

Reference  Al  also  showed  that  the  conductivity  of  a  plasma 
sprayed  AX^Oj  layer  increases  by  heating  it  to  near  the  melt  tem¬ 
perature  and  then  cooling  it  to  room  temperature.  After  five  such 
cycles,  the  conductivity  was  just  20%  less  than  that  of  non-porous 
bulk  alumina. 


The  zirconia  thermal  data  are  presented  in  Table  A2.  The 
conductivity  of  porous  zirconia  was  not  known  and  so  equation  Al 
was  used  to  establish  it. 


Table  A3  presents  the  stainless  steel  substrate  and  the 
nichrome  bond  thermal  properties. 
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Table  AI.  Plasma  Sprayed  Alumina  Thermal  Properties 


Density:  230  lb/ft3  (non-porous) 

Heat  of  Fusion:  461  Btu/lb 
Melt  Temperature:  4173  R 


Thermal  Conductivity  and  Specific  Heat 


Temperature 

°R 

Non-Porous 
Conductivity 
(BTU/f t-hr-°R) 

15%  Porous 
Conductivity 
(BTU/f t-hr-°R) 

Specific 

Heat 

(BTU/lb-°R) 

400 

24.1 

2.5 

0.18 

800 

12.6 

2.5 

0.23 

1200 

7.2 

1.8 

0.26 

1600 

4.7 

1.6 

0.29 

2400 

3.1 

1.6 

0.33 

4177 

3.1 

1.6 

0.42 

Data  From  Reference  Al. 
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Table  A2.  Plasma  Sprayed  Zirconia  Thermal  Properties 


Density:  359  lb/ft^  (non-porous) 

Heat  of  Fusion:  303.8  BTU/lb 
Melt  Temperature:  5310°R 


Thermal  Conductivity  and  Specific  Heat 


Temperature 

(°R) 

Non-Porous 

Conductivity 

(BTU/ft-hr-OR) 

Specific 

Heat 

(BTU/lb-°R) 

400 

0.88 

0.1 

800 

0.97 

0.13 

1600 

1.12 

0.15 

2400 

1.30 

0.16 

3200 

1.  48 

0.15 

4000 

1.82 

0.15 

Table  A3,  Stainless  Steel  and  Nichrome  Bond  Thermal  Properties 


303  Stainless 

Nichrome  Bond 

Density  (lb/ft^) : 

494 

525 

Melt  Tempetature  (°R) : 

3010 

3012 

Temperature  (°R) 

Conductivity 

(BTU/ft-hr-°R) 

400 

8.5 

4.7 

1000 

10.9 

6.8 

2000 

13.6 

10.1 

3000 

13.6 

19.4 

Temperature  (°R) 

Specific  Heat 

(BTU/lb-°R) 

400 

.12 

.11 

1000 

.13 

.11 

2000 

.16 

.11 

3000 

.16 

.11 
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A2 .  Silicone  Rubber 

The  thermal  properties  for  the  Sylgard  184  silicone 
resin  used  as  the  infiltrant  and  Silastic  E  silicone  rubber 
used  as  an  ablator  are  presented  in  Table  A4.  The  decomposi 
tion  of  the  silicone  resin  was  modeled,  using  a  heat  of 
decomposition  determined  from  Dow  Corning  data  based  on 
plasma  torch  tests.  The  kinetics  scheme  was  approximated 
using  the  decomposition  rate  data  for  natural  rubber,  three 
reinforced  polymeric  materials,  and  carbon  phenolic,  as 
shown  in  Table  A4.  Based  on  this  comparison,  it  was  decided 
that  the  three-legged  carbon  phenolic  decomposition  scheme 
would  be  an  adequate  representation  of  silicone  resin  infil¬ 
trant.  The  frequency  factors,  A^  and  the  activiation 
energy,  ,  needed  in  equation  4-2  of  Section  4.1.1,  were 
evaluated  from  this  data.  At  was  determined  to  be  5.1xl04 

I  .  1 

sec  ,  however,  this  resulted  in  decomposition  at  room 
temperature.  This  behavior  is  unrealistic  and  a  result  of 
an  empirical  smoothing  of  actual  kinetics  data  at  low  tem¬ 
peratures.  To  avoid  this  problem.  At  was  set  equal  to 
5.1x10  secA,  which  essentially  eliminates  the  low  temper¬ 
ature  reaction. 


Table  A4,  Silicone  Rubber  Thermal  Properties 


Propertv 


Sylgard  184 
Silicone  Resin 


Silastic  E 
Silicone  Rubber 


Density  (lb/ft  )  65.5 

Specific  Heat  @  77°F 
(BTU/lb-F)  0.34 

Thermal  Conductivity 
(BTU/ft-hr-  F)  0.09 

Heat  of  Decomposition 
(BTU/lb)  4500 


69.3 


0.34 


0.12 


ftfMt  I 

OCCV*'PQS272(K  R2MET2CS 


Reaction  Kinetics  ** 

Frequency  Factor,  Aj_  (Sec  x) 
5 . Ixl0~6 ,  85.4,  1.01125xl06 


Activation  Energy,  E^  (BTU/lb-mole)  1 


1208,  25200,  56041 


1*  M  )•  I  •  « 

O/TI  B  ir4  1*111 


Properties  from  Dow  Corning  data  sheets 

Modified  three-legged  carbon  phenolic  decomposition  scheme 


Bulk  Materials  Thermal  Properties 


A  series  of  bulk  materials  intended  to  serve  as  RF  window 
materials,  structural  materials,  or  thermal  protection  materials 
were  investigated.  Their  thermal  properties  are  presented  in 
Table  A5  and  Figure  Al. 


Table  A5.  Bulk  Materials  Thermal  Properties 
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APPENDIX  B 


SUMMARY  OF  CAVITY  RESULTS 
FOR  ALL  TEST  RUNS 


This  appendix  gives  the  results  of  each  test  run  completed 
for  each  sample.  Table  B»1  gives  the  results  of  S— band  tests ( 
and  Table  B.2  gives  the  X-band  results.  Where  four  or  more  runs 
are  indicated,  an  average  (Avg.)  and  standard  deviation  (S.D.) 
are  provided.  In  computing  these  statistics,  the  high  and 
low  values  of  er  and  tanc  have  been  discarded.  Where  only  two 
runs  are  shown,  the  results  are  averaged.  A  total  of  seventy- 
eight  runs  are  tabulated. 

The  thickness  of  each  sample  is  given  in  Table  B.3.  In 
the  table,  the  thickness  value  is  listed  in  the  column.  The 
height  of  the  end-plate  shoulder  is  denoted  as  D2"  Figure  B-l 
provides  the  cross-section  of  a  specimen  end-plate  showing  the 
location  of  these  parameters. 

Two  cavity  heights  were  used  for  the  S-band  tests,  while 
only  one  was  used  at  X-band.  The  nominal  heights  at  S-band 
were  0.200  inch  and  1.00  inches.  The  run  numbers  that  correspond 
to  the  larger  S-band  cavity  are  preceded  by  an  asterisk C*) 


Table  B.l.  Tabulation  of  All  S-Band  Cavity  Measurements 
for  SAI  Materials 


Table  B.l.  Tabulation  of  All  S-Band  Cavity  Measurements 
for  SAI  Materials  (Continued) 
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Table  B  .3 .  Ceramic  Sample  Thicknesses 


SAI 

Frequency 

Sample 

Sample  If 

Band 

D1 

°2 

ZrO,  -MO 

2  g 

2 

S 

0.056 

0.053 

ZrO?  -  12ZY20^ 

3 

S 

0.060 

0.053 

A12°3_  2%  Ca02 

4 

s 

0.059 

0.054 

Zr02 

5 

S 

0.058 

0.056 

Spinel 

6 

s 

0.058 

0.057 

Zr02 -  Ca02 

201 

s 

0.082 

0.041 

ZrO-M  0 

2  g 

210 

s 

0.069 

0.051 

Spinel 

356 

s 

0.084 

0.040 

Zr02  -  Ca02 

201 

X 

0.059 

0.054 

Zr°2  -Y203 

202 

X 

0.064 

0.051 

ZrO,  -MO 

2  g 

210 

X 

0.065 

0.052 

A!203-  Ti02 

333 

X 

0.065 

0.052 

Al^O^-TiO^urnec 

i  333 

X 

0.056 

0.052 

Zr02~  Y203burned 

202 

X 

0.052 

0.010 

*  The  parameters  and  are  indicated  in  Figure  B-l. 


Dimenions  are  in  inches. 
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APPENDIX  C 


ANALYSIS  OF  A  CYLINDRICAL  CAVITY  MEASUREMENT 
PROCEDURE  FOR  THE  DETERMINATION  OF  THE  DIELECTRIC 
CONSTANT  AND  THE  LOSS  TANGENT  OF  A 
THIN  DIELECTRIC  SLAB  MOUNTED  ON  A  METAL  DISC 
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1.0 


INTRODUCTION  AND  SUMMARY 


The  process  for  producing  the  plasma-sprayed  ceramic 
materials  required  the  use  of  a  metal  support  or  backing  plate. 
In  addition,  it  was  desired  to  keep  the  ceramic  materials  as 
thin  as  possible  so  as  to  minimize  deleterious  effects  on  the 
array  performance.  At  X-band  (9  GHz),  for  example,  the  ceramic 
thicknesses  were  expected  to  be  less  than  one-eighth  of  a  wave¬ 
length.  Hence,  conventional  bulk  measurement  techniques  were 
not  applicable. 

The  Engineering  Experiment  Station  of  The  Georgia 
Institute  of  Technology,  acting  as  a  subcontractor  to  SAI  for 
the  electrical  properties  measurement  tasks,  suggested  that  a 
microwave  resonant  cavity  approach  could  prove  useful.  EES 
further  suggested  the  use  of  a  TMq^q  excitation  mode.  Although 
cavity  techniques  have  been  previously  used  for  the  measurement 
of  electrical  properties,  no  documented  investigations  could  be 
found  in  the  technical  literature  which  exactly  addressed  the 
problem  of  measuring  thin  dielectrics  mounted  on  metal 
facing-plates.  Furthermore,  the  method  had  to  be  reliable  and 
inexpensive.  For  these  reasons,  an  analytical  investigation  of 
the  cylindrical  cavity  measurement  procedure  was  conducted.  In 
this  Appendix,  we  present  the  results  of  that  analysis. 

The  conclusions  can  be  briefly  summarized  as: 

1.  A  TM  cylindrical  cavity  measurement  procedure  is  a 
reliable,  cost-effective  approach. 

2.  The  theoretical  precision  of  the  dielectric  constant 
and  loss  tangent  measurements  for  the  plasma-sprayed 
ceramic  materials  are  as  good  with  the  TM  cavity 
excitation  mode  as  they  are  with  the  more  common  TE 
excitation  mode. 

3.  The  construction  difficulties  associated  with  a  TM- 
excited  cavity  are  significantly  less  than  those 
associated  with  a  TE-excited  cavity. 
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In  the  remaining  portions  of  this  appendix,  we  examine  in  detail 
the  cavity  measurement  technique  for  the  determination  of  the 
electrical  properties  of  thin,  dielectric  materials  that  are 
mounted  on  metal  discs.  The  absolute  precision  of  the  TE  and 
TM  excitation  modes  is  investigated.  These  results  indicate 
that  the  suggestion  by  Georgia  Tech  of  the  ™q^q  m°de  for  the 
cavity  excitation  was  a  reasonable  choice.  Based  upon  their 
recommendation  as  well  as  the  analysis,  it  was  agreed  by  SAI  and 
Georgia  Institute  of  Technology  that  the  TMgj_o  resonant 
cavity  mode  would  be  used  in  the  measurement  procedure. 

A  sensitivity  analysis  was  performed  for  both  the 
dielectric  constant  extraction  process  and  the  loss  tangent 
extraction  process,  for  a  resonant  cavity  with  a  TM  excitation 
mode.  This  analysis  indicated  that  for  measurement  errors  at 
about  the  level  of  those  found  in  the  experimental  data  collection 
process,  a  true  dielectric  constant  of  seven  can  be  measured  with 
a  tolerance  of  +0.50  at  the  95%  confidence  level.  Similarly, 
loss  tangents  larger  than  0.002  can  be  determined  with  a  relative 
error  less  than  25%.  True  loss  tangent  values  less  than  0.002 
cannot  be  determined  with  an  error  bound  less  than  0.002.  The 
actual  loss  tangent  error  bounds  are  a  function  not  only  of  the 
nominal  loss  tangent  value,  but  also  of  the  nominal  dielectric 
constant  value.  The  actual  error  sensitivities  are  given  later 
in  graphical  form  for  some  representative  cases  of  interest. 

2 . 0  GEOMETRY 

The  geometry  of  this  problem  is  illustrated  by  Figure  C.l. 
A  right  circular  cylinder  of  radius  "a"  and  height  "d"  is 
partially  filled  from  the  bottom  to  a  height  "b"  with  a  material 
having  a  permittivity  e  =  ereo  and  a  conductivity  a.  The 
dielectric-filled  portion  of  the  cylinder  is  called  region  1, 
while  the  remaining  free  space  portion  is  called  region  2.  The 
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coordinate  system  is  oriented  such  that  the  z  -  axis 
cident  with  the  longitudinal  axis  of  the  cylinder. 


is  coin- 


Figure  C.l,  Circular  Resonator  With  Dielectric  Sample  Filling 
Region  1  (0  <  z  £  t>) 


/ 
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REVIEW  OF  CIRCULAR  RESONATOR  THEORY 


We  begin  with  Maxwell's  curl  equations 

VxH  =  J  +  yE  U 

V  x  E  =  -M  -  z  H  L2 


where 


y  =  jue  +  a  (3] 

A 

z  =  juy  (4] 

and  where  J  and  M  are  electric  and  magnetic  current  sources, 
respectively.  We  use  the  fact  that  all  boundaries  conform  to 
cylindrical  coordinate  surfaces  and  have  circular  symmetry  to 
conclude  that 


3 

J$ 


0 


(5) 


provided  that  any  excitation  is  ^-independent.  Expansion  of 
(1)  and  (2)  in  cylindrical  coordinates  and  application  of  (.5) 
results  in  two  independent  equation  sets,  one  TEz  and  one  TMz 
viz: 


I 


\lr 


3 

Jr 


32E. 


+  k2  E , 


z  J, 


4.0  TE  EXCITATION  MODES 

Since  the  two  mode  sets  do  not  cross-couple,  we  can  analyze 
each  set  independently.  We  shall  briefly  consider  the  TE  mode  set. 
This  approach  has  received  considerable  attention  in  the  literature 
(see  references  Cl  through  C5)  .  The  modes  which  are  excited  are 
found  by  solving  the  homogeneous  form  of  the  differential  equation 
in  (6)  to  give 


E4>  =  Ji(krr)  [A  sin  (kz2)  +  B  COS  (kzz)J 

H  ■  -r— —  J.  (k  r)  ["-A  sin  (k  z)  +  B  cos  (k  z)l 
r  1  r  L  z  z  -* 

=  1  3  (rE4>) 

”z  “ jwyr  3r 


(9) 

(10) 

(ID 


/ 
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where 


and 


(12) 


(13) 


Note  that  is  the  first  root  of  the  Bessel  function  J(kra)  ■  0; 
namely  =  3.832. 

This  mode  set  can  be  excited  by  J^,  and  Mz*  An  appli¬ 
cation  of  the  equivalence  principle 

n  x  H  =  J  (14) 

for  H  I  leads  to 

Z  | 

J.$  =  n  x 
<P 

-  Nr.a  1  1151 

Thus,  the  surface  current  density  J  on  the  wall  is  azimuthal  and 
excitation  coupling  could  be  obtained  with  an  axially  directed 
slot  in  the  cylindrical  cavity  wall. 


4.1  TE0in  Mode  Set 

The  TE  cavity  is  usually  operated  in  the  TEQ^n  mode.*  For 
this  mode,  the  cavity  is  above  cutoff  for  many  other  modes.  These 

*  For  circular  resonators,  the  modes  are  designated  TEk,m,n  and 
TMk,m,n*  Integer  k  determines  periodicity  in  <J>,  m  denotes  the 
number  of  electric-field  zeros  in  the  radial  direction  (not 
counting  zero  on  axis) ,  and  n  is  the  number  of  half -wavelengths 
in  the  axial  direction. 
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other  modes  will  obscure  the  electrical  properties  of  the 
material  under  test  if  they  are  allowed  to  exist  in  the  cavity. 
Figure  C  .2  provides  a  pictorial  view  of  the  modes  which  are 
physically  excitable  together  with  their  relative  cutoff  fre¬ 
quencies. 


TEn  TEqi  TEu  TEa i  TEit 

*  4  1  TV 

TMU  TMti  TMat 


Figure  C.2.  Relative  Cutoff  Frequencies  of 
Waves  in  a  Cylindrical  Resonator 

The  multimode  contamination  of  the  cavity  is  the  major  problem 
associated  with  TEQln  excitation.  Substantial  efforts  have 
been  described  in  the  literature  which  are  directed  at  minimizing 
the  excitation  of  unwanted  modes.  Unfortunately,  the  success 
of  these  approaches  is  directly  related  to  mechanical  tolerances 
which  have  finite  limits. 


4.2  Determination  of  the  Dielectric  Constant 

The  usual  perturbation  theory  approach  is  valid  for 
this  situation.  Consequently,  for  an  empty  cavity  resonant 
frequency  of  w0/2tt  and  a  dielectrically-filled  cavity  resonant 
frequency  of  co/2ir,  the  normalized  shift  in  the  resonant  fre¬ 
quency  of  the  filled  cavity  (u>-(Uq) /o)q  is 


X  A*  E,  av 
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where 


Ae  *  c  - 


(17) 


and  the  volume  V  is  composed  of  both  region  1  and  region  2.  The 
quantity  Eq  is  the  E— f ield  in  the  empty  cavity ,  while  E^n£  is 
the  E-field  in  the  sample.  For  samples  which  perturb  the  empty 
cavity  field  only  a  small  amount,  Eint  -  EQ  and  thus 


w-w. 


u> 


0 


(18) 


The  small  perturbation  assumption  also  allows  us  to 
approximate  the  permittivity  e  associated  with  the  energy  lost 
by  the  free  space  permittivity  eQ.  Then, 

-  -  k*H—s  &) 

“V  2J\h\  dv  v  0 1 

V 

where  41V  "  is  the  volume  of  the  material  sample  and  "V"  is  the 
s 

total  volume  of  the  empty  cavity.  But 


Af 


f-f0  _ 
f- 


(20) 
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Consequently,  we  are  left  with  the  job  of  computing  the  integral 
of  J Eq | ^  over  two  different  volumes.  From  equation  (9),  we 
know  that  EA  =  £,<£.  That  is,  the  unperturbed  E-field  in  the 

U  <y 

cavity  has  only  an  azimuthal  component,  viz: 


*  J1(krr)  B  sin(kzz)  <J> 


(21) 


and  as  before 


k  =  mr/d 
z 


2  2  2 
*  -  kr  +  kz 


Then, 


Af  . 


(Er-1) 


2ttB2/*  f  J?  (k  r)  sin2  (k  z)  r  dr  dz 

J  0J0  1  r  z _ 

3  a 


2ttB 2f  f  J2  (k 
J0J0  1  1 


2 

r)  sin  (k  z)  r  dr  dz 
2 


fhai  2 

/er-M  Jo  d 

"  V  2  I  -Ja  .  2  nlrz 


dz 


2  mrz  . 
sin  —3—  dz 


and  finally. 


4.3 


Sensitivity  in  the  Determination  of  the  Dielectric 
Constant  Value 

It  is  easily  seen  that  the  dielectric  constant  value  may 
be  found  from  the  normalized  resonant  frequency  shift  of  the 
cavity.  To  illustrate  the  precision  of  the  extraction  process, 
the  normalized  resonant  frequency  shift  is  plotted  as  a  function 
of  the  relative  dielectric  constant  in  Figure  C-3.  Variation 
of  the  cavity  fill  factor  b/d  produces  the  family  of  curves 
shown.  The  value  of  n  has  been  set  to  1  for  simplicity.  Hence, 
the  mode  of  operation  is  TEq^.  Note  that  the  curves  are  linear 
with  a  slope  s,  viz: 

s  =  ^  £l  -  sine  (2nit  b/d)j  (23) 

4 . 4  Fill  Factor 

In  order  to  interpret  Figure  C*3  correctly,  one  must 
know  the  fill  factor,  which  is  a  function  of  the  sample  thick¬ 
ness  b,  as  well  as  the  cavity  height  d.  As  an  example,  we  con¬ 
sider  two  representative  values  of  b;  these  are  0.050  inches 
and  0.125  inches.  The  cavity  height  can  be  found  from 


or 


^2 

(D  '  U 


oeo 


■  H^)2*  m 


(24) 
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Solving  for  d,  one  finds  that 


(.25) 


where  f  is  the  desired  resonant  frequency  in  Hertz  and  a  is  the 
radius  of  the  cavity.  Figure  C*4  illustrates  the  b/d  ratio  for 
b  =  0.050  inches  as  a  function  of  the  cavity  radius  in  inches. 

Three  frequencies  are  shown;  (1)  3.0  GHz  (S-Band) ,  (2)  5.0  GHz 
(C-Band)  and  (3)  10  GHz  (X-Band) .  Note  that  once  the  cavity 
radius  becomes  larger  than  about  four  inches,  further  increases 
do  not  significantly  affect  the  fill  factor.  We  continue  to 
restrict  our  attention  to  the  teqh  mode.  Figure  C.5  provides 
the  same  information  as  Figure  C.4,  except  that  in  Figure  C.5 
the  thickness  of  the  material  has  been  increased  to  0.125  inches. 

The  fill  factor  has  been  increased  by  the  larger  material 
thickness.  Note,  however,  that  a  cavity  radius  of  four  inches 
is  still  approximately  the  value  for  which  increasing  the  cavity 
radius  provides  little  additional  increase  in  the  fill  factor. 

To  briefly  summarize,  we  have  found  that  the  largest 
fill  factor  possible  for  a  0.050  inch  thick  sample  is  about 
0.08.  For  a  sample  that  is  about  0.125  inches  thick,  the 
greatest  fill  factor  which  can  be  achieved  is  about  0.20. 

Returning  to  Figure  C.3,  it  can  be  observed  from  the  non-zero 
slope  that  these  fill  factors  provide  some  sensitivity  to  the 
value  of  the  dielectric  constant.  However,  the  greatest 
sensitivity  occurs  at  the  higher  fill  factor  values.  Note 
that  the  cavity  configuration  under  consideration  has  a  large  radius 
but  a  very  small  height.  The  machining  of  the  axial  slot  for 
excitation  may  be  quite  difficult  due  to  the  small  cavity  height. 
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CAVITY  b  -  0.050 


$ 


CAVITY  RADIUS  (a)  [in.] 

C.4.  Fill  Factor  as  a  Function  of  Cavity  Radius  for  TE, 
Mode  and  Material  Thickness  of  0.050  Inches 
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In  the  next  sections,  we  will  determine  the  relationships 
for  the  TMQln  mode  corresponding  to  those  we  have  just  presented 
for  the  TEq11  mode.  Thus,  comparisons  between  the  two  excita¬ 
tion  modes  may  be  made  and  the  appropriate  one  selected  for 
further  investigation. 

5.0  TM  EXCITATION  MODES 

The  TM  excitation  modes  will  be  considered  in  the  same 
manner  as  the  TE  mode  set.  That  is,  we  shall  investigate  the 
cavity  fields  and  determine  the  effects  of  placing  a  material 
sample  at  one  end  of  the  cavity.  Although  the  analysis  is 
straightforward,  it  is  more  complicated  than  that  for  the  TE 
excitation  mode.  This  is  due  to  the  fact  that  the  perturbation 
approach  is  no  longer  valid,  and  a  more  complete  treatment  must 
be  undertaken. 

We  begin  by  referring  to  the  TM  set  of  equations  (7) . 

This  set  can  be  excited  by  sources  M , ,  J  ,  or  J  ,  provided  that 

**  $  i  z 

circular  symmetry  is  preserved.  Note  that  J  and  J  are  a 
continuous  current  that  flows  between  the  cylinder  walls  and 
the  endplates.  An  application  of  the  equivalence  principle 

-n  x  E  =  M  (26) 

at  the  endplate  at  z  =  d  results  in 


=  -M. 


z  =  d 


which  shows  that  an  annular  slot  in  the  endplate  is  a  possible 
excitation  method.  To  determine  the  modes  excited  we  may  solve 
the  homogeneous  form  of  the  differential  equation  in  (7)  to  give 


H*  -  Jj.  <v> 


cos  k  z 
z 


sin  k  z 
z 
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where 


*2  -  4 ♦  *z 


(31) 


Equations  (27)  through  (31)  form  the  basis  for  the  consideration 
of  both  the  unloaded  and  loaded  resonator . 


5.1 


Unloaded  Resonator 

In  the  unloaded  resonator,  the  boundary  conditions  are 


=  E 


z  =  0 


=  0 


(32) 


z  =  d 


=  0 


(33) 


r  =  a 


Application  of  these  conditions  in  (28)  through  (31)  results  in 


Hj>  =  A  Ji  (krr)  cos  kzz 


(34) 


E„  = 


k 

A  (krr)  sin  k^z 

y 


(35) 


L  -  A  Jn  (k  r)  cos  k  z 


(36) 


where 


a)2^oco 

(37) 

H55)2 

(38) 

/ 

fir 

(39) 

Consideration  of  the  resonance  conditions  in  expressions  (31) 
and  (37)  through  (39)  shows  that  the  dominant  mode  occurs  for 
n  =  0/  and  gives  Er  =  0  and 


H  =  A  j  (k  r) 
<p  1  r 


(40) 


E 

z 


(k  r) 
r 


(41) 


with  (31)  we  obtain 


~  =  2.405c 

r0  2  ira 


where  fg  is  the  resonant  frequency  of  the  dominant  mode  and 
c  is  the  velocity  of  light  in  free  space. 

The  Q  of  the  unloaded  resonator  can  be  obtained  from 
the  basic  definition 


Q  = 


coW 

P 


(43) 


where  W  is  the  energy  stored  in  the  resonator  and  P  is  the 

wer  dissipated  in  the  walls.  The  derivation  follows  classical 
and  will  not  be  repeated  here.  The  result  is 
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1.202  n, 


Q„  - 


U  Rs  (1  +  S> 


(44) 


where  is  the  unloaded  Q,  is  the  intrinsic  impedance  of 
free  space ,  and  Rg  is  the  resistivity  of  the  walls. 

5 . 2  Loaded  Resonator 

In  the  case  of  the  resonator  with  dielectric  loading 
from  z  =  0  to  z  =  b  (Figure  C.l),  Equations  (28)  through  (31) 
are  written  for  0  <  z  <  b  (region  1)  and  for  b  <  z  <  d  (Region  2) 
viz : 


V  =  C  J1  <krr)  cos  kzlz 


Lzl 


Jrl 


C  J1  (krr)  sin  kzlz 


0  <  z  <  b 


(45) 


E  i  =  —  C  JQ  (k  r)  cos  kzlz 

Yl 


%2  =  B  Ji  (krr)  COS  kz2  (d  "  Z) 


kz2 


'r2 


B  (krr)  sin  kz2  (d  -  z) 


b  <  z  <  d  (46) 


Ez2  =  ^  B  JQ  (krr)  cos  kz2  (d  -  z) 
y2 
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where 


and 


(47) 


(48) 


2 

w  y0 


(e 


2 

(1)  u 


oeo 


(49) 

(50) 


In  (47)  and  (48),  kr  is  given  by  (38)  but  kzl  and  kz2  must  be 
determined  in  conjunction  with  the  boundary  conditions  at 
z  =  b.  Equations  (45)  and  (46)  contain  the  boundary  conditions 
at  z  =  0  and  z  =  d.  Application  of  the  requirement  for  con¬ 
tinuity  of  the  tangential  E  and  H  fields  at  z  =  b  yields  the 
additional  relations 

C  cos  k  b  =  B  cos  k  -  (d  -  b)  (51) 

zl 


C  sin  kzlb 

*1 


M 

— —  B  sin  k  2  (d  -  b) 
~  z 

*2 


(52) 


tan  kzl  b  - - #-  tan  kz2  (d  -  b)  (53) 

yL  y2 

The  Q  of  the  loaded  resonator  can  be  obtained  from 
equation  (43),  viz: 
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<d(W,  +  W  ) 


where  and  W2  are  the  energy  stored  in  the  dielectric  and 

free  space  regions,  respectively;  also,  PD  is  the  power  loss 

in  the  dielectric  and  P„  is  the  power  loss  in  the  resonator 

w 

walls.  By  inverting  (54) ,  we  may  obtain 


i  =  _L  +  _L 

ql  qd  qw 


where  Q_  and  QTT  are  the  Q's  due  to  the  dielectric  loss  and 
D  W 

wall  loss,  respectively.  Explicitly, 


qd  - 


w(W1  +  W2) 


QW 


u)^  +  W2) 


where 


.b  .2tt  _a 


H,  -  \  f  f  f  si  |Er,|2+  |Ezl|2)  rdrdodz 


00  0 


W„  -  i  f*  r  f  (|E^I2  +  M2'  r  dr  d* dz  ( 


.b  2tt  a 


PD-i  f  f  f  -  (|Erll2  *  !Ez1l2) 


r  dr  d$  dz 


'Q  0  0 


J  TBS  »  ■  ■ 

w  2^ 


2  r  dr  d(j>  (61) 
z  =  d 


In  (58)  through  (61),  the  various  field  components  are  given 
by  (45)  and  (46). 

5.3  Determination  of  the  Dielectric  Constant 

Expression  (53)  can  be  solved  to  give  a  relationship 
among  the  resonant  frequency  f  of  the  loaded  resonator,  the 
unloaded  resonant  frequency  f q,  the  dielectric  constant  the 

loss  tangent  (tan  6)  and  the  resonator  dimensions.  Since  the 
equation  is  transcendental  and  has  complex  roots,  numerical 
methods  are  required  for  its  solution.  Fortunately,  for  the 
cases  of  present  interest  in  the  Georgia  Tech  experimental  pro¬ 
gram,  the  cavity  dimensions  are  small  enough  so  that  the  tangents 
can  be  replaced  by  their  arguments  with  the  result 

A 

>4 b  =  -  ^  4  (d  -  b»  (62) 


where  from  (47)  through  (50)  and  (3), 


1c  er  (1  -  j tan  6) 
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(64) 


and  where 


~  *  £r  (1  -  jtan  6) 

y2 


(65) 


k2  =  o)2y0e0  (66) 

tan  6  =  -2-  (67) 

U)£ 


Substitution  of  (63)  through  (67)  into  (62)  yields  a  complex 
algebraic  equation,  the  real  part  of  which  produces 


or 


(68) 


(69) 


To  compare  these  results  with  those  from  perturbation  theory, 
we  subtract  unity  from  both  sides  to  give 


(.70) 


/ 
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where 


Af  -  f  -  f  (71) 

o 

If  the  second  term  in  brackets  in  (70)  is  small  enough,  a 
Taylor  series  expansion  gives  the  approximation 

Af  _  1  (er  1  1  b  (72) 

(8) 

which  is  the  perturbation  theory  result  ,  with  quasi-static 
modification. 

5 . 4  Determination  of  the  Loss  Tangent 

The  loss  tangent  can  be  related  to  the  Q  of  the  dielectric 
by  consideration  of  Equation  (56) .  Substitution  of  (45)  and 
(46)  into  (58)  through  (60)  followed  by  substitution  of  the 
result  into  (56)  yields 

k^2  sin2  kz2  (d  -  z)  +  k2  cos2  kz2  (d  -  z) 

2  2 
k2,  sin  k2)zj  +  kp  |  cos  k^z  |  )  dz 

(73) 


Again,  we  shall  assume  that  the  cavity  dimensions  are  small 
enough  so  that  the  trigonometric  functions  in  (73)  may  be 
replaced  by  the  first  terms  in  their  Taylor  series  expansions. 
In  this  case,  the  dominant  contribution  from  both  integrals  is 


/ 


i 
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the  leading  (unity)  term  in  the  cosine  and  we  obtain 

°d  dhr  1  + 


where  we  have  used 


^|2  =  er  [X  +  (tan  6)2]*£r 


cos  kz^  b 


cos  k 


Since  Qq  is  not  an  observable,  we  also  require  a  calculation  of 
Qw  in  (57).  Substitution  of  (45)  and  (46)  into  (58),  (59), 
and  (61)  and  subsequent  substitution  of  that  result  into  (57) 
gives 


'W  ~  R 


^  /(|k*' s1"  k*,z|  +  k*  |cos  kz,2|  )  dz 
c|2  fa  +  2  f  J.os  kz)zjZ  dzj  +  |b|2  a  +  zj 


cos  kz2 


2 

1-1  y  [k22  sin2  kz2  (d  -  z)  +  k2  cos2  krf  (d  -  z)J  dz 


r2 1  b 


r  2  r  r 

2  a  +  2  f  I  cos  kz^z  dz  +  IbJ2  a  +  2  / 

0  J  '  L  b 


cos*  kz2  ^  “  z) 
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If  the  same  approximations  are  made  in  (77)  as  were  made  in  (73)  , 
the  result  is  simply 


°w  -  °u  (i) 

where  Qu  is  given  by  (44).  Therefore,  we  have  in  Equation  (55) 


M3 


where  is  an  observable  and  QD  and  Qu  are  given  in  (68)  and 
(44) ,  respectively. 

5.5  Sensitivity  in  the  Determination  of  the  Dielectric 

Constant 

An  indication  of  the  sensitivity  in  the  dielectric  con¬ 
stant  determination  is  given  by  Equation  (70) .  We  note  that 
the  b/d  ratio  is  set  by  the  resonator  design  and  that  the 
normalized  frequency  change  Af/fQ  is  an  observable  in  the 
experiment.  We  plot  normalized  frequency  change  versus  dielectric 
constant  with  the  b/d  ratio  as  a  parameter  (Figure  C.6) .  Note 
that  any  errors  in  the  measurement  of  Af/fQ,  while  relatively 
unimportant  in  the  determination  of  low  dielectric  constants, 
become  increasingly  significant  as  the  dielectric  constant 
increases.  Also,  the  smaller  the  b/d  ratio,  the  more  difficult 
the  task. 

5.6  Sensitivity  in  the  Determination  of  the  Loss  Tangent 

Equation  (79)  is  used  to  plot  Q^,  which  is  an 
observable  in  the  experiment,  as  a  function  of  the  logarithm  of 
the  loss  tangent,  with  dielectric  constant  as  a  parameter 
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FILL  FACTOR 


DIELECTRIC  CONSTANT 

Figure  C. 6,  Normalized  Frequency  Shift  vs.  Dielectric  Constant 


(Figure  C.7,  Figure  C.8,  and  Figure  C.9).  In  all  three  figures, 

the  b/d  ratio  is  0.25.  In  Figure  C.7,  for  a  conductivity  equal 

to  that  of  brass  (a,.  =  1.57  x  1 07 mhos/meter ) ,  we  note  that 

brass 

accurate  determination  of  the  loss  tangent  becomes  difficult  for 
tan  6  <  10"3.  This  limitation  becomes  more  stringent  in 
Figure  C.8  and  Figure  C.9,  which  give  the  results  for  conduc¬ 
tivities  one-half  and  one-quarter  that  of  brass,  respectively. 

5.7  Analysis  of  the  Approximations  in  a  Typical  Case 

Consider  the  Georgia  Tech  resonator,  viz: 
a  =  1.5" 
d  =  0.2" 
b  =  0.050” 
b/d  =  0.250 

For  this  case,  expression  (42)  yields 

f  =3.014  GHz 
o 

For  a  dielectric  constant  of  10,  Equation  (69)  gives 

j-  =  0.8803 
o 

and  thus 


we  obtain 


f  =  2.653 


kzjL  b  =  0.2083 


kz2  (d  “  b)  =  j  0.1141 


/ 
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Figure  C.9.  Loaded  Q  Versus  Loss  Tangent,  o  =  0.25  Brass 
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Although  these  results  are  obtained  from  the  approximation  in 
(62) ,  they  may  be  substituted  into  the  exact  Equation  (53) 
with  consistency  to  the  third  decimal  place  (0.3466  on  the 
left  side  of  equation;  0.3403  on  the  right  side).  Furthermore, 
we  may  calculate  the  B/C  ratio  in  Equation  (51)  and  obtain 

|  «  0.9720 

which  compares  well  with  our  approximation  of  unity.  The 
remaining  approximations  involve  the  expression  1  -  j  tan  6 
and  are  valid  as  long  as  tan  6  <<  1. 

It  should  be  pointed  out  that  the  relative  simplicity 
of  the  results  of  this  analysis  depends  on  the  trigonometric 
approximations.  For  cavity  sizes  where  the  approximations  are 
no  longer  valid.  Equations  (53) ,  (73) ,  and  (77)  must  be  used 
instead  of  (61),  (74),  and  (78).  Although  the  integrals  are 
easily  evaluated,  the  algebra  becomes  quite  unwieldy.  As  the 
resonator  is  presently  constructed,  the  field  patterns  are 
altered  only  slightly  by  addition  of  the  dielectric  sample 
except  in  one  readily  recognizable  fashion.  First,  the  additional 
Er  component  is  small.  Second,  because  of  the  small  "d"  dimen¬ 
sion,  the  additional  cosine  variation  in  z  in  the  and  E z 
components  is  small;  that  is 

cos  k^z  a:  1  (0  <  z  <_  b) 
cos  kz2  (d  =  z)  =r  1  (b  <_  z  <  d) 

This  fact  is  also  reflected  in  the  near-unity  value  of  the  B/C 
ratio.  The  one  significant  change  is  the  jump  (nearly  equal  to 
the  dielectric  constant)  in  the  value  of  normal  E  at  the  boundary 
z  =  b.  But  this  is  precisely  the  change  postulated  in  the 

(C8 ) 

quasi-static  modification  to  the  perturbation  theory  results 


-178- 


This  section  provides  an  analysis  of  the  sensitivity  of 
the  electrical  properties  test  procedure  to  errors  in  the  esti¬ 
mates  of  the  dielectric  constant  value  and  the  loss  tangent  value. 
The  expressions  are  complex  functions  of  a  number  of  variables, 
among  them  the  dielectric  constant,  the  Q  of  the  cavity  due  to 
wall  losses,  and  the  various  error  parameters. 

6.1  Errors  in  the  Dielectric  Constant 

Since  the  dielectric  constant  value  is  influenced  only 
by  the  resonant  frequency  shift,  the  simple  expression 


Aer  = 


2er2[l  -  b/d  (1  -  l/er)]1/2  A (fQ  -  f) 


has  been  found  which  relates  the  error  bounds  of  the  dielectric 
constant  to  the  error  in  the  determination  of  the  resonant  fre¬ 
quency  shift.  The  symbol  A  should  be  interpreted  as  "the  error 
in."  For  example,  if  the  error  in  the  normalized  resonant  fre¬ 
quency  shift  is  about  one  percent  [A(fg  -  f)/fQ  =  and 
dielectric  constant  is  found  to  be  8,  then  the  error  bound  for 
the  dielectric  constant  value  is  +0.50. 

Equation  (80)  was  obtained  by  differentiating  the  nor¬ 
malized  resonant  frequency  shift  of  Equation  (70)  with  respect 
to  the  dielectric  constant.  Some  algebraic  manipulation  leads 
directly  to  Equation  (80) . 


Analysis  of  the  95  percent  confidence  bounds  (3o)  on  the 
dielectric  constant  values  from  the  experiments  leads  to  a  con¬ 
clusion  that  the  corresponding  error  in  determining  the  resonant 


frequency  shift  is  about  +  3.0  MHz  for  an  unloaded  resonant 
frequency  of  3.0  GHz  [A(fg  -  f)/fg  =  Substitution  of  this 

value  into  Equation  (80)  for  the  various  dielectric  constant 
values  yields  Figure  C.10.  The  fill  factor  was  set  at  the 
typical  value  of  0.25  for  the  construction  of  the  figure. 

The  horizontal  axis  is  the  nominal  value  of  the  dielectric 
constant  of  the  material.  The  bounds  on  the  dielectric  constant 
value  are  given  along  the  vertical  axis.  The  nominal  value  trace 
is  the  center  of  the  three  lines.  The  upper  bound  at  the  95 
percent  confidence  level  on  the  value  of  the  dielectric  constant 
is  given  by  the  top  line.  The  lowest  of  the  three  lines  details 
the  lower  bound  on  the  dielectric  constant  value  at  the  95  per¬ 
cent  confidence  level.  In  order  to  fix  ideas,  a  material  with 
a  true  (nominal)  dielectric  constant  of  12,  for  example,  will 
have  its  value  experimentally  determined  to  be  lower  than  13  and 
greater  than  11  with  a  high  probability.  The  error  bounds  for 
other  nominal  values  of  the  dielectric  constant  can  be  found 
from  Figure  C-.10  or  from  the  use  of  Equation  (80)  . 

6.2  Errors  in  the  Loss  Tangent 

The  expression  for  the  loss  tangent  is  a  function  of 
both  the  measured  cavity  Q  and  the  dielectric  constant  of  the 
material.  Hence,  a  simple  expression  cannot  be  found  which 
relates  changes  in  both  factors  simultaneously  to  the  change  in 
the  loss  tangent  value.  The  effects  are  determined  individually 
and  compared  in  a  series  of  examples .  The  errors  in  the  measure¬ 
ment  of  the  cavity  Q  will  be  shown  to  dominate  the  loss  tangent 
error  function  at  least  for  the  size  of  the  error  in  the  dielec¬ 
tric  constant  value  found  in  the  experimental  work. 

When  the  analysis  results  are  applied  to  a  sample  case 
with  parameters  on  the  order  of  those  expected  for  the  plasma- 
sprayed  ceramic  materials,  the  minimum  error  bound  for  the  loss 
tangent  value  is  about  1  x  10  ^  to  2  x  10  Thus,  values  of 
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the  loss  tangent  whose  magnitude  is  less  than  this  figure  will 
have  confidence  bounds  that  are  so  great  that  the  loss  tangent 
value  will  be  meaningless.  A  complete  treatment  is  provided  in 
the  following  pages. 

6.2.1  Loss  Tangent  Errors  Due  to  Errors  in  Determining  the 
Dielectric  Constant 

The  Q  of  the  loaded  cavity,  Q  ,  can  be  related  to  the  loss 

Li 

tangent  and  dielectric  constant  of  the  test  material  by  sub¬ 
stituting  the  results  from  Equations  (56)  and  (57)  into  Equation 
( 55) ,  viz: 


i—  =  tan  6 


R  [1  +  a/d] 


|l_+',r(L^Zaj]+ 1-202  (120,)-  Cl  _  b/d  (1  _  1/£r)] 


17? 

(81) 


Rearranging  and  solving  for  the  loss  tangent  produces 


tan  5 


Rs  »  -  a/di 

Ql  1.202  (120,) 


[1 


1  +  e. 


-  b/d  (1 


(82) 


After  taking  the  partial  derivative  of  Equation  (82)  with 
respect  to  the  dielectric  constant  and  performing  some  algebraic 
manipulation,  the  following  form  can  be  generated,  viz: 


3  tan  6  _  (1  -  b/d) 

~TTr - (b/S~ 


tan  <5 


b) 


1  b/d 

N  ,  c  (i  -  b/a>l 

<2e  2[l  -  b/d  (1  -  1/e  )] 

L1  +  er  b7d J 

(83) 
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If  we  let  the  partial  derivatives  become  some  small  change. 


3  tan  6  A  tan  <5 


3  e_ 


A  e. 


then. 


*  tan  « 


tan  6  +  t 


(b/d) 


QW  2er2[l  -  b/d  (1  -  1/Er)} 


(84) 


where 


t» » - 1  +  «r  <85» 

If  the  value  for  the  error  in  the  dielectric  constant  is 
found  according  to  Equation  (80)  at  approximately  the  95  percent 
confidence  value  and  inserted  into  Equation  (85),  then  the  error 
in  determining  the  loss  tangent  due  to  errors  in  the  value  of 
the  dielectric  constant  is  depicted  by  Figure  C-ll.  The  fill 
factor  is  again  specified  as  the  0.25  value  used  in  Figure  C-10. 
The  nominal  loss  tangent  values  are  given  along  the  horizontal 
axis,  while  the  error  bounds  can  be  found  along  the  vertical  axis. 
The  nominal  values  of  the  dielectric  constant  (e  ■  8,  16,  24) 
are  shown.  The  upper  three  curves  (in  descending  order)  are  the 
upper  error  bounds  for  er  =  24,  er  *  16,  and  er  ■  8  respectively. 
The  lower  three  curves  (in  descending  order)  are  the  lower  bounds 
for  loss  tangent  errors  for  nominal  values  of  the  dielectric 
constant  equal  to  8,  16,  and  24,  respectively.  Note  that  this 
complete  analysis  has  assumed  that  there  are  no  errors  made  in 
determining  the  Q  of  the  loaded  cavity. 

From  Figure  C.ll  it  can  be  seen  that  the  error  bounds  on 
the  loss  tangent  value  become  tighter  as  the  dielectric  constant 
decreases.  The  minimum  detectable  loss  tangent  can  be  viewed  as 
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the  value  at  which  the  upper  bound  curves  intersect  the  vertical 

axis.  For  such  a  definition,  the  minimum  detectable  loss  tangent 

-4 

value  is  on  the  order  of  10  or  less.  Such  a  value  is  an 

extremely  small  figure  for  the  cavity  measurement  procedure. 

Unfortunately  we  shall  see  in  the  next  part  that  the  driver  in 

the  loss  tangent  error  is  an  error  in  determining  the  Q  of  the 

loaded  cavity.  When  this  factor  is  considered,  the  minimum 

detectable  loss  tangent  value  will  be  seen  to  be  much  greater 
-4 

than  the  10  value  due  to  dielectric  constant  errors. 


6.2.2  Loss  Tangent  Errors  Due  to  Errors  in  Finding  the  Q  of 
the  Loaded  Cavity 

We  begin  by  noting  that  Equation  (55)  can  be  rewritten 
to  produce  an  explicit  expression  for  as  a  function  of  Qw 
and  Qd,  viz: 


Q 


L 


QW  QP 
QW  +  °D 


(86) 


Taking  the  partial  derivative  of  Equation  (86)  with  respect  to 
the  loss  tangent  leads  to 


3  °L  -  ■  » 

T”tanT  UW  9  tan  6 


Qw  +  qd 


(87) 


Ww  T 

(Qw  +  Qd) 2  (tan  5) 2 


(88) 


where  x  is  given  by  Equation  (85) . 

If  the  relative  change  in  the  parameters  is  equated  to 
the  partial  differential  as 


8  .  *  °L 

9  tan  6  ~  A  tan  6 


(89) 
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then, 


i  tan  6  »  -  (1  +  (-nT5-2-  A  Q  (90) 

Qw  T  L 

The  parameters  A  tan  6  and  A  QL  can  be  interpreted  respectively 
as  the  error  in  the  loss  tangent  and  the  error  in  the  Q  of  the 
loaded  cavity.  Substitution  of  Equations  (56)  and  (57)  into 
Equation  (90)  above  leads  directly  to  the  result 


A  tan  6 


(tan  3)2  + 

T 


2  tan  5  + 
QW 


*  QL 


(91) 


For  the  unloaded  resonant  frequency  f q  and  the  3  db  band¬ 
width  of  the  cavity  resonance  f^^/  it  is  well  known  that 


(92) 


and  thus. 


d  Ql  = . 

d  f3db  f3db 


(93) 


The  usual  substitution  of  the  small  change  notation  for  the 
differential  produces  the  desired  result,  viz: 


A  °L 


A  f 


3db 


(94) 


Equation  (94)  specifies  the  relationship  between  the 
error  in  the  Q  of  the  loaded  cavity  and  the  error  in  measuring 
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the  3  db  bandwidth  of  the  resonant  response.  A  simple  rearrange 
ment 


A  °L 

Ql 


A  f 


3db 


'3db 


(95) 


shows  that  the  relative  error  of  the  two  parameters  are  identical. 
Because  the  determination  of  A  f exP®rimental  •  a 
reasonable  error  appears  to  be  on  the  order  of  5  to  6  percent. 

If  we  take  the  maximum  value  for  Q_  in  order  to  produce  a  worst- 

L 

case  estimate,  then  QL  -  Qw  =  850  and  A  *  50.  Note  that  the 
Qw  *  850  value  is  representative  of  the  experimental  results 
and  corresponds  to  a  cavity  constructed  from  a  metal  whose 
conductivity  is  equal  to  one-half  that  of  brass  ((Jj3rass  = 

1.57  x  10^  mhos/meter). 

The  major  result  previously  given  as  Equation  (91)  has 
been  plotted  in  Figure  C.12  for  the  three  values  of  the  dielec¬ 
tric  used  in  section  6.2.1.  The  error  in  the  Q  of  the  loaded 
cavity  is  assumed  to  be  50.  The  conductivity  of  the  metal  in 
cavity  is  assumed  to  be  one-half  that  of  brass. 

It  is  quite  apparent  that  the  minimum  detectable  loss 

tangent  value  due  to  this  representative  error  in  is  much 

greater  than  the  minimum  detectable  value  due  to  errors  in  the 

dielectric  constant.  In  particular,  it  appears  that  for  e  =  8 

-3  . 

the  minimum  detectable  loss  tangent  is  about  2  x  10  ,  while 
for  er  =  24  the  value  may  be  on  the  order  of  6  x  10  The 
criteria  for  the  minimum  detectable  loss  tangent  is  again  the 
intersection  of  the  vertical  axis  and  the  upper  bound  curves. 

Although  these  values  indicate  that  the  cavity  measure¬ 
ment  procedure  is  not  as  sensitive  as  some  other  measurement 
procedures,  the  cavity  approach  is  usable  for  very  thin  samples 
while  the  other  methods  generally  are  not.  In  addition,  since 
the  samples  are  so  thin,  the  calculated  losses  attributed  to 
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the  material  are  as  insignificant  for  tan  6  =  10  as  they  are 
for  tan  6  <  10~4.  Of  course,  this  is  only  true  provided  the 
frequency  is  low  enough  so  that  the  material  thickness  is  much, 
much  smaller  than  a  wavelength.  For  the  S-Band  and  X-Band 
frequencies  of  interest  here,  the  material  thickness  constraint 
is  always  satisfied. 
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1 . 0  INTRODUCTION 

In  order  to  provide  protection,  a  phased  array  antenna  is 
often  covered  with  a  dielectric  slab.  The  resulting  structure  has 
been  shown  to  support,  at  certain  angles  in  the  radiation  pat¬ 
tern,  a  surface  wave  mode  which  propagates  along  the  face  of  the 
antenna.  The  existence  of  this  mode  produces  a  resonance  which 
results  in  a  null  in  the  active  element  pattern.  This  phenomenon 
is  known  as  "array  blindness."  In  this  appendix,  we  follow  the 
analysis  of  Knittel,  Hessel,  and  Olinex^02^  to  produce  a  set  of 
curves  relating  slab  thickness,  dielectric  constant,  and  blind¬ 
ness  angle  for  a  simple  array  model. 

2 . 0  ARRAY  MODEL 

Consider  a  two  dimensional  slot  array  covered  by  a  uniform 
slab  of  dielectric  material.  To  analyze  the  surface  guided  waves, 
assume  that  the  slots  occupy  a  sufficiently  small  percentage  of 
the  total  area  of  the  antenna  face  that  the  dielectric  slab 
essentially  rests  on  a  uniform  ground  plane.  In  addition,  assume 
that  the  array  is  large  enough  that  the  great  majority  of  elements 
"see"  essentially  the  same  physical  environment.  There  is  there¬ 
fore  little  error  in  assuming  that  the  array  is  infinite  in  sur¬ 
face  area.  Finally  assume  that  the  array  aperture  distribution  is 
separable  so  that  the  scanning  phenomena  in  each  principal  plane 
can  be  separately  analyzed. 

3.0  DIELECTRIC  SLAB  ANALYSIS 

Consider  a  slab  of  thickness  T  and  dielectric  constant 

er  resting  on  a  perfectly  conducting  ground  plane  (Figure  D-l) . 

The  dominant  surface  guided  mode  supported  by  this  structure 

is  TM  .  We  may  express  the  x-dependence  of  the  fields  by 

isin  u£)  f  <  T  and  b  exp(-vx)  for  x  >  T.  In  addition,  the 

cos  ux  J  r 
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z— dependence  can  be  given  by  exp  (+ik  z)  where  the  relations  among 
u,  v,  and  k  are  given  for  the  TM  modes  by 

Z 

uT  tan  uT  =  er  vT  (2). 


u2  +  kz2  -  er  k2 


-v2  +  kz2  -  k2 


Elimination  of  v  and  k2  gives 


tan  uT  =  er  (er 


'(n)2  - 


2  2 

where  k  =  u  u  e  .  Equation  (.4 )  is  the  characteristic  equation 
o  o  (D5) 

for  TM  mode  propagation  in  the  slab.  It  is  not  hard  to  show 

that,  for  propagating  modes,  u,  v,  and  kz  are  all  real. 

4.0  EXCITATION  OF  TM  SLAB  MODES  BY  ARRAY 

The  TM  modes  in  the  slab  are  excited  by  the  slots  located 
at  a  uniform  spacing  d  and  polarized  with  the  electric  field 
vector  in  the  z-direction  (Figure  D*2) .  (Note  that  E-plane  scan 
produces  the  TM  slab  modes  whereas  H-plane  scan  would  produce  TE 
slab  modes) .  To  produce  a  main  beam  at  angle  9  requires  a 
progressive  phase  shift  a  between  radiators  of 

a  =»  kd  sin  9  (5) 

Let  k  .  be  the  z-directed  wave  number  along  the  antenna  face. 
zA 

Then, 
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k  .  *  -y  =  k  sin  0  (6) 

zA  a 

For  excitation  of  the  surface  modes,  kzd  and  kzAd  may  differ  only 
by  multiples  of  2mir,  where  m  is  an  integer,  or 

k  -  k  sin  e  +  251  (7) 

z  a 

Substitution  into  (2)  gives 

(i)2+  (-ne+^,2  =  cr  (8) 

which  is  -he  constraint  equation  imposed  on  the  slab  by  the  array 
in  the  direction  of  E-plane  scan. 

5.0  ARRAY  BLINDNESS 

The  combination  of  Equations  (.4)  and  (S)  give  the 
criteria  for  array  blindness.  Indeed,  substitution  of  (8)  into 
(4)  followed  by  some  rearrangement  yields 


(9) 


where 


2 

(sin  0  +  ttt) 


T  * 


(10) 


In  Figure  D.3/  we  plot  slab  thickness  as  a  function  of  /t~  with 
the  array  blindness  angle  as  a  parameter  for  the  lowest  mode 
(m  =  -1  in  Equation  (10)).  From  the  design  point  of  view, 
given  the  array  scanning  requirements  and  the  dielectric  constant 
of  the  slab  material,  the  curves  give  maximum  allowable  thickness 
of  the  slab.  Figure  D-.4  as  an  alternative  gives  the  slab  thick¬ 
ness  plotted  against  blindness  angle  with  dielectric  constant 
as  a  parameter. 
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(DIELECTRIC  CONSTANT) 

Slab  Thickness  vs  Dielectric  Constant 
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APPENDIX  E 


1.0  MECHANICAL  PROPERTIES  OF  HARDENED  RADAR  CANDIDATE 

MATERIALS 

An  extensive  search  of  the  mechanical  properties  of 
Hardened  Radar  Candidate  Materials  has  been  performed.  The 
available  data  has  been  reviewed  and  compiled  as  attached. 

2.0  MECHANICAL  PROPERTIES  OF  CANDIDATE  MATERIALS 

The  materials  under  consideration  are  grouped  into 
three  categories,  namely  nonmetallic  bulk  material,  metal  and 
plasma  sprayed  ceramics. 

Most  data  for  nonmetallic  bulk  materials  are  found 
from  references  El  through  E5.  Mechanical  properties  for  some 
bulk  materials  are  manufacturing  process  dependent;  however, 
limited  information  is  provided  by  the  suppliers.  Except  for 
columbium  the  data  for  metal  are  rather  complete.  The  main 
sources  are  references  El,  E4,  E6,  E7) .  A  number  Qf  technical  art¬ 
icles  funished  by  the  Metals  and  Ceramics  Information  Center 
provides  the  only  source  of  information  for  plasma  sprayed  cer¬ 
amics.  The  properties  of  plasma  sprayed  ceramics  strongly  de¬ 
pend  on  spray  conditions  (speed,  distance,  angle,  temperature) , 
particle  size  and  coating  thickness.  Because  so  many  variables 
are  involved,  mechanical  property  tests  are  required  to  complete 
and  re-evaluate  the  data  provided  herewith. 

The  mechanical  properties  in  the  survey  include  melt¬ 
ing  point,  specific  gravity,  Poisson's  ratio,  porosity,  linear 
thermal  expansion  coefficient,  modulus  of  elasticity,  tensile, 
compressive  and  bond  shear  strengths.  Temperature  independent 
properties  (the  first  four  items)  are  tabulated  as  shown  in 
Table  El.  Other  mechanical  properties  are  plotted  vs.  temper¬ 
ature  in  Figures  El  through  E4  for  nonmetallic  bulk  materials. 
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APPENDIX  F  -  THERMAL  ANALYSIS  RESULTS 


A  detailed  thermal  analysis  of  the  ceramic  coatings 
was  made  using  the  thermal  model  of  Section  4.1,  the  mater¬ 
ial  data  of  Appendix  A,  and  the  expected  heating  environment 
of  the  Thermal  Radiation  Simulator.  The  purpose  of  the  ana¬ 
lysis  was  to  size  the  material  specimens  for  the  tests,  to 
predict  the  materials  expected  performance  in  the  tests,  to 
provide  temperature  profiles  for  the  structural  analysis,  and 
to  parametrically  investigate  the  effect  of  various  material 
property  and  design  variables. 

The  expected  TRS  environments  used  are  showr  in  Fig¬ 
ure  FI.  The  radiation  only  environment  has  a  maximum  flux 
2 

of  350  cal/cm  -sec.  When  the  specimen  is  exposed  to  the  fire- 

2 

ball,  the  maximum  flux  used  is  1000  cal/cm  -sec.  Of  this, 

2 

350  cal/cm  -sec  is  due  to  radiation,  and  the  remainder  is 
assumed  to  be  due  to  a  combination  of  convection  and  conden¬ 
sation  of  the  gaseous  aluminum  oxide  on  the  samples.  Assum¬ 
ing  a  sample  absorptivity  of  0.2  the  maximum  flux  input  into 

2 

the  sample  is  then  0.2  X  350  =  70  cal/cm  -sec  radiation  and 
2 

650  cal/cm  -sec  convection/condensation.  The  fluence  follows 
these  same  relationships. 

The  radiation  only  heating  is  what  the  samples  pro¬ 
tected  by  the  pyrex  glass  should  receive,  while  the  total 
fireball  heating  is  the  maximum  that  the  samples  exposed  to 
the  fireball  could  receive.  It  was  felt  that  these  two  heat 
profiles  should  define  the  lower  and  upper  heating  boundaries 
for  sizing  the  samples  and  instrumentation.  Also,  the  fire¬ 
ball  heating  is  similar  to  the  actual  nuclear  environment  when 
the  material  has  an  absoptivity  of  0.2  (80%  of  the  nuclear 
heating  would  be  reflected.  This  similarity  is  shown  in 
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in  Figure  F2  in  which  the  heat  received  by  the  sample  with  an 
absorptivity  of  0.2  is  plotted  for  the  two  environments. 

The  thermal  analysis  was  performed  for  both  the  alum¬ 
ina  and  zironia  coatings,  with  and  without  the  silicone  resin 
infiltrant.  A  schematic  of  the  thermal  model  is  shown  in  Fig¬ 
ure  F3.  Parameters  investigated  include  coating  thickness, 
steel  substrate  thickness,  and  porosity.  The  effect  of  absorp¬ 
tivity  was  investigated  using  the  nuclear  environment  (pre¬ 
sented  in  Section  4. 1.1. 2)  and  using  the  solar  furnace  environ¬ 
ment  (presented  in  Section  4.4.3). 

The  maximum  coating  surface  temperature,  maximum  in¬ 
terface  temperature,  and  maximum  substrate  backface  temperatures 
calculated  are  presented  in  Table  A1  for  the  various  models  con¬ 
sidered.  All  cases  used  a  surface  absorptivity  and  emissivity 
of  0.2. 

Examples  of  the  in-depth  temperature  profiles  and  sur¬ 
face  recessions  calculated  are  presented  in  Figures  F4,  F5, 
and  F6  for  .06  inch  alumina,  .09  inch  alumina,  and  .06  inch 
zirconia  coatings  respectively,  exposed  to  the  fireball  en¬ 
vironment. 

As  indicated  in  Figures  F4-F6,  the  ceramic  coating  sur¬ 
face  heats  to  its  melt  temperature  immediately  and  the  surface 
begins  to  recede.  The  in-depth  temperature  profiles  show  the 
temperatures  in  the  material  as  a  function  of  time.  As  the 
surface  recedes,  the  origin  of  the  coordinate  system  follows 
the  receding  surface,  so  that  the  0.0  on  the  abscissa  is  always 
at  the  new  material  surface. 

The  alumina  and  zirconia  behave  quite  similarly  ther¬ 
mally,  except  the  zirconia  has  a  higher  melt  temperature. 
Comparing  cases  6  to  16  and  8  to  18,  the  recession  for  each 
is  similar.  The  steel  substrate  temperature  is  a  little 
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IN-DEPTH  TEMPERATURE  PROFILES 


Figure  F4 .  0,060  Inch  Alumina  on  0,1875  Inch  Stainless 
Steel  Exposed  to  TRS  Fireball 
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Figure  F6.  0.060  Inch  Zirconia  on  0,1875  Inch  Stainless  Steel 
Exposed  to  TRS  Fireball 
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higher  for  the  zirconia.  Even  though  the  zirconia  would  re¬ 
radiate  away  more  heat,  its  lower  thermal  mass  appears  to 
amount  for  this  higher  substrate  temperature. 

Considering  the  coating  thickness  effects,  comparing 
cases  6  to  8  and  16  to  18,  the  thicker  coating  exhibits  more 
ablation.  The  thinner  coating  is  less  of  an  insulator  and 
permits  more  of  the  heat  to  be  conducted  into  the  steel  sub¬ 
strate  during  the  ablation  and  the  melting  of  the  surface 
ceases  at  an  earlier  time  during  the  exposure.  Of  course, 
even  though  the  thicker  material  ablates  more,  it  still  has 
more  material  remaining  after  the  exposure.  Hence,  the  thicker 
coatings  may  be  necessary. 

Considering  the  steel  substrate  thickness,  the  steel 
acts  as  a  thermal  heat  sink  with  a  very  low  temperature  grad¬ 
ient.  The  thicknesses  considered,  .1875  and  0.25  inches  for 
the  ablating  cases,  have  sufficiently  large  thermal  masses 
that  the  more  resistive  ceramic  coating  behaves  independently 
of  the  steel  substrate.  This  can  be  seen  by  comparing  cases 
6  to  7  and  8  to  9 .  The  steel  temperature  varies  according  to 
its  thermal  mass. 

The  effect  of  changing  the  porosity  in  the  non-infil- 
trated  coatings  effectively  changes  the  conductivity  and  ther¬ 
mal  mass.  This  is  shown  by  comparing  case  2  to  case  5.  In 
case  5  the  higher  conductivity  of  bulk  material  and  its  higher 
mass  was  used.  Heat  is  conducted  through  the  coating  faster, 
and  the  coating  absorbs  more  heat  at  a  given  temperatre,  re¬ 
sulting  in  a  lower  coating  surface  temperature.  Case  5  has 
a  slightly  higher  steel  temperature,  due  to  less  heat  being 
reradiated  away  at  the  ceramic  surface. 

Next,  the  effect  of  the  silicone  resin  was  investi¬ 
gated.  The  results  are  summarixed  in  Table  F2  in  which  the 
case  numbers  correspond  to  the  non-inf iltrant  cases  discussed 
earlier. 


Comparison  of  Performance  With  and  Without  Res 


Considering  the  low  level  heating  of  the  TRS  radia¬ 
tion  only  environment,  the  temperatures  do  not  get  high  enough 
over  a  sufficient  length  of  time  to  cause  resin  decomposition. 
The  effects  on  temperature  are  negligible  and  only  due  to  the 
resin's  effect  on  conductivity  and  thermal  mass. 

At  the  high  level  heating  environment,  the  resin  decom 
poses  and  reduces  the  ceramic  ablation  and  gives  a  somwhat  re¬ 
duced  steel  substrate  temperature.  By  comparing  Case  8  to  6, 
it  can  be  seen  that  the  thicker  coating,  .090  inch  in  Case  8, 
which  thus  contains  more  infiltrant,  has  a  greater  reduction 
in  ablation.  However,  the  resin  effects  are  small  as  there 
is  still  significant  ablation.  The  resin  is  much  less  effec¬ 
tive  in  this  radiation  environment  than  in  a  convective  en¬ 
vironment,  such  as  high  speed  flight,  where  blowing  due  to  the 
decomposing  resin  aids  the  reduction  in  heat  transfer. 


■ 


APPENDIX  G  -  FRENCH  SOLAR  TEST  DATA 

This  appendix  contains  the  French  Solar  Furnace  Test 
data.  Table  Gl  presents  the  samples  pre-test  and  post-test 
dimensional  and  mass  data.  Figures  Gl  through  G6  present 
the  pre-test  and  post-test  photos  of  each  sample.  Figures 
G7  through  G12  present  the  measured  front  face  and  back  face 
temperatures.  The  instrumentation  and  temperatures  are 
discussed  in  Section  4.4. 


-227- 


Table  61.  Sample  Dimension  and  Weight  Data 
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Figure  Gl.  Alumina  Coating/Silicone  Infiltrant,  Sample  F2 
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Figure  G2.  0.25  Inch  Silastic  E,  Sample  F3 
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Figure  G3,  0,5  Inch  Silastic  E,  Sample  F4 
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POST-TEST 

Figure  G^.  Boron  Nitride  (HBR),  Sample  F8 
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Figure  G 5 .  AS-3DX>  Sample  F9 
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Figure  G6,  Magnesia-Spinel  Coating/Silicone  Infiltrant,  Sample  F16 
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Alumina  Coating/Silicone  Infiltrant 
Sample  F2 ,  Temperature  Response 
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Front  Face 


Figure  G9,  0,5  Inch  Silastic  E/0,125  Inch  Steel,  Sample  F4, 
Temperature  Response 


TIME 


Figure  G12,  Magnesia  Spinel  Coating/Silicone  Infiltran 
Sample  F16,  Temperatupe  Response 


